o 
o 

(N 

G 

S : Chiral Perturbation Theory on the 

Lattice and its Applications 

T 1 , 

^H ■ 

o 
\o 

Tt '• Daniel Arndt 

o 

June 2004 

6-;, 

> 



Abstract 

Chiral perturbation theory (xPT), the low-energy effective theory of QCD, can be used 
to describe QCD observables in the low-energy region in a model-independent way At 
any given order in the chiral expansion, %PT introduces a finite number of parame- 
ters that encode the short-distance physics and that must be determined from exper- 
iment or numerical lattice QCD simulations. In this thesis, we calculate a number of 
hadronic observables in the quenched and partially quenched versions of xPT: 

Chiral corrections to B^ — ► Z?W at zero recoil are investigated in quenched ;\PT. 
We study in detail the charge radii of the meson and baryon octets, electromagnetic 
properties of the baryon decuplet, and the baryon decuplet to octet electromagnetic 
transitions in both, quenched and partially quenched xPT. We further show how ef- 
fects due to the finite size of the lattice can be accounted for in heavy meson %PT and 
calculate, as explicit examples, neutral B meson mixing and the heavy-light meson 
decay constants. We also demonstrate how one can account for effects due to finite 
lattice spacing in the low-energy theories, considering as an example electromagnetic 
meson and baryon properties. 

The results of our calculations are crucial to extrapolate quenched and partially 
quenched lattice data from the heavier light quark masses used on the lattice to the 
physical values. 
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Glossary 



QCD: Quantum Chromodynamics. 
LQCD: Lattice QCD. 
QQCD: Quenched QCD. 
PQQCD: Partially QQCD. 
%PT: Chiral Perturbation Theory. 
QxPT: Quenched xPT. 
PQxPT: Partially QxPT. 
HMxPT: Heavy Meson xPT. 



EFT: Effective Field Theory. 

HQET: Heavy Quark Effective Theory. 

CKM: Cabibbo-Kobayashi-Maskawa. 

LEC: low-energy constant. 

LO: leading order. 

NLO: next-to-LO. 

NNLO: next-to-NLO. 
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Chapter 1 

Introduction 



Quantum chromodynamics (QCD), a part of the very successful standard model of 
particle physics, was formulated over 30 years ago. It is the theory that describes 
the interaction of quarks and gluons, which are the building blocks of hadrons. In 
principle, it not only enables the calculation of properties of protons and neutrons, that 
make up the nuclei of atoms, but of the nuclei themselves. In short, QCD describes all 
hadronic properties of matter. 

Unfortunately, even though QCD is simple enough to be written down in the form 
of a few partial differential equations, solving it to calculate even basic hadronic prop- 
erties, such as the mass or the magnetic moment of the proton, is very complicated 
and still poses a challenge. In a similar theory, quantum electrodynamics, the fun- 
damental coupling constant is small at low energies and observables can be arranged 
in the form of a series, called a perturbative expansion. Therefore, in calculating a 
property, such as the anomalous magnetic moment of the electron, one only needs to 
consider the first few terms in this series that dominate — usually a cumbersome but 
easy task — whereas subsequent terms can be neglected because they are small. 

In QCD the picture is different. The fundamental coupling constant in QCD, as, 
depends upon the energy exchanged in the process under consideration (see Fig. ll.lt in 
a different way: as is small at large energies, such as occur during a particle collision 
in a large particle accelerator or in a quark-gluon plasma. Here, perturbative tech- 
niques are applicable. At energies < 1 GeV, however, the coupling constant becomes 
large, as ~ 1, and a perturbative expansion in powers of as fails. The series does not 
converge, instead of becoming smaller, subsequent terms get bigger and bigger. This 
is what makes solving QCD so complicated in the low energy region which is relevant 
for hadronic properties because that is where quarks and gluons bind together into 
composite states. 

A way to solve this problem is to use lattice QCD (LQCD). Here, one simulates 
QCD with the help of computers on a finite-sized 4-dimensional grid (or lattice) that 
represents points in discretized space-time (see Fig. I1.2D . Since even for today's most 
powerful computers this task is very time-consuming, theorists use additional approx- 
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Figure 1.1: Measured as plotted against the momentum /j, at which the measurement was 
made. The three lines show the central value and the ±la of the Particle Data Group's average. 
The data points are from several experimental measurements. Figure taken from Ref. [1]. 



imations to simplify the calculation: among others, they neglect (or partially neglect) 
contributions of quark-antiquark pairs that constantly pop out of and disappear into 
the vacuum (the so-called "quenched" and "partially quenched" approximations) and 
that are very costly to calculate; and they simulate with light quarks (of the up and 
down flavors) that are several times more massive than in nature. 

Because of such approximations, lattice theorists need to know how to connect their 
results to QCD of the real world. In particular, for each property they measure on the 
lattice they need to know how to extrapolate from the heavier quarks they use down to 
the quark masses of nature. A model-independent way to do this extrapolation is to use 
a low-energy effective theory that exploits the symmetries of QCD and is formulated in 
terms of the relevant degrees of freedom in the low-energy region, mesons and baryons, 
rather than quarks and gluons: chiral perturbation theory (xPT). Since the quark 
mass, m q , dependence is explicit in xPT (the low-energy constants are independent 
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Figure 1.2: Lattice size L and lattice spacing a. In general, these can be different for each of 
the four dimensions. However, often only the time direction has different L and a. 



of m q ) it is the only rigorous tool for extrapolating LQCD results down to physical 
quark masses. Since most simulations today use either the quenched or the partially 
quenched approximations of QCD, one has to use the quenched or partially quenched 
versions of %PT to do the appropriate extrapolations. 

The research contained in this thesis involves the calculation of a number of 
hadronic properties using xPT as well as quenched %PT (QxPT) and partially 
quenched %PT (PQxPT). Our results for QxPT are necessary to extrapolate existing 
quenched lattice data of these properties to the physical regime. 1 Moreover, because 
of the conceptual advances in lattice computing algorithms made in the past few years 
and because of the availability of faster computers, many of these properties will be 
simulated with improved precision in partially quenched QCD, although it will be a 
long time before real simulations with light physical quarks become feasible. There- 
fore these lattice results need to be extrapolated to real-world QCD using our results 
for PQxPT 

Besides (partial) quenching and simulating at heavier light quarks, there is a num- 
ber of further artifacts that come about by using LQCD and that must be taken into 
account and included in the (P)QxPT treatment. As steps in this directions, we have 
included 0(a) effects (due to the non-zero lattice spacing, a, the "graininess" of the 
discrete lattice) in the calculation of baryon properties and finite L effects (due to the 
finite size of the lattice box, L) in the calculation of properties of heavy quark systems. 



1 Note that our results in the present form cannot be used to extrapolate staggered lattice simulations. 



This thesis contains work carried out over the past two years and it is laid out as 
follows. In Chapter El we give a brief introduction into QCD, xPT, LQCD, QxPT, and 
PQxPT that is needed for all subsequent chapters. We elaborate on the implications of 
the quenched and partially quenched theories for the extrapolation of lattice QCD sim- 
ulations carried out at an unphysical regime to the physical regime. The subsequent 
chapters deal with the calculation of a number of hadronic properties. Chapters |3] and 
HI involve heavy mesons. In Chapter |3j we study the semileptonic JBM — ► £)(*) decays 
in the heavy quark limit and calculate the lowest order chiral corrections from the 
breaking of heavy quark symmetry at the zero recoil point in QxPT [2]. In Chapter HI 
we incorporate finite volume effects in the calculation of properties of heavy quark sys- 
tems. In particular, we investigate how the scale A, which comes from the breaking of 
heavy quark symmetry, influences finite volume effects. This work was carried out in 
collaboration with David Lin [3]. In Chapters HHZl we calculate a number of hadronic 
properties in the baryon sector in QxPT and PQxPT. Whereas Chapter |5] involves the 
baryon octet, Chapter |6] and Chapter |7] deal with the baryon decuplet and the baryonic 
octet-decuplet transition, respectively. In Chapter HO we extend the calculation of the 
subsequent three chapters by incorporating finite a effects. These chapters are work 
done in collaboration with Brian Tiburzi [4-7]. Finally, in Chapter |9] we summarize 
and conclude. Several appendices contain supplemental material that has been taken 
out of the main text in order to improve readability. 

Most of the work contained in this thesis has been published previously: 

• Daniel Arndt, Chiral 1/Mq Corrections to B^*> — > JDW at Zero Recoil in Quenched 
Chiral Perturbation Theory, Phys. Rev. D 67, 074501 (2003). 

• Daniel Arndt and Brian C. Tiburzi, Charge Radii of the Meson and Baryon Octets 
in Quenched and Partially Quenched Chiral Perturbation Theory, Phys. Rev. D 
68, 094501 (2003). 

• Daniel Arndt and Brian C. Tiburzi, Electromagnetic Properties of the Baryon De- 
cuplet in Quenched and Partially Quenched Chiral Perturbation Theory, Phys. 
Rev. D 68, 114503 (2003), Erratum-ibid. D 69, 059904 (2004). 



• 



• 



Daniel Arndt and Brian C. Tiburzi, Baryon Decuplet to Octet Electromagnetic 
Transitions in Quenched and Partially Quenched Chiral Perturbation Theory, 
Phys. Rev. D 69, 014501 (2004). 

Daniel Arndt and Brian C. Tiburzi, Hadronic Electromagnetic Properties at Finite 
Lattice Spacing, Phys. Rev. D 69, 114503 (2004). 



• Daniel Arndt and C.J. David Lin, Heavy Meson Chiral Perturbation Theory in 
Finite Volume, Phys. Rev. D in press, [hep-lat/0403012]. 



Chapter 2 

QCD, Chiral Perturbation 
Theory, and the Lattice 



In this chapter we introduce the part of the standard model that describes the in- 
teractions of quarks and gluons, QCD. Since QCD can only be solved perturbativly 
at energies well above the energy scale relevant for hadronic properties, we describe 
;\PT, which is QCD's low-energy effective theory. %PT — as an effective field the- 
ory (EFT) — introduces a number of unknown parameters that encode the underlying 
short-distance physics and that have to be fixed by comparison to either experimen- 
tal measurements or to results from numerical lattice QCD (LQCD) simulations. We 
briefly describe LQCD and the quenched and partially quenched approximations that 
are used frequently and introduce the low-energy chiral effective theories that can be 
used to extrapolate results from lattice simulations that employ these approximations: 
QxPT and PQ%PT. Lastly, we comment on how reliable it is to predict QCD properties 
from lattice simulations that use the quenched and partially quenched approxima- 
tions. 

2.1 QCD and Chiral Symmetries 

The Lagrangian of QCD is given by 

£qcd = ~\g%G a ^+ Yl Va(HP-rn q ) ab q b (2.1) 

a,b=u,d,s 

where the eight gauge bosons A A are contained in the gluon field strength tensor G A V 
that is given by 

G A v = d»A A -d v A A -gf ABC A B A c v . (2.2) 

The structure constants f ABC are defined by 

[T A ,T B ]=if ABC T c (2.3) 



and the Ta are the eight generators of color SU(3). The Lagrangian also contains the 
triplet of quarks q = (u, d, s) of the up, down, and strange flavors 1 with mass matrix 

m q = diag(m u ,m d ,m s ). (2.4) 

The quarks are minimally coupled to the gluon fields via 

D^ = d^ + igA*T A . (2.5) 

Assuming that one is in a regime where perturbation theory is applicable, one can 
calculate how the strong coupling constant 

asW = —. — (2.6) 

depends on the renormalization scale /i. From the QCD /3 function calculated to 0(g 3 ) 
one finds 

127T 

as ^ = (33-2N f) io gi »yA qcB iy (2 - 7) 

This means that, as long as the number of quark flavors Nf is smaller than 16, as 
becomes larger with decreasing \i, a behavior known as asymptotic freedom. Moreover, 
if n —* Aq CD then as blows up. Of course, in that case the theory is not perturbative 
in the first place. However, Aqc D , which can be determined from fitting Eq. ( 12 .71 1 to 
experimental measurements to be about 200 MeV (see Fig. Il.lt . can still be viewed as 
the scale where QCD becomes strongly coupled. 

In the limit of vanishing quark masses (m q — » 0) the quark part of Eq. d2.lt becomes 
simply 

£ = qilfa = q L ilf)q L + q R ilpq R (2.8) 

which exhibits an exact global SU(3)l x SU(3)r symmetry. This means that the left- 
and right-handed quark fields, 

1 - 7 5 1 + 7 5 

qh = — g — q and qR = — 2 — g ' (2 ' 9) 

transform under independent SU(3) flavor space rotations, 

q L ->■ Lq L , q R ->• Rq R , (2.10) 

with L G SU{3) L and R £ SU(3) R . 

In nature, the masses of the light quarks are not zero. If they are turned on then 
the term 

qm q q = q L m q q R + q R m q q L (2.11) 

1 Although the standard model has six flavors of quarks that, in principle, should all be included, only 
the three lightest flavors are relevant for the calculation of hadronic properties at energies < 1 GeV. The 
quarks of the charm, top, and bottom flavors have masses that are typically much larger than 1 GeV. 



appears in the Lagrangian which is only invariant if L = R. In that case the symmetry 
is broken down to its diagonal subgroup: SU{3)l x SU(3)r — » 5C/(3)y. Although the 
masses of the light quarks are not zero, they are nevertheless small compared to the 
scale Aq CD . One would therefore expect nature to exhibit at least an approximate 
SU(3)l x SU(3)r symmetry. However, such a symmetry is not seen. What one does 
see is evidence for just a single SU(3). One therefore assumes that the symmetry 
SU(3)l x SU(3)r is spontaneously broken down to the observed SU(3). This symmetry 
breaking is accomplished by the formation of scalar quark bilinears qq that have a 
non-zero vacuum expectation value: 

(q R q{) = \&. (2.12) 

Under an SU(3)l x SU(3)r transformation this vacuum expectation value becomes 

(&<&) - X(LR^, (2.13) 

which means that for A / the vacuum expectation value is only unchanged if L = R 
and the chiral symmetry SU(3)l x SU(3)r is spontaneously broken down to its diagonal 
subgroup SU(3)v- The eight broken subgroups cause the appearance of eight massless 
Goldstone bosons that are the fluctuations along the directions where the potential is 
constant. These eight Goldstone bosons are believed to be realized in nature as the 
pseudoscalar meson octet. The fact that the pseudoscalar mesons are light but not 
massless reflects the fact that the SU(3)l x SU(3)r is only an approximate symmetry 
of the Lagrangian. The Goldstone bosons can be represented by a 3 x 3 matrix £ that 
transforms as 

£ -► L£i? f (2.14) 

and can be written as 

S = exp^V (2.15) 

where $ is a traceless hermitian matrix given by 

* = I *- -^7T° + -^n K» I (2.16) 

K~ K° 

and / is a constant with dimensions of mass, known as the pion decay constant. 

At energies below A x , 2 these Goldstone bosons are the only degrees of freedom 
and one can write down an effective Lagrangian that describes their interactions. In 
principle, any term with the correct dimensions that obeys all the symmetries of the 
QCD Lagrangian in Eq. ( 12.11 1 can be included in such an effective Lagrangian. How- 
ever, since the number of such terms is infinite, one has to use a truncation scheme 

z Actually, energies should be below the mass of the p, m p — 770 MeV, since the p is not included in the 
EFT. This is accomplished by treating m p ~ A x . 




that limits the number of terms to be included: Because higher order terms contain 
more and more derivatives, they are suppressed by the scale A x . For the theory with 
massless quarks, the lowest order term is 



f-tr (d^duE) . (2.17) 



If m q is non-zero, then the term qm q q in Eq. 1 12.11 ) is not invariant under SU(3)l x 
SU(3)r. This can be fixed by treating m q as an independent field, a so-called spurion, 
that is assumed to transform as 

m q -► Lm q R j . (2.18) 

Under simultaneous chiral transformations of the quark and spurion fields the mass 
term in Eq. ( 12.1ft is invariant. Using the spurion technique one can then include terms 
into the EFT Lagrangian that involve m q . Doing so yields the lowest order Lagrangian 
of xPT that is of order p 2 /A 2 x [8, 9] 



C = 4-tr (d^djl) + Atr fm S + mls t 



+ . . . . (2.19) 

This Lagrangian includes all possible terms up to order p 2 /A 2 ,; terms that are of higher 
order in the p/A x chiral expansion (more derivatives, more powers of m q ) have been 
neglected. 

The Lagrangian in Eq. ( 12.19ft is valid as long as p <C A x and p/A x can be used as a 
small expansion parameters. It can be systematically expanded to include terms that 
are of higher order in p/A x . Each term is accompanied by a constant (/ and A in the 
above lowest order Lagrangian) that is a priori unknown. Observables receive con- 
tributions from both long-range and short range physics; the long-range contribution 
arises from the (non-analytic) structure of pion loop contributions, while the short- 
range contribution is encoded in these low-energy constants that appear in the chiral 
Lagrangian and are unconstrained in xPT. These constants must be determined from 
experiment or lattice simulations. 

By expanding Eq. ( 12. 19ft to lowest order in the meson fields <5 one can calculate the 
masses of the Goldstone bosons in terms of the quark masses, m q , and the constants / 
and A. For the masses of off-diagonal meson, that are made up of the (anti-)quarks q 
and q', one finds 

(m q + m' q ) . (2.20) 



2 _4A 

m qq> - p 



For example, this gives 



4A . . 4A 

-2 (m u + m d ) , m K ± = -^ 



m 7T± '- z 72 ( m u + m d ) , m K ± = -72 (m u + m s ) , (2.21) 



and 



4A 
m K o = m^o = -jz (m d + m s ) . (2.22) 



Clearly, in the isospin limit, where m u = ma, the charged and neutral kaons have equal 
mass. Similarly one finds for the masses of the mesons on the diagonal 



4A . . 4A 

-j2 ym u + m d ) and m v = — ~ 2 



m TT° = : 72 ( m « + md ^ and m v = TJ2 ( m « + m d + 4m s ) , (2.23) 



so that m n o = m n ± in the isospin limit. In nature, the kaons are much heavier than 
the pions, which reflects the fact that m s » m u , m d . 

2.2 Lattice QCD 

If the coupling as is small then one can use perturbative methods to calculate the 
vacuum expectation value of an operator O from the path integral 3 

(O) = 7F / VA V ^ V 1 Oexpf-f d A x CJ (2.24) 

with the generating functional Z defined as 

VAVqVqexp(- f d 4 xC) (2.25) 



by expanding the exponential in powers of the interacting part of the Lagrangian and 
solving the functional integral analytically for the first few terms in the series. This 
approach fails in the strong coupling region (for energies smaller < Aq CD ) because the 
expansion parameter as becomes large. A way to solve QCD in the strong coupling re- 
gion has been proposed by Wilson [10] and it involves putting QCD on a 4-dimensional 
discrete space-time lattice and solving it numerically using computers. This method, 
known as lattice QCD, basically involves two steps: 

1. The infinite dimensional functional integral in Eq. ( 12.241 1 needs to be discretized 
so that it can be calculated in a finite number of steps. This is accomplished by 
discretizing space-time and putting QCD in a 4-dimensional space-time lattice. 
In Wilsons formulation of LQCD, the fermionic fields (quarks) live on the lattice 
sites whereas the gauge fields (gluons) are defined on the links which are the 
lines that connect neighboring lattice sites. 

2. Even after this discretization solving the functional integral means summing 
over an enormous number of paths in configuration space. However, for most of 
these paths the exponential is tiny; the integral is dominated only by a small 
number of paths. In lattice simulations one tries to exploit this by sampling only 
a small number of gauge configurations that minimize the action using Monte 



3 LQCD is formulated in Euclidean space, which can be accomplished by a Wick rotation from 
Minkowski space t — > —He. We will use Euclidean space in this subsection only. 




Figure 2.1: Quark line diagram representing the pion correlation function. The hatched blobs 
represent sources for the pion at the indicated positions. The quarks are represented by solid 
lines. While only the quarks connected to the sources are shown, this diagram must be evalu- 
ated with all quark and gluon contributions. 

Carlo methods. Then one can approximate O(A) as the average over this finite 
ensemble of gauge configurations 



w))4e°^) (2 - 26) 



N 

c 

where N is the number of configurations in the ensemble. 

As an example, consider a pion that, being a pseudoscalar, can be represented by 
7v(x) = d(x)j 5 u(x) (see Fig. 12. ID . The correlation function for this pion can be written 
as 

(0|7f(x)vr(0)|0) = (0|n(x)7 5 d(x)d(0)7 5 u(0)|0) 

= — / VAVqVqu(x)~/ 5 d(x)d(0)~f 5 u(0) exp ( - / d 4 xC 

= ^JvAtr [G d (x, 0) 7 5 G U (0, ^) 7 5 ] exp (S G [A}) det[^+ m q ) 

\ L d\ i II u \ > ) i \ I sampled over gauge field configurations 

where Sq[A] is the pure Yang-Mills part of the gauge field action and G^ 1 = If)+ m q is 
the inverse propagator for a quark of flavor q [the discretized version of which appears 
in the last line of Eq. (I2.27M . In an LQCD simulation one approximates the expectation 
value by the average over the weighted samples. 

The computing power available today puts severe restrictions on what can be simu- 
lated: Typically, the size of the lattice, L, is limited to a few fermi (~2-4 fm); obviously 
it should be at least as big as the Compton wave length of the lightest particle one 
wants to simulate. Moreover, the lattice spacing, a, should be as small as possible so 
that discretization artifacts are kept to a minimum; typically a ~ L/5, so that a typical 
box size would be 5 x 5 x 5 x 5. 
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But even with these constraints it turns out that lattice simulations with realistic 
quark masses (m u ^ d ~ 5 MeV, m s ~ 100 MeV) are not feasible with the computational 
power that is available today. 

2.3 Quenching and Partial Quenching 

The fermion determinant in Eq. ( 12.271 1 is very expensive to compute since it typically 
scales ~ m~ 2 - 5 . In contrast, the propagators are much less costly to calculate as they 
scale ~ rriq 1 - 

The mass that appears in the fermion determinant is the mass for quarks that 
are generated in the gauge field background, i.e., it is only assigned to quarks that are 
generated dynamically from vacuum polarization in the gluonic background. These so- 
called "sea" quarks are not connected to the sources of the correlator. The quarks that 
are connected to the sources, and that have their mass appearing in the propagators, 
are called "valence" quarks. Since in a lattice simulation the calculation of the fermion 
determinant (that involves only sea quarks) is independent of the calculation of the 
propagators (involving solely valence quarks) one has the freedom to vary the masses 
of the sea and valence quarks independently. 

As an extreme way to save computing time one can omit calculating the fermion 
determinant completely. This is called the quenched approximation of QCD (QQCD). 
Effectively, this is a theory without sea quarks as they are treated as being infinitely 
heavy. Although simulating QQCD is much less costly than simulating full QCD (by 
a factor ~ 1000) it turns out that there exists, as will be explained shortly, no known 
connection between QQCD and QCD. Although there are hints that quenching might 
not make much difference for certain observables, it does introduces uncontrolled sys- 
tematic errors. 

A less severe approximation is partially quenched QCD (PQQCD). Unlike in QQCD, 
where the sea quark masses are set to infinity, they are kept finite in PQQCD. Sea 
quarks are thereby retained as dynamical degrees of freedom and the fermion deter- 
minant is no longer equal to one. However, by efficaciously giving the sea quarks 
larger masses, the fermion determinant becomes much less costly to calculate than in 
full QCD. The main advantage of PQQCD, compared to QQCD, is that there does exist 
a known analytic connection to QCD: By setting the sea quark masses equal to the 
valence quark masses one recovers QCD. 

How would one in practice do a perturbative QQCD or PQQCD calculation? The 
obvious approach is to write down all QCD Feynman diagrams that contribute to a cer- 
tain order in perturbation theory. Then, for QQCD, one simply disregards all diagrams 
that contain virtual quark loops since these consist of sea quarks. For PQQCD, one as- 
signs the sea quark mass to the quarks that appear in virtual loops. This method has 
been used in, for example, Refs. [11, 12]. Although dropping or modifying individual 
diagrams is very illustrative, this methods is somewhat artificial. 

A more systematic way, that does not require modification of individual diagrams, 
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is to include ghost quarks (that have bosonic statistics) in the quenched theory and 
to introduce ghost and sea quarks in the partially quenched theory. In the next two 
subsections we will introduce the field theoretical formulation of QQCD and PQQCD. 
We will also explain their effective low energy theories, QxPT and PQxPT, that are 
needed to properly extrapolate lattice data from the heavier light quark masses used 
on the lattice to realistic masses. 

Note that, although in general the number of valence and sea quark flavors need 
not be identical, we use the case of flavor 577(3) and work with three valence and 
three sea quark flavors throughout most of this thesis. The case of flavor SU(2), with 
two valence and two sea quark flavors, is very similar and will be explained when 
appropriate. 

2.3.1 QQCD and QxPT 

In QQCD the quark part of the Lagrangian is written as [13] 

a,b=u,d,s a,b=u,d,s j,k=u,d,s,u,d,s 

(2.28) 
Here, in addition to the fermionic light valence quarks u, d, and s their bosonic coun- 
terparts u, d, and S have been added. These six quarks are in the fundamental rep- 
resentation of the graded group 577(3|3) [14-16] and have been accommodated in the 
six-component vector 

Q = (u,d,s,u,d,s) (2.29) 

that obeys the graded equal-time commutation relation 

Qf (x)Qf (y) - (-l)^'Qf (y)Qf (x) = rt^ 3 (x - y), (2.30) 

where a and j3 are spin and i and j are flavor indices. The graded equal-time com- 
mutation relations for two Q's and two Q+'s can be written analogously. The grading 

factor 

J 1 for A: = 1,2,3 

% = t forfc = 4,5,6 (2 " 31) 

takes into account the different statistics for fermionic and bosonic quarks. The quark 
mass and charge matrices are given by 

m Q = diag(m u , m d , m s ,m u , m d , m s ) (2.32) 

and 

e = diag (l-i--i4 -!•-£)• <2 - 33) 

respectively, so that diagrams with closed ghost quark loops cancel those with valence 
quarks as illustrated in Fig. 12.21 
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Figure 2.2: Cancellation of valence and ghost quark loops in QQCD. Dashed lines represent 
ghost quarks. Since loops containing valence and ghost quarks of the same flavor have the 
opposite sign, the first two diagrams cancel completely, effectively removing any diagram with 
closed quark loops. 



For massless quarks, the Lagrangian in Eq. d2.28l > exhibits a graded symmetry 
SU(3\3)l <S> SU(3\3)r <S> U(1)v that is assumed to be spontaneously broken down to 
SU(3\3)v <8> U(l)v- The low-energy effective theory of QQCD that emerges by ex- 
panding about the physical vacuum state is QxPT. The dynamics of the emerging 
36 pseudo-Goldstone mesons can be described at lowest order in the chiral expansion 
by the 0{E 2 ) Lagrangian 4 [17-21] 

C = ystr (Ve^e) + A str (m Q E + m^J + ad»$ d^ - figsg (2.34) 



where £ is defined in Eq. ( 12.151 1 and 



<I> 



71" X 
X 7T 



(2.35) 



Here the ir, tt, and x are 3x3 matrices of pseudo Goldstone bosons with quantum 
numbers of qq pairs, pseudo Goldstone bosons with quantum numbers of qq pairs, and 
pseudo Goldstone fermions with quantum numbers of qq pairs, respectively: 



7] u 7T + K + 
71 = | 7T _ Tfd K° 

K~ K° t] s 



TT 



flu 7T+ K + 

R- K° fj s 



Xvu Xtt+ Xk+ 
and x= ! X-k- X Vd Xro 

XK- Xk° Xr] s 

(2.36) 
The pion decay constant is / = 132 MeV, and we have defined the gauge-covariant 
derivative Z? M S = <9 M £ + ieA^lQ, E]. The str() denotes a supertrace over flavor indices 
defined as 



str(X) = £(-l)wx«. 



(2.37) 



i=l 



*Here, E ~ p, m n where p is an external momentum. 
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Upon expanding the Lagrangian in H2.34H one finds that to lowest order the mesons 
with quark content QQ' are canonically normalized when their masses are given by 

m QQ' = "72 ( m Q + m Q')- ( 2 - 38 ) 

One also finds that the propagator for off-diagonal (flavored) Goldstone mesons com- 
posed of (ghost-) quarks Q and Q' is given by 

GQQ'(P) = 2-—-- ( 2 - 39 ) 

P ~ m QQ> + le 

The flavor-singlet field <J>o is defined as 

<S> = -L s tr(<S>) = -L( v '-fi'). (2.40) 

$o is invariant under SU(3\3)l <8> SU(3\3)r (g) U(l)v and thus arbitrary functions of it 
can be included in the Lagrangian. To lowest order in the chiral expansion only the two 
operators included in Eq. d2.34D with parameters a and fio remain and are understood 
to be inserted perturbativly [19]. Notice that this singlet field <J>o is not heavy as in 
%PT and therefore cannot be integrated out. It introduces a new vertex, the so-called 
hairpin with the propagator 

G VaVa = 2 \ . + - ^~ ap2 \ 2 (2.41) 

nana p 2_ m 2 a+K ^2_ m 2^ +ie y 

that exhibits a double pole which causes quenching artifacts and is ultimately respon- 
sible for the sick behavior of the quenched theory. 

2.3.2 PQQCD and PQxPT 

The physics for flavor off-diagonal mesons in PQQCD is very similar to the QQCD 
case. The quark part of the Lagrangian is extended once again by including three light 
fermionic sea quarks j, I, and r and can be written as [22-29] 

£ = ^2 Qa(Hp-m q ) ab q b + ^ 9a( i] P- m q)ab% + ^2 ?Bea,o(*#- ^sea)abQsea,6 

a,b=u,d,s a,b=u,d,s a,b=j,l,r 

J2 Qj(ip-m Q ) jk Q k . (2.42) 

j,k=u,d,s,u,d,s,j,l,r 

These nine quarks are in the fundamental representation of the graded group 
SU (6\3) [14-16] and have been accommodated in the nine-component vector 

Q = (u,d,s,j,l,r,u,d,s) (2.43) 
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Figure 2.3: Cancellation between valence and ghost quark loops in PQQCD. As in Fig. 12.21 
dashed lines represent ghost quarks whereas fat solid lines represent the (heavier) sea quarks. 
Like in QQCD, valence quark loops are canceled by their ghostly counterparts. The inclusion of 
the (fermionic) sea quarks effectively replaces the valence quark masses in loops by sea quark 
masses. Obviously, in the QCD limit, where valence and sea quarks have equal mass, one 
recovers QCD. 



that obeys the graded equal-time commutation relation in Eq. A2.30D . Now, however, 
the grading factor is 

1 for k = 1,2,3,4,5,6 

for k = 7, 8, 9 

The quark mass matrix is given by 



'Ik 



(2.44) 



m Q = diag(m u , m d ,m 



s , mj,mi,m r ,m u , m d , m s ) 



(2.45) 



so that, in a perturbative expansion, diagrams with closed ghost quark loops cancel 
those with valence quarks just like in QQCD. Effects of virtual quark loops are, how- 
ever, present due to the contribution of the finite-mass sea quarks (see Fig. 12.31 1. 

It has been recently realized [30] that the light quark electric charge matrix Q is 
not uniquely defined in PQQCD. The only constraint one imposes is for the charge 
matrix Q to have vanishing supertrace. Thus, as in QCD, no new operators involving 
the singlet component are subsequently introduced. Following [31] we use 



Q = diag 



>Qj,QhQr,Qj,Ql,q r 



m u , mi -* m d . 



and 



111; 



111 



(2.46) 



independently 



so that QCD is recovered in the limit nij 
of the q's. 

For massless quarks, the Lagrangian in Eq. ( 12.421 1 exhibits a graded symmetry 
SU(6\3)l <S> SU(6\3)r <8> U(l)y that is assumed to be spontaneously broken down to 
SU(6\3)y <S> U(l)y. The low-energy effective theory of PQQCD that emerges by ex- 
panding about the physical vacuum state is PQxPT. The dynamics of the emerging 
80 pseudo-Goldstone mesons can be described at lowest order in the chiral expansion 
by the Lagrangian given in Eq. H2.34H with £ as denned in Eq. ( 12.151 ) but <I> now being 
extended to include mesons that contain sea quarks 



$ 



M X t 
X M 



(2.47) 
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The M, M, and x are matrices of pseudo-Goldstone bosons with quantum numbers 
of qq pairs, pseudo-Goldstone bosons with quantum numbers of qq pairs, and pseudo- 
Goldstone fermions with quantum numbers of qq pairs, respectively. Explicitly they 
are given by 



/ 



M 



>1u 



IT 



K+ J° L+ R+ \ 



and 



7T 


>ld 


K° 


j- 


L° 


R° 






K~ 


J+ 


Vs 

J? 


Vj 


Y + 


R° s 

Y + 


' 


M -- 


L~ 


L u 


L» 


Y ii 


Vi 


Y° 






R- 


R° 


R° s 


Y~ 


1 lr 


Vr ) 














f X V 


u A.7T + 


Xk+ 








X = 


1 Xn 

\ Xk 


Xr ld 

Xro 


Xk» 

Xris 




(2.48) 



(2.49) 



Meson masses and non-singlet propagators are similar to the quenched case as given 
in Eqs. d2~38t and J2~39l >. 

The flavor singlet field given by $o = str($)/v / 6 is, in contrast to the QxPT case, 
rendered heavy by the U(1)a anomaly and can therefore be integrated out in xPT. 
Analogously, its mass fio can be taken to be on the order of the chiral symmetry break- 
ing scale, no — ► A x . In this limit the flavor singlet propagator becomes independent of 
the coupling a and deviates from a simple pole form [22, 23]: 



G 



iS 1 



ab 



(q 2 - mf,)(9 2 - m 2 rr ) 



3d> 



>h i}b 



q 2 - m\ a +ie 3 (q 2 - m 2 a + ie) (q 2 



m 



m, 



+ ie) (q 2 — m 2 x + ie) 



(2.50) 



This can be more compactly written in a form that only contains single poles: 
G r . 



'VaVb 



6 a bp{m } l K-QKr-Q p(m , 



+ - 



I (mjj - m 2 b )(m 2 r 



m. 



Vb ' 



3 (m 2 a - m 2 b ){m 2 b - m 2 x 



p ( m %)-- 



1 (m x - m 2 d ){m 2 x - m 2 r ) 



3 (m 2 x - m 2 a )(m 2 x - m 2 b ) 



P{m x ) 



where 



and mx is given by m 2 x 



P(m) = - 
(rrijj + 2m 2 r )/3. 



m l + ie 



(2.51) 



(2.52) 
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2.3.3 Inclusion of the Baryon Octet and Decuplet in Q^PT and PQxPT 

Just as there are mesons in QQCD [PQQCD] 5 with quark content QiQj that contain 
valence (v) and ghost (g) [v, g, and sea(s)] quarks, there are baryons with quark com- 
positions QiQjQk that contain these two [three] types of quarks. Restrictions on the 
baryon fields Bij k come from the fact that these fields must reproduce the familiar octet 
and decuplet baryons when i, j, k = 1-3 [11,31,32]. To this end, one decomposes the 
irreducible representations of SU(3\3)y [SU(6\3)y] into irreducible representations of 
SU(3) V ® S77(3) g <g> U{\) [S77(3) v ® SU(3) S ® SU(3) g ® 17(1)]. 

Baryon Octet 

The method to construct the octet baryons is to use the interpolating field 

fl&* ~ (QTQfQY - QTQTQf] e abc (c 75 u, (2.53) 



ijk y*i **j 

which when restricted to i, j, k = 1-3 has non-zero overlap with the octet baryons. 



Under SU(3\3)y [SU(6\3)y], where Qi — ► UijQj and Qi — > QjUL, Bijk transforms as 



B ijk — {-) Vm ^ + ^ )+{Vm+ ^ ){m+m) U im U jn U k iB mnl . (2.54) 

Using the commutation relations in Eq. 12.30B one sees that B^k satisfies the symme- 
tries 

&ijk I, ) ^ikji 

= B ijk + (-) 1+ ™iB jik + (-) 1+ *w +r >^ +r >MB kji . (2.55) 

The spin- 1/2 baryon octet B^ k = Bij k , where the indices i, j, and k are restricted to 
1-3, is contained as an (8, 1) [(8, 1, 1)] of SU(3) Y ® SU(3) g [SU(3) V ® SU(3) S SU(3) g ] 
in the 70 [240] representation. The octet baryons, written in the familiar two-index 
notation 



iA+ i£° E+ 



P 



B=\ E" ^A-^E° n I (2.56) 



20 -* A 



are embedded in B^k as [11] 



Bijk = —7= (eijiBki + eikiBji) . (2.57) 

As explained in Ref. [11], it is convenient to switch to the three-index "quark flow" 
notation B^k as opposed to the familiar two-index notation of the octet baryons. The 

5 Here, we explain the inclusion of baryons for the quenched case; the partially quenched case is very 
similar and included in square brackets. 
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Table 2.1: Embedding of the baryon octet and decuplet into SU(3\3) V for QQCD. 





Octet 




Decuplet 






SU(3) y <g> SU(3) g 


dim 


SU(3\®SU(3) g 


dim 


qqq 


(8,1) 


8 


(10,1) 


10 


qqq 


(6,3)0(3,3) 


27 


(6,3) 


18 


qqq 


(3,6)8(3,3) 


27 


(3,3) 


9 


qqq 


(1,8) 


8 


(1,1) 


1 




70 


38 



Table 2.2: Embedding of the baryon octet and decuplet for SU(6\3) V for PQQCD. 





Octet 




Decuplet 






SU(3) v ®SU(3) B ®SU(3) g 


dim 


SJ7(3) V <g> SU(3) B <8> 


SU(3) g 


dim 


qqq 


(8,1,1) 


8 


(10,1,1) 




10 


qqq s 


(6,3,1)0(3,3,1) 


27 


(6,3,1) 




18 


qq s q s 


(3,6,1)0(3,3,1) 


27 


(3,6,1) 




18 


q s q s qs 


(1,8,1) 


8 


(1,10,1) 




10 


qqq 


(6,1,3)0(3,1,3) 


27 


(6,1,3) 




18 


qq s q 


(3,3,3)0(3,3,3) 


54 


(3,3,3) 




27 


q s q s q 


(1,3,6)0(1,3,3) 


27 


(1,6,3) 




18 


qqq 


(3,1,6)0(3,1,3) 


27 


(3,1,3) 




9 


q s qq 


(1,3,6)0(1,3,3) 


27 


(1,3,3) 




9 


qqq 


(1,1,8) 


8 


(1,1,1) 




1 




240 


138 



reason that the two-index notation is possible at all is due to the fact that a 3 x 3 matrix 
contains 8 elements plus an overall constant. 

Besides the conventional octet baryons that contain valence quarks, qqq, there are 
also baryon fields with other types of quarks contained in the 70 (240). Since we 
are only interested in calculating one-loop diagrams that have octet baryons in the 
external states, we will need only the Bijk that contain at least two valence quarks. 
We use the explicit construction in [13,31]. For example, baryons that consist of two 
valence and one ghost quark are denoted by the tensors ahc and ahc that transform as 
a 27 = (6, 3) + (3, 3) of SU(3) V ® SU(3) g [27 = (6, 1, 3) + (3, 1, 3) of SU{3) y ® SU(3) S ® 
SU(3) g ]. For completeness, we list the transformations for octet baryons containing 
any combination of quarks in Table O for QQCD and in Table O for PQQCD. In 
Appendix|A|we list the transformations of the doublet and quartet baryons for the two 
flavor case SU{2) (Tables lAH and IX2b. 
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Baryon Decuplet 

Similarly, the familiar spin-3/2 decuplet baryons are embedded in the 38 [138]. Here, 
one uses the interpolating field 

T$ ~ (QTQfQY + Q? ' b Q]' c Q a k ' a + QT c QT a Qk b ) *** (CV% 7 (2.58) 

that describes the 38 [138] dimensional representation of S77(3|3)y [S77(6|3)y] and 
has non-zero overlap with the decuplet baryons when the indices are restricted to i, j, 
k = 1-3. Due to the commutation relations in Eq. H2.30D . %j k satisfies the symmetries 

%jk = {~) l+ri ^T 3ik = (-Y^m %kr (2.59) 

The decuplet baryons are then readily embedded in T by construction: T^ k = %j k , 
where the indices i, j, k are restricted to 1-3. They transform as a (10, 1) [(10, 1, 1)] un- 
der SU (3) v <g> SU (3) g [S77(3) v ® SU{3) S ® SU (3) g l Because of Eq. <2~59T ). T ijk is a totally 
symmetric tensor. Our normalization convention is such that Tin = A ++ . For the spin- 
3/2 baryons that contain two valence quarks — the only ones relevant for our purpose — 
we use the states constructed in [13,31]. For example, spin-3/2 baryons consisting of 
two valence and one ghost quark transform as (6, 3) [(6, 1, 3)] under S77(3) v <g> S77(3) g 
[S77(3) v (g> S77(3) s <3 S77(3) g ]. For completeness, we list the transformations for the 
remaining decuplet baryons in Table O(QQCD) and Table O(PQQCD); the trans- 
formations for the two flavor case are given in Appendix lAl 

Free Lagrangian for Baryons 

At leading order in the heavy baryon expansion, the free 577(3) Lagrangian for the B^ k 
and T ijk is given by [11] 

C = i(Bv- VB) + 2a M (BBM+) + 2f3 M (BM+B) + 2a M (BB) str (M+) 

-i (T% • PT M ) + A (T%) + 2 lM {T^M+T^) - 2a M {T%) str (M+) , 

(2.60) 

where M+ = \ (tfmcjtf + £toq£) with £ 2 = S. The covariant derivatives of Bij k and 
%j k both have the form 

(V»B) ijk = d»B ijk + (V»)*B ljk + n Vi{Vj+Vm) (V% m B imk + (-)(vi+n i )(v k + rin )(y*) knBi . nm 

(2.61) 
The brackets in 1 12. 60H are shorthands for field bilinear invariants originally em- 
ployed in [11] 

(BTB) = B kji T^B ijk ^, (2.62) 

(BTYB) = B a kji TiY u B ljk ^ (BTBY) = (-)^ + ^^ + ^B a kji T^Y kn B ijn ^, (2.63) 

{T^TT^) =T il k ^Ti% ]k ^ (2.64) 

(T^TYT,) =Tt%TiY tl T ljk ^ and (Brr%) = B^r£ (Y») u T ljk ^ (2.65) 
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which ensure that the contraction of flavor indices maintains proper transformations 
under chiral rotations. To lowest order in the chiral expansion, Eq. ( 12.601 1 gives the 
propagators 

v ■ k v ■ k — A 

for the spin- 1/2 and spin-3/2 baryons, respectively. Here, v is the velocity and k the 
residual momentum of the heavy baryon which are related to the momentum p by 
p = Mbv + k. Mb denotes the (degenerate) mass of the octet baryons and A the 
decuplet-baryon mass splitting. The polarization tensor 

p^ v = ( v ^v u - g^) - -S»S U (2.67) 

reflects the fact that the Rarita-Sch winger field T? k contains both spin- 1/2 and spin-3/2 
pieces; only the latter remain as propagating degrees of freedom (see [33], for example). 

2.4 Extrapolation of Lattice Data 

If unquenched lattice simulations with light enough quarks were possible today then 
one could simply use x?T to extrapolate to the physical quark masses. Unfortunately, 
now and in the foreseeable future this is not the case and one is bound to simulate 
using the quenched or partially quenched approximations and to extrapolate to the 
physical quark masses using the appropriate low-energy effective theories, QxPT and 
PQxPT. The next question then is: What statements about QCD can be made from 
extrapolated QQCD or PQQCD lattice data? 

Since PQQCD retains a U(1)a anomaly, the equivalent to the singlet field in QCD 
is heavy (on the order of the chiral symmetry breaking scale A x ) and can be inte- 
grated out [22, 23] — -just like in QCD. Therefore, the low-energy constants appearing 
in PQxPT are the same as those appearing in xPT. By fitting PQxPT to partially 
quenched lattice data one can determine these constants and make physical predic- 
tions for QCD. The advantage of PQQCD is that, since one can vary the sea quark 
masses independently from the valence quark masses, one has an enlarged parameter 
space (more adjustable "knobs") and can hope to determine the low-energy constants 
with greater accuracy by fitting to a larger number of partially quenched lattice results 
(see Fig, mi). 

For example, since the valence and ghost quarks have equal masses, the contribu- 
tion of valence quarks in disconnected quark loop diagrams is eliminated by the ghost 
quarks. The effects of disconnected loop diagrams are solely due to sea quarks and the 
physics of the sea sector can be explored by varying the sea quark masses. Further- 
more, in the limit where the masses of the sea quarks become equal to those of the 
valence and ghost quarks, one recovers QCD. 

In processes that involve electroweak gauge fields, the "theory space" of PQQCD is 
enlarged once more since one can chose arbitrary values for the charges of the ghost 
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Figure 2.4: Domains for QCD and PQQCD in quark mass space. Full QCD is represented 
by the point on the diagonal. xPT "lives" on the diagonal. Partial quenching opens up the 
parameter space to the 2-dimensional gray region that has a large overlap with the region 
where PQxPT is applicable. The low-energy parameters for PQxPT (which are the same as 
for %PT) can be determined from this overlap with improved precision. Figure from Ref [34] 
courtesy Ruth Van de Water. 



and sea quarks, qj, qi, and q r in Eq. H2.46I ). For example, if one choses qj = qi = q r = 
then photons can only couple to valence quarks. In the case qj = 2/3, q\ = —1/3, and 
q r = —1/3 contributions of the valence and ghost sectors cancel and photons can only 
couple to sea quarks [32]. 

In QQCD the answer to the question raised above is different. The problem with the 
quenched approximation is that the Goldstone boson singlet is no longer affected by the 
U(1)a anomaly as in QCD. In other words, the QQCD equivalent of the rf that is heavy 
in QCD remains light and must be included in the QxPT Lagrangian. This requires the 
addition of new operators and hence new low-energy constants. In general, the low- 
energy constants appearing in the QxPT Lagrangian are unrelated to those in x?T 
and extrapolated quenched lattice data is unrelated to QCD. Although there is some 
empirical evidence that the difference between QQCD and QCD for some observables 
is small at large quark masses both theories deviate considerable in the low quark 
mass region. In fact, several examples show that the behavior of meson loops near 
the chiral limit is frequently misrepresented in QxPT [2, 5, 13, 35]. In the following 
chapters, we find this is additionally true for a number of meson and baryon hadronic 
properties. 
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Chapter 3 

Chiral 1/M 2 corrections to 

#(*) _^ £)(*) at Zero Recoil in QxPT 



In this chapter, we study the semileptonic B^ — ► Z)M decays in the limit that the 
heavy quark masses are infinite. We calculate the lowest order chiral corrections, 
which are of 0(1/M 2 ), from the breaking of heavy quark symmetry at the zero recoil 
point in QxPT. These results will aid in the extrapolation of quenched lattice calcula- 
tions from the light quark masses used on the lattice down to the physical ones. 

3.1 Introduction 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix describes the flavor mixing among 
the quarks; its elements are fundamental input parameters for the standard model. 
Their precise knowledge is not only crucial to determine the standard model but also 
to shed light on the origin of CP violation. The matrix element that parametrizes the 
amount of mixing between the b and c quarks, V c b, can be extracted from the exclusive 
semileptonic B meson decays B — ► Dlv and B — ► D*lu, where I = e,[i. Heavy quark 
effective theory (HQET) (for a recent review, see [36]), which is exact in the limit of 
infinite masses M for the heavy quarks, predicts the width of the process B — ► D*lu as 

dT_ {B ^ Dl = $0_ k{u):Fb ^ (w)2| (g.D 

where u = v-v' is the scalar product of the 4-velocities v and v' of the B and D* mesons, 
respectively. K.(u) is a known kinematical factor and F(u) is a form factor whose value 
at the kinematical point w = 1 is JF(1) = 1 in the M — ► od limit. There are, however, 
perturbative and nonperturbative corrections to .F(l), 

Tb^d* (1) = f]A + ^i/m 2 + • ■ ■ , (3.2) 

where the parameter t]a ~ 0.96 is a QCD radiative correction known to two-loop or- 
der [37] and d 1 / M 2 are non-perturbative corrections of 0(1/M 2 ) to the infinite mass 
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limit of HQET. Note that, according to Luke's theorem [38], there are no 0(1/M) cor- 
rections at zero-recoil. One chooses the zero-recoil point because, for uj = 1, Tb^,d* can 
be expressed in terms of a single form factor Ka\ given by 

(D'Mjcr^Bw = _ aA (iy ., (33) 

^m B m D * 

This is in contrast with the general case uj > 1 for which Fb^d*{u) is a linear com- 
bination of several different form factors of B — ► D*lv mediated by vector and axial 
vector currents. 

Several experiments, most recently by CLEO [39], have determined the product 
[Fb->d* (l)|V^t,|] 2 by measuring dT b^d* / 'duj and extrapolating it to the zero-recoil point. 
The mixing parameter |V^,| can then be extracted once the value Tb^d* (1) that en- 
codes the strong interaction physics has been evaluated. The uncertainty in |V^| is 
therefore determined by the experimental errors and by theoretical uncertainties in 
the determination of Tb^d* (!)• Presently, the theoretical uncertainties dominate. 1 

A model-independent way of calculating ;F(1) is provided by numerical lattice QCD 
simulations. Recently, such calculations have been performed [40-43] for the decays 
B — ► D^*'lv using QQCD. Several systematic uncertainties, such as from statistics and 
lattice space dependence, contribute to the error of these calculations. Another con- 
tribution to the uncertainties comes from the chiral extrapolation of the light quark 
mass. This extrapolation can be done by matching QQCD to QxPT and calculating 
the non-analytic corrections 5\/m^ m Eq. 03.21 ) in QxPT. The formally dominant con- 
tributions to these corrections come from the hyperfine mass splitting between the 
heavy pseudoscalar and vector mesons that stems from the inclusion of heavy quark 
symmetry breaking operators of 0(1/M) in the Lagrangian. 

In QCD, the corrections due to D meson hyperfine splitting have been calculated in 
xPT by Randall and Wise [44]. A more complete treatment, involving additional cor- 
rections due to B meson hyperfine splitting, 0(1/M) axial vector coupling corrections, 
and 0(1/M) corrections to the current, has been given in [45]. Recently, the D meson 
hyperfine splitting corrections have also been determined in PQxPT [46] for PQQCD. 

In this chapter, we calculate the 0(1/M 2 ) corrections in Eq. iS.'Zi due to D and 
B meson hyperfine splitting in QxPT. These corrections are — upon expanding in pow- 
ers of the hyperfine splitting A — of order Aqcd 3 ™/ 2 / (M n mq ) for n > 2 and for- 
mally larger than those coming from the inclusion of 0(1/M) heavy quark symmetry 
breaking operators in the Lagrangian and current which are suppressed by powers of 
Aq CD /A/. This argument is similar to the one that applies to x?T [36]. Our QxPT 
calculation can be used to extrapolate lattice results [42] that use the quenched ap- 
proximation down to the physical light quark masses. So far, this extrapolation has 

Similarly, one can use the decay B — > Dlu to extract [JPb^d(1)|V C 6|] 2 from the measured dTB^D/duj. 
However, dYs—D/dco is more heavily suppressed by phase space near oj = 1 than dTs^D* /du. In ad- 
dition, the B — > D channel is experimentally more challenging. Thus the extraction of | V c b \ from this 
channel is less precise but serves as a consistency check. 
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been based upon the ^PT calculation [44]. Using Q^PT should therefore give a better 
estimate of the uncertainties related to the chiral extrapolation. 

A central role in the lattice calculation of B — ► D* [42,43] is played by the double 
ratios of matrix elements 

_ {D\c 7 °b\B)(B\h°c\D) 
+ (D\afc\D)(B\b7%\BY 

_ (D*\c 7 b\B*)(B*\h c\D*) 
1 (D*\c~i°c\D*)(B*\b~i%\B*)' 

and 

_ (D*\c^ l5 b\B)(B*\b^ l5 c\D) 

Al (D*\cj^ 5 c\D){B*\b^-f 5 b\B)' 

In these ratios, statistical fluctuations are highly correlated and cancel to a large de- 
gree. The 0(1/M 2 ) correction to the double ratios can therefore be calculated fairly 
accurately and used to derive the C(l/M 2 ) correction to the matrix elements them- 
selves. For this reason, we also calculate C(l/M 2 ) corrections to the decay B* — ► D* 
in addition to the experimentally accessible decays B — ► D and B — ► D*, and thus the 
corrections to 1Z + , TZ±, and IZax- 

3.2 Quenched Heavy Meson Chiral Perturbation Theory 

The D mesons with quantum numbers of cQ can be written as a six-component vector 

D = {D u , D d , D s , Du, Dg, D- s ). (3.7) 

Heavy quark symmetry is provided by combining creation and annihilation operators 
for the pseudoscalar and vector mesons, D and D*, respectively, together into the field 
H D : 

H D = ^(W+ il5 D), (3.8) 

S D = l o H m 7 o = (r t +i75l) t ) i+i ) (3.9) 

where v denotes the velocity of a heavy meson. In HQET the momentum of a heavy 
quark is only changed by a small residual momentum of O(Aqcd)- Hence, v is not 
changed and H is usually denoted by an index v which we have dropped here to un- 
clutter the formalism. In the heavy quark limit, the dynamics of the heavy mesons are 
described by the Lagrangian [47, 48] 

C D = -itr[H°v^8 ab + iV£)H?] + gtr(H?H° 7ul5 A» ba ) + 7 tr(H?H? 7/i75 ) strA» 

(3.10) 
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where the traces tr() are over Dirac indices and supertraces str() over the flavor indices 
are implicit. The additional coupling term involving $o ~ strA^ is a feature of QxPT 
and not present in ^PT. The light-meson fields are 

\ = \(£}d^ - &$) = -l d ^ + o{&) (3 n) 

and 

V, = ^dpt + td^) = ^,d^] + 0(^). (3.12) 

Expanding the Lagrangian Co to lowest order in the meson fields leads to the (deriva- 
tive) couplings DD*d(p and D*D*d(f> whose coupling constants are equal as a conse- 
quence of heavy quark spin symmetry. At leading order in the 1/M expansion, the 
DDdcft coupling vanishes by parity. 

An analogous formalism applies to the fields B and B* which are combined into 
H B . Note that the axial coupling g is the same for H D and H B mesons at this order in 
the 1/M expansion as dictated by heavy quark flavor symmetry. 

We do not include terms of order m q ~ ^m n in the Lagrangian as explicit chiral 
symmetry breaking effects are suppressed compared to the leading corrections. The 
presence of these terms is implied by the nonzero meson masses m qq . 



3.3 Matrix Elements of B& _> £>(*)/ 



v 



The non-zero hadronic matrix elements for B^ — ► Z)W can be defined in terms of the 
16 independent form factors h±, hy, hA 123 , and /11...10 as [36,49] 

(D(v')\c^b\B(v)) = Mw)(w + v , r + Mw)(u _ ^ (3 13) 



(D*(v',e')\c^b\B(v)) _ ^vap,* , , 



(3.14) 



<£>V, OI^T^M = _^ i(w)(w + 1)e „ + ^ >)( , . e> „ + ih {u){v . £> ^ 



^m B m D * 

(D*(v',e')\c^b\B*(v,e)) 
^m B *m D * 



(3.15) 

-(e'* • e)[h!(u)(v + v'T + h 2 (u)(v - v'f] + h 3 (u)(e'* ■ v)e" 

+h A {io){e ■ v')e'*» - (e • v')(e'* ■ v)[h 5 (u)v» + h 6 (u)v >fX ], 

(3.16) 
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and 

(gVif p^g M l = i£ ^ {e :e 5 [h 7 (u)(v + vr + hs(u)(v-vr] 
^rriB*mD* 

+ v' a v p [h 9 (co)(e'* ■ v)e' M + h 1Q (u)(e ■ v')e'*»] } . 

(3.17) 

Here, u = v ■ v' and v (e) and v' (e') are the velocities (polarization vectors) of the 
initial state B^ meson and final state D^ meson, respectively. Note that we will 
not explicitly calculate matrix elements of B* — > D as these can be easily related to 
the B — » D* calculation by a Hermitian conjugation of the matrix elements and an 
interchange of the c and b quarks, i.e., fiW <-► £)(*). 

In the heavy quark limit the matrix elements in Eqs. d3.13tt - d3.17tt are reproduced 
by the operator 

erf {I - 75 )k - -£(u;)tr[ZT„V(l - 7s)^]. (3.18) 

Here, £(w) is the universal Isgur-Wise function [50,51] with the normalization £(1) = 1. 
To lowest order in the heavy quark expansion one finds 

h + (uj) = hy(uj) = iiA^u) = Ha z {oj) = h\(uj) = h%(uj) = h±(u) = hj(u) = £(w) (3.19) 

and the remaining 8 form factors vanish. 

The discussion of the B^ — > D^*>W matrix elements is similar for different flavors 
of the light quark q content of the i?M and JDM mesons; it applies equally to q = u, d, 
or s as the theory splits into three similar copies of a one-flavor theory. In the limit of 
light quark SU(3)y flavor symmetry the matrix elements (and in particular the Isgur- 
Wise function) are therefore independent of the light quark flavor. However, in nature 
the masses of the u, d, and s quarks are different and SU(3)v is not an exact symmetry. 
Therefore our results will include terms that depend upon m q via the meson masses 



/?? 



qq 



defined in Eq. ( EOSl 



3.4 1/M 2 Corrections 

The lowest order heavy quark symmetry violating operator that can be in- 
cluded in the Lagrangian Co in Eq. A3. 101 ) is the dimension-three operator 
(Ac 2 /Afo)tr ( H a a iW H.^a ilv J . It violates heavy-quark spin and flavor symmetries and 

comes from the QCD magnetic moment operator ca tlu G A u T A c, where G A is the gluon 
field strength tensor and T A with A = 1, ... ,8 are the eight color ST/ (3) generators. 
This operator gives rise to a mass difference between the D and D* mesons of 

A D = m D * -m D = -8^. (3.20) 

Md 
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Figure 3.1: Graphs contributing to wavefunction renormalization for heavy (a) pseudoscalar 
and (b) vector mesons. A thin (thick) line denotes a heavy pseudoscalar (vector) meson, a 
dashed line denotes the $o, while a dashed-crossed line denotes the insertion of a hairpin. A 
full (empty) vertex denotes a g (7) coupling. 

This effect can be taken into account by modifying the D and D* propagators which 
become 

and -^^-v.vu) 



2(v ■ k + 3A D /4 + it) 2(v ■ k - A D /4 + ie) 

respectively, so that in the rest frame, where v = (1, 0, 0, 0), an on-shell D has residual 
energy of — 3A^/4 and an on-shell D* has residual energy of A_d/4. A similar effect 
due to the inclusion of a QCD magnetic moment operator for the b quark applies to the 
B(*> mesons. 

There are no corrections to the matrix elements for the semileptonic decays B^*> — > 
Z)(*W of <D(1/M) at zero-recoil according to Luke's theorem [38]. The leading correc- 
tions enter at C(l/M 2 ). In addition to tree-level contributions from the insertion of 
0(1/M 2 ) suppressed operators into the heavy quark Lagrangian or the current there 
are one-loop contributions from wave function renormalization and vertex correction. 
These one-loop diagrams have a non-analytic dependence on the meson mass m qq and 
depend on the subtraction point /i. This dependence on /i is canceled by the tree-level 
contribution of the C(l/M 2 ) operators. 

Because of the absence of disconnected quark loops in QQCD, which manifests it- 
self as a cancellation between intermediate pseudo Goldstone bosons and pseudo Gold- 
stone fermions in loops in QxPT, the only loop diagrams that survive are those that 
contain a hairpin interaction or a 7 coupling [see Eq. d3. 10D 1. 

The wave function renormalization contributions for the pseudoscalar and vec- 
tor meson, Z D / B and Z* D , B , respectively, come from the one-loop diagrams shown in 
Fig. 13. II and have been calculated in [35,47,48]. Including the a coupling we find, for 
these diagrams, 



P 



1 + ^tfi(A) - ^H 2 (A) + ^F X (A) (3.22) 
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Figure 3.2: QxPT graphs which contribute to the vertex correction of the form factors (a) h+(l), 
(b) h,A 1 {l), and (c) hi(l). A full (empty) square denotes the insertion of the operator C7 M 7 5 & 

(cy&). 



and 



l + ^I ffl (-A)-i^ ff2 (-A) + ^F l( -A) 



+ -^H l (0)-^H 2 (0) + ^?F l (0). (3.23) 



The functions Hi, Hi, and Fi come from loop integrals and are given in Appendix iBl 
Note that in the heavy quark limit where A = one recovers Z = Z* , as required by 
heavy quark symmetry. 

The vertex corrections come from one-loop diagrams. The non-vanishing contribu- 
tions are shown in Fig. 13.21 Combining the wave function renormalization and vertex 
corrections and including a local counterterm to cancel the dependence on the renor- 
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malization scale fx, we find the following corrections for the form factors: 



X+W + n + S 



Qig-y 



l j^[^ H 5 (A B ,A D ) - aH 8 (A B ,A D )} - ^G 5 (A B ,A D ) 



x+(v) + 



(IA 



(4vr/) 2 \3m 2 
+0({A B ,A D } 3 ), 



%-*» 



m 



F 



log ^- ] (A B - A D ) 2 



(3.24) 



Sh Al (l) 



and 

Stn(l) 



XaM + 



Z b -1 



+ 



J D 



1 



3/ 2 

3/2 



2 2 

2 2io^ 

[/ig F 5 (A B , -A D ) - a H 8 (A B , -A D )} - -^Lg 5 (A b , -A d ) 



P 



4ig~/ 



[t4 H 5 (A B ,0) - a H 8 (A B , 0)] - -^G 5 (A B , 0) 



X Al (») + 



(i& 



(4tt/) 2 \9m 2 
+0({A B ,A D } 3 ), 



g 2 a 



2gj 

3 



log ^j] (3A| + A 2 D + 2A B A D ) 



(3.25) 



Xi(n) + ^- r 



1 + ^JL 



3/2 

_2iy2 

S/ 2 



[nl H 5 (-A B , -A D ) - a H 8 (-A B , -A D )] 



^G 5 (-A B ,-A D ) 



[n 2 H 5 (0,0)-aH 8 (0,0)) 



1 (<?& 
(4tt/) 2 V9m 2 

+0({A B ,A D } 3 ), 



4ig7 
/ 2 

ff 2 « _ 2jn 
9 3 



G 5 (0,0) 



log^-lCAB-A,,) 2 



(3.26) 



which are defined by h + (l) = 1 + 5h + (l) and analog expressions for Sh Al (1) and 5hi(l). 
The functions i/5, -ffs> and G5 come from loop-integrals that are listed in Appendix iBl 
and we have defined m = m qq . The insertions of tree-level C(l/M 2 ) operators are 
represented by the functions X+(/x), X Al (//,), and Xl(//) which are independent of m 
and exactly cancel the \x dependence of the logarithm. These functions can be extracted 
from lattice simulations by measuring the zero-recoil form factors for a varying mass 
of the light quark. 

Experimentally, Ad ~ 142 MeV and A B ss 46MeV so that the ratios A^/ra and 
A B /m, which enter the form factor corrections through the function R(A/m) (defined 
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in Appendix [Bll, are 0(1). On the lattice, however, one can vary all quark masses. Ex- 
panding first in powers of A and then taking the chiral limit m — » one finds the formal 
limits given in Eqs. H3.24t - d3.26t where we have only kept the pieces non-analytic in to. 
This demonstrates that the terms linear in A^ and As, although present in wave func- 
tion renormalization and vertex corrections, cancel as required by Luke's theorem [38]. 
The leading order corrections are 0({A#, A^} 2 ). 

As a consistency check one can restore heavy quark flavor symmetry by taking 
As = Ad- Since the 0(1/M 2 ) corrections to h+(l) and hi(l) are proportional to (A# — 
Ad) 2 they disappear as they should since the charge associated with the operators cj^c 
and 67^6 is conserved. This argument does not apply for the B —> D* transition matrix 
element in the limit A# = Ad since there is no conserved axial charge associated with 
the operators 07^750 and 67^756. 

In the chiral limit, the term proportional to /i 2 , has a 1/to 2 singularity and domi- 
nates over the terms proportional to a and 7 that are only logarithmically divergent. 
This is analogous to a term of the form (m 2 q — m 2 j) / 'm qq found by Savage [46] for 
PQxPT (here, m qq and rrijj are valence and sea quark masses, respectively). In the 
limit rrijj — ► m qq this term, however, vanishes as PQxPT goes to xPT where the dom- 
inant term is logm qq . In QxPT, on the other hand, the l/m 2 q pole persists, revealing 
the sickness of QQCD where the hairpin interactions give a completely different chiral 
behavior than in QCD. 

The size of /xo can be estimated from the 77-7/ mass splitting [18], large Nc ar- 
guments [52, 53] (Nc being the number of colors), or lattice calculations. These es- 
timates imply /x « 500 — 900 MeV; for the purpose of dimensional analysis we use 
/io ~ 0(Aqcd)- Taking g ~ 0(1) we therefore find that Sh+, 5Ha 1 , and 8h\ are of order 
A n /m n ~ Aqcd 3 ™^ 2 / (M n m q ) for n > 2 and thus larger than tree-level heavy quark 
symmetry breaking operators that are suppressed by Aq CD /M. 

To show the dependence of the zero-recoil form factors on the mass of the light 
spectator quark it is necessary to know the numerical values of the parameters /xo, 9, 
a, and 7. In determining reasonable values for these couplings we follow the discussion 
by Sharpe and Zhang [48]. Assuming that g is similar to the xPT value we use g 2 = 0.4. 
The hairpin coupling a is proportional to 1/Nc, and thus assumed to be small; we use 
two values, a = and a = 0.7. The coupling 7 is known to be suppressed by 1/Nc 
compared to g, the sign is undetermined. We take —g<j<g (see [48] and references 
therein). 

With these parameters, the dependence of h + (l) and hA 1 (l) on the mass of the 
light spectator quark m q is show in Figs. 13.31 and 13.41 respectively. The graphs are 
plotted against m q in units of the strange quark mass m s with m q /m s = m 2 /m 2 where 
m 2 s = Im 2 ^ . The behavior of h + (l) in QxPT is dominated at small m by the 1/m 2 pole 
that is non-existent in xPT. Lattice calculations of h + (l) [41] show a small downward 
trend for decreasing m q down to the chiral limit that is similar to the downward trend 
seen from the xPT calculation (dashed line). The same behavior (down to m q ki 0.1m s ) 
can also be seen for QxPT for a certain choice of parameters (e.g., 7 positive). The 
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Figure 3.3: Dependence of h + (l) on the mass m q of the light spectator quark in QxPT. For com- 
parison, the xPT result from [44] is also shown (dashed line). The result has been normalized 
to unity for m q = m s . We have chosen ^o = 700 MeV and g 2 = 0.4. 



case of hA 1 (l) is different as there is a pole at m = Ap which is close to the physical 
pion mass. Here, both D* and it can be on-shell and the decay B —> D*tt becomes 
kinematically allowed. Lattice calculations of Ua x (1) [42] for m q = (0.6 . . . l)m s show a 
small downward trend for decreasing m q similar to the downward trend seen from the 



xPT calculation (dashed line in Fig. 13.41 ). A similar trend down to m q « 0.2m s can also 
be seen in the QxPT calculation for a relatively large positive value of 7. 

Although the downward trend in the lattice data for the two cases seems significant 
as the statistical errors are highly correlated, the uncertainty is still relatively high 
(typically ±0.01) and the existing lattice data can be accommodated by a wide range of 
values for the parameters in the QxPT Lagrangian. 

As can be seen in the figures, the variation of the quenched result is primarily due 
to the parameter 7 as the sensitivity of the result to the value of a is very small. We 
have also checked how the result depends on the parameter g in the reasonable range 
0.1 < g 2 < 0.5 and found that the change from the value g 2 = 0.4 is at most 25% for 
g 2 = 0.5, still well within the statistical errors of the lattice data. 

Finally, we calculate the double ratios defined in Eqs. ( I3.4l >-( l3~6b using the results 
in Eqs. d3^4t -( l3~26t . We find 



K+ = l + 26h+(l), 



(3.27) 
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Figure 3.4: Dependence of h Al (1) on the mass of the light spectator quark in QxPT. The dashed 
line denotes the %PT result [44]. The numerical values for the parameters are those used in 
Fig. 1331 



and 



Ki = l + 2Shi(l) 



(3.28) 



l + X Al (ji) 



' (l {n\ [H B (A B , -A d ) + H 5 (A D , -A B ) - H 5 (A D , -A D ) - H 5 (A B , -A B )] 
- a [H 8 (A B , -A D ) + H 8 (A D , -A B ) - H S (A D , -A D ) - H 8 (A B , -A B )]} 



"3/2 



2ig-f 



P 
l+^i(M) 



[G 5 (A B , -A D ) + G 5 (A D , -A B ) - G 5 (A D , -A D ) - G 5 (A B , -A B )} 



1 



2, ,2 



Wj4 

9m 2 



2g 2 a 4g7 
9 3 



(4vr/)2 
+0({A B ,A D } 3 ), 

where X Al (//) is the counterterm associated with K Al 



log^^A^A,,! 2 



(3.29) 



3.5 Conclusions 

Knowledge of the B(*> — > D^> form factors at the zero-recoil point is crucial to ex- 
tract the value of V& from experiment. In the limit that the heavy quarks are 
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infinitely heavy, HQET predicts that the form factors h+, h^, and h\ are equal, 
/i + (l) = /iAi(l) = ^i(l) = £(!)• The formally dominant correction due to breaking of 
heavy quark symmetry comes from the inclusion of a 0(1/M) dimension-three opera- 
tor in the Lagrangian that leads to hyperfme-splitting between the heavy pseudoscalar 
and vector mesons. These leading order corrections are 0({A#, A^} 2 as required by 
Luke's theorem. 

Recent lattice simulations using the quenched approximation of QCD have made 
a big step forward in determining these zero-recoil form factors. Presently, however, 
the simulations use light quark masses that are much heavier than the physical ones 
and therefore rely on a chiral extrapolation down to the physical quark masses. In this 
chapter we have calculated the dominant corrections to the form factors h+ , Jia 1 , and 
h\ in QxPT and determined the non-analytic dependence on the light quark masses 
via the light meson masses m qq . Using these results, instead of the %PT calculation, to 
extrapolate the QQCD lattice measurements of these form factors down to the physical 
pion mass should give a more reliable estimate of the errors associated with the chiral 
extrapolation. We have also calculated the corrections to certain double ratios that are 
used in lattice QCD calculations of the decay B — ► D*. 
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Chapter 4 

HMxPT in a Finite Volume 



In this chapter, we study finite volume effects in heavy quark systems in the frame- 
work of heavy meson chiral perturbation theory (HMxPT) for QCD, QQCD, and 
PQQCD. A novel feature of this investigation is the role played by the scales A and 
S s , where A is the mass difference between the heavy-light vector and pseudoscalar 
mesons of the same quark content, and S s is the difference of the masses of the u and 
d, and the mass of the s quark that is due to light flavour 51/(3) breaking. The primary 
conclusion of this chapter is that finite volume effects arising from the propagation of 
Goldstone mesons in the effective theory can be altered by the presence of these scales. 
Since A varies significantly with the heavy quark mass, these volume effects can be 
amplified in both heavy and light quark mass extrapolations. 

As an explicit example, we present results for B parameters of neutral B meson 
mixing matrix elements and heavy-light decay constants to one-loop order in finite 
volume HM^PT for full, quenched, and Nf = 2 + 1 partially quenched QCD. Our cal- 
culation shows that for high-precision determinations of the phenomenologically in- 
teresting 51/(3) breaking ratios, finite volume effects are significant in QQCD and not 
negligible in PQQCD, although they are generally small in QCD. 

4.1 Introduction 

Numerical calculations of hadronic properties using LQCD have provided significant 
inputs to particle physics phenomenology. In particular, the joint effort between ex- 
periment and theory to investigate the unitarity triangle in the CKM matrix from 
B meson decays and mixing has made impressive progress [54], in which LQCD has 
played an important role. Nevertheless, current lattice calculations are still subject to 
various systematic errors. In this chapter, we address finite volume effects which arise 
in lattice calculations for heavy-light meson systems from the light degrees of freedom. 
Our framework is HMxPT with first order 1/M and chiral corrections. We assume the 
mass hierarchy 

m QQ/ < A x < M, (4.1) 

34 



where rngg/ is the mass of any Goldstone meson given in Eqs. ( 12.201 ) and ( I2.38l > and M 

is the mass of the heavy-light meson. Under this assumption, we discard corrections 

of the size 

rnqy (42) 

M 

Concerning the finite volume, we work with the condition that 

m Q Q,L > 1, (4.3) 

where L is the spatial extent of the cubic box. Therefore, given that f n L/\/2 (f n r> 
132 MeV) will be close to one in lattice simulations in the near future, one can still 
neglect the chiral symmetry restoration effects resulting from the Goldstone zero mo- 
mentum modes [55,56] when Eq. ( 14.31 1 is satisfied. 

The main task of this work is to study the volume effects due to the presence of the 
scales 

A = M P * - M P , (4.4) 

and 

5 S = M Ps - M P , (4.5) 

where P* and P are the heavy-light vector and pseudoscalar mesons containing a u or 
d anti-quark, 1 and P s is the heavy-light pseudoscalar meson with an s anti-quark. The 
scale A appears due to the breaking of heavy quark spin symmetry that is of 0(1/M) 
and S s comes from light flavour SU(3) breaking in the heavy-light meson masses. Un- 
der the assumption of Eq. d4.ll >. A is independent of the light quark mass, and 5 S does 
not contain any 1/M corrections, at the order we are working. 

In the real world, both A and 5 S are not very different from the pion mass. In 
fact [1], M Bs - M B = 91 MeV, M Ds - M D = 104 MeV, M B * - M B = 46 MeV, and 
Mp>* — Mp, = 142 MeV. In current lattice simulations, these mass splittings vary be- 
tween and ~ 150 MeV. Therefore it is important to include them in the investigation 
of finite volume effects. Equation J4.3t > implies that the Compton wavelength of the 
Goldstone meson is small compared to the size of the box. Therefore finite volume 
effects mainly result from the propagation of the Goldstone mesons to the boundary. 
However, as shown in Section 14731 A and 5 S can, in a non-trivial way, alter these effects. 
In particular, since A varies with the heavy quark mass, finite volume effects can be 
significantly amplified in heavy quark mass extrapolations. 

This chapter is organized as follows. In Section l4~2l we summarize the ingredients 
of HM^PT relevant to this work. We mainly expand the treatment in Section 13.21 to 
the partially quenched case. Section l4~3l is devoted to the discussion of HMxPT in fi- 
nite volume, emphasising the role played by S s and A. We then present an explicit 
calculation of neutral B meson mixing and heavy-light decay constants in Section 14.41 
and discuss the phenomenological impact that finite volume effects can have. We con- 
clude in Section 14.51 Some mathematical formulae and results are summarized in 
Appendix O 

J We work in the isospin limit. 
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4.2 Partially Quenched Heavy Meson Chiral Perturbation 
Theory 

The chiral Lagrangian for the Goldstone mesons in the three theories QCD, QQCD, 
and PQQCD has been discussed in detail in Chapter |2] and we will not repeat it here. 
The inclusion of the heavy-light mesons into the quenched version of HMxPT has been 
discussed in Section 13.21 Here we will expand this discussion and also include the 
cases of QCD and PQQCD. 

HMxPT was first proposed in Refs. [57-59], with the generalization to quenched 
and partially quenched theories given in Refs. [47,48]. The 1/M and chiral corrections 
were studied by Boyd and Grinstein [60] in QCD and by Booth [35] in QQCD. The 
spinor field appearing in this effective theory is 

HP = ±±* te«V - PPhs) , (4.6) 



a, (j, 

where Pa and Pa,]?' annihilate pseudoscalar and vector mesons containing a heavy 
quark Q and a light anti-quark of flavour a. 

The HM^PT Lagrangian, to lowest order in the chiral and 1/M expansion, for D 
mesons containing a heavy quark c and a light anti-quark of flavour a is then given 
in Eq. ( I3.10D . The Lagrangian for a heavy B meson containing a b quark is analogous. 
In QCD, the term proportional to 7 is non-existent as the rf is integrated out. We do 
not formally distinguish the coupling g in the three theories with the understanding 
that the numerical values are different. The HMxPT Lagrangian for mesons contain- 
ing a heavy anti-quark Q and a light quark of flavour a is obtained by applying the 
charge conjugation operation to the Lagrangian in Eq. A3. 10l > [61]. At this order, the 
propagators for P a and Pa mesons are 

and -'(ayVfr), (4.7) 



2(vk + ie) 2(v-k + ie) 

respectively. 

The effects of chiral and heavy quark symmetry breaking have been systemati- 
cally studied in full [60] and quenched [35] HMxPT Amongst them, the only relevant 
feature necessary for the purpose of this work, i.e., the investigation of finite volume 
effects, are the shifts to the masses of the heavy-light mesons. These shifts are from 
the heavy quark spin breaking term 

and the chiral symmetry breaking terms 



at) 



) 



bb 

(4.9) 
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We choose to work with the effective theory in which the heavy-light pseudoscalar 
mesons that contain a heavy quark and auord valence anti-quark are massless. 
Notice that the term proportional to X[ in Eq. d4.9l > causes a universal shift to all the 
heavy-light meson masses. This means that the masses appearing in the propagators 
of heavy vector mesons and any meson containing an s anti-quark (valence or ghost) 
are shifted as follows: 

-i(g„u - Vf,v u ) i and -i( gfiu - yy) (4lQ) 



2{v ■ k - /S. + ie) 1 2(vk-6 s + ie)' " 2(v ■ k - A - S s +ie)' 

for P* , P s , and P*, respectively. The mass shifts can be written in terms of the cou- 
plings in Eqs. <4~8t and <4~9l l. A = -8A 2 /M , and 

5 S = 2X 1 B (m s - m u ). (4.11) 

In PQQCD, there are two additional mass shifts because the sea quarks have dif- 
ferent masses from the valence and ghost quarks: 

S s = M Ps -M p = 2X 1 B (m r -m j ), (4.12) 

and 

<5 sea = Mp-M P = 2X 1 B (m j - m u ). (4.13) 

where P (P s ) is the heavy-light pseudoscalar meson with a d (s) sea anti-quark. The 
propagators of the heavy mesons containing sea anti-quarks are: 



-i(j9fw ~ Vf.Vy) 



2{vk- <5 S ea + ie) 2(v ■ k - A- 5 sea + ie) 



(4.14) 



1 —„ , and -^ ~ V ^\ (4.15) 

2{v ■ k - (5 sca - 5 s +ie) 2(v ■ k - A - 5 sca - 5 S + ie) 

for P, P* (vector meson with a d sea anti-quark), P s , and P* (vector meson with an s 
sea anti-quark), respectively 

4.3 Finite Volume Effects 

In this section, we discuss generic features of finite volume effects in HMxPT. For 
clarity, we use the symbol A for one of (A, 5 S , S s , 5 sca ) or any sum amongst them. 

In the limit where the heavy quark mass goes to infinity and the light quark masses 
are equal, all the heavy mesons in HM^PT become on-shell static sources, and there is 
a velocity superselection rule when the momentum transfer involved in the scattering 
of the heavy meson system is fixed [62]. For illustration, consider the vertex with cou- 
pling g in the Lagrangian in Eq. ( I3.10D . The heavy-light meson P can scatter into PT-. 
by emitting a Goldstone meson with mass itiqq/ through this vertex. The momenta of 
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the mesons P and P/*% are Mpv^, and Mp* v ll + k^ = Mp-u^ + feu, where the velocity 
Vfj, = (1,0,0,0) in the rest frame of the heavy mesons, and k^ is the soft momentum 
carried by the Goldstone meson. The infinitely heavy P and P/*% mesons do not prop- 
agate in space. Therefore, when such a system is in a cubic spatial box, finite volume 
effects result entirely from the propagation of the Goldstone meson to the boundary 
with momentum k ~ mQQ>- In this case the volume effects behave like exp (—rnqqiL) 
multiplied by a polynomial in 1/L. 

The breaking of heavy quark spin and ST/ (3) light flavour symmetries in HM^PT 
can induce a mass difference Mp* = Mp + A, which complicates the above picture. 
In this scenario, the P/\ is still regarded as a static source, but it is off-shell with 
the virtuality A. The period during which the Goldstone meson can propagate to the 
boundary is limited by the time uncertainty conjugate to this virtuality, i.e., 

St ~ -J-. (4.16) 

This means that finite volume effects, which arise from the propagation of the Gold- 
stone mesons in such a system, will decrease as A increases. Eq. ( I4.16l > also indicates 
that the suppression of the volume effects by a non-zero A is controlled by the param- 
eter 

IW (4 . 17) 

To see explicitly how this phenomenon appears in a calculation, we consider a typ- 
ical sum in one-loop HMxPT, with a Goldstone propagator and a heavy-light vector 
meson propagator in the loop, in a cubic box with periodic boundary conditions: 

1 ^ f dk 1 



J(m, A) = -i T3 l^J -^ ik 2_ m 2 + ie){v . k _ A + ie y W.18) 

k 

where the spatial momentum k is quantized in finite volume as k = (2ir/L)i, with i 
being a three dimensional integer vector. Using the Poisson summation formula, it is 
straightforward to show that J(m, A) = J(m, A) + J^vim, A), where 

/d 4 k 1 
7 7T77^ n w i 2 T> (4.19) 

(2vr) 4 (k 2 - m 2 + ie){v ■ k - A + ie) ' 

is the infinite volume limit of J(m, A), and 

T / ^ 1 V^ t°° n^ $ sm(\k\\n\L) , Ann . 

J FV (m,A) = - ni J2 d\k\ ' J— " -J \ (4.20) 

47r z z — iJ ujyjj + A) \n\L 



with uj = y | k\ 2 + m 2 , is the finite volume correction to J(m, A). In the asymptotic limit 
where mL » 1 it can be shown that (with n = \n\) 

J F y(m, A) = Y^ £ f exp(-nmL) A, (4.21) 
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where 

,_ ., „ 1 



.A = exp(z 2 )[l-erf(z)] + 



nmL 
1 



' ' /? '""' N +^-exp(z 2 )[l-erf(z)] 



vr 14 2 / 2 



1 /9z 5z 3 7z 5 z 7 \ fz & z 8 . , r „ . 



(nmL)'- 

+° (t^t^ ) ■ (4 - 22) 



with 



A I nmL 
m 



(4.23) 



The quantity A is the alteration of finite volume effects due to the presence of a non- 
zero A. It multiplies the factor exp(— nmL), which results from the propagation of the 
Goldstone meson to the boundary. It is possible to analytically compute the higher 
order corrections of A in powers of l/(nmL). This way, one can achieve any desired 
numerical precision. Here it is clear that this alteration of volume effects is controlled 
by the ratio in Eq. d4.17l >. 

Next, we consider different limits of A at fixed m and L. When A = 0, clearly A = 1. 
If A is very small compared to m, such that z <C 1, A is dominated by the [l/(mL)]° 
term, i.e., A ~ exp(z 2 )[l — erf(z)]. Since erf(z) grows much faster than exp(z 2 ) in this 
regime, A will decrease as A increases. When A is of 0(m) or larger, z » 1, and 
one can perform an asymptotic expansion of the error function. It can be shown that 
in this situation, A ~ 1/z. That is, A also decreases as A increases. We have also 
numerically checked that this is true when z « 1. This means that the asymptotic 
formula in Eq. ( I4.21l > reproduces the physical picture outlined in the beginning of this 
section for any A. To demonstrate how fast the asymptotic form in Eq. 04.221 1 converges 
to Eq. H4.20I ). we define 

JT r ax JFV nUm KA)- JFv as y m P(m,A) . Aa .. 

djFv(m < A) = W) ' (4 - 24) 

where Jpv num is the function Jpy evaluated numerically [Eq. ( 14.201 11. and JFv asymp is 
the asymptotic form in Eq. A4.22I ). In Fig. 14.11 we plot dJpy as a function of m with 
three choices of A. It is clear from this plot that Jpy is approximated well (to < 3%) 
by the asymptotic form when mL > 2.5. We use the asymptotic forms for integrals 
of this type throughout this work. Also, in this paper we only include the terms with 
\n\ = l,\/2, \/3, \/4, and \/5 in the Poisson summation formula. We have confirmed that 
truncating the sum at \n\ = y/5 is a very good approximation (to ~ 3%) when mL > 2.5. 
The function J-py(m,A) is plotted against m in Fig. 14.21 with L = 2.5 fm and three 
choices of A. It is clear from this plot that A can significantly alter the finite volume 
effects in J(m, A). 
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Figure 4.1: dJyv(m,A) plotted against m for three different A. This function indicates the 
deviation (in percent) of the asymptotic form of Jfv from the definition in Eq. <4.20l l. In this 
plot, L =2.5 fm and m = 0.197 GeV corresponds to ml = 2.5 whereas m = 0.32 GeV corresponds 
to mL = 4. The curve for A = 0.1 GeV is hidden by that for A = 0.2 GeV. 



Another typical sum that appears in one-loop HM^PT in finite volume is 

i x-^ f dk 1 



/C(m,A) 



L 3 



E 



2tt (k 2 - m 2 + ie){v ■ k - A + ie) 2 ' 



It is straightforward to show that )C(m, A) = K(m, A) + K^tm, A), where 

d 

K(m,A) = —J(m,A) 



is the infinite volume limit of K,(m, A) and 







KYv(m,A) = — JFv(m,A) 



(4.25) 



(4.26) 



(4.27) 



is the finite volume correction to K(m, A). The function K^yim, A) is plotted against 
m in Fig. 14.31 with L = 2.5 fm and three choices of A. As expected, \KYvim, A)\ also 
decreases when A increases at fixed m and L. 

4.4 Neutral B Mixing and Heavy-Light Decay Constants 

The study of neutral B meson mixing allows the extraction of the magnitude of the 
CKM matrix element Vtd, and hence the determination of one of the sides of the uni- 
tarity triangle. The frequency of the B^—B^ oscillations, which is given by the mass 
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Figure 4.2: J$y(m, A) plotted as a function of m, for three different A. Again, L =2.5 fm. 

difference, Arrid, in this mixing system has been well measured by various experimen- 
tal collaborations [54]. It is also calculable in the standard model via an operator prod- 
uct expansion in which the top quark and W boson are integrated out. Resumming the 
next-to-leading order (NLO) short-distance QCD effects, one obtains 



Am d = ^M^\V td V* b \ 2 VB S {xt)C B {^ 



8vr 2 



\(B d \Q* B = 2 (n)\B d )\ 

2Mb 



(4.28) 



where /x is the renormalisation scale, x t = mf/M^, and So(x t ) ~ 0.784x° 76 (to better 
than 1%) is the relevant Inami-Lim function [63]. The coefficients tjb = 0.55 and Cb{^) 
are from short-distance QCD effects [64, 65]. The matrix element of the four-quark 
operator 

0% B = 2 = [67^(1 - 75MMI " 75)4 (4.29) 

between B d and B d states contains all the long-distance QCD effects in Eq. H4.28H . and 
has to be calculated non-perturbatively Since \V t b\ = 1 to good accuracy and Am d has 
been well measured, a high-precision calculation of (B d \0^ B=2 (n)\B d ) enables a clean 
determination of |T4d|- 

The frequency of the rapid B s —B s oscillations can be precisely measured at the 
Tevatron and LHC [54]. Therefore an alternative approach is to consider the ratio 

(4.30) 



Am, 



Am d 



Vt 



ts 



V 



td 



M Bd 



M Ba 



(B s \0 



(B d \Of B = 2 \B d ) 



in which many theoretical uncertainties cancel. Here Am s is the mass difference in 
the B s -B s system and Of B=2 = [67^(1 - 75 )s] [67^(1 - 75 )s]. The unitarity of the CKM 
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Figure 4.3: Kw(m, A) plotted as a function of to for three different A. Again, L = 2.5 fm. 

matrix implies \V ts \ ~ \V c b\ to a few percent, and \V c b\ can be precisely extracted by 
analysing semileptonic B decays [54]. Therefore a clean measurement of Am s /Amd 
will yield an accurate determination of \Vtd\- 

The matrix elements in Eq. H4.30I I are conventionally parameterized as 



(B q \0^\B q ) = | 



Mlf 2 B B B U), 



(4.31) 



where the parameter Bs q measures the deviation from the vacuum-saturation approx- 
imation of the matrix element, and q = d or s. The decay constant fB q is defined by 

(0\h^ 5 q\B q (p)) = ip„f Bq - (4.32) 

LQCD provides a reliable tool for calculating these non-perturbative QCD quanti- 
ties from first principles. 2 Since Am s /Amd will be measured to very good accuracy, 
it is important to have clean theoretical calculations for [the SU(3) breaking ratios 
of] the matrix elements, decay constants and B parameters involved. Current lattice 
calculations have to be combined with effective theories in order to obtain these ma- 
trix elements at the physical quark masses. This procedure can introduce significant 
systematic errors and dominate the uncertainties in the 51/(3) breaking ratio [75, 76] 



£ 



Ib 3 \/Bb s 

fB\[B~B 



(4.33) 



2 Some selected reviews in the long history of lattice calculations for the B mixing system can be found 
in Refs. [66-74]. 
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Figure 4.4: Wave-function renormalization diagrams for heavy-light mesons 

which is the key theoretical input for future high-precision determination of \Vu\ via 
the study of neutral B mixing. 3 However, the use of effective theory also offers a 
framework for studying finite volume effects in lattice calculations [18,26,77-89]. We 
will demonstrate in this section that finite volume effects might turn out to exceed the 
current quoted systematic errors for quantities such as £. 

4.4.1 One-Loop Calculation in a Finite Volume 

Here, we discuss one-loop calculations for the B parameters and heavy-light decay con- 
stants mentioned above in finite volume HM^PT including the appropriate mass shifts 
to the first non-trivial order of the chiral and 1/M expansion. The inclusion of other 
first-order corrections in these quantities is straightforward. It simply introduces ad- 
ditional low-energy constants (LECs) which account for short-distance physics and do 
not give rise to finite volume effects at this order, so we will not discuss this issue here. 
We have performed the calculation for QCD, QQCD and PQQCD with the mass shifts 
given in Eqs. d4 . 1 0I >— ( f4~T5b . 

For the purpose of this work, the axial current 67^75^ is 



Aft = |tr D 



wX Q) zl b 



(4.34) 
and the four-quark operator o KPa=2 (when P a = B dtS , Q APa=2 becomes Off =2 ) is 



Qaa = A p 



£P*w)tV Up*(Q)nY + UpmY UpiQ) 



(4.35) 



in HMxPT [61], where n and f3 are the low-energy constants which have to be deter- 
mined from experiments or lattice calculations. Notice that the index a in Eq. H4.35K 
is not summed over. Again, the inclusion of the chiral and 1/Mp corrections in these 
operators simply introduces additional LECs and we do not investigate this aspect 
here. We assume that n and f3 are the same in QCD, QQCD, and PQQCD. Also, A^ and 
Qaa can C0U pi e t the rf in QQCD, but the couplings are 1/N C suppressed [48], and we 
neglect them. 

The diagrams contributing to fp, s) and Bp {s) are presented in Figs. 14 .41 - 14 .61 Note 
that only diagrams with an intermediate heavy meson depend on the heavy meson 
mass shifts. 



3 Notice that the symbol £ as defined in Eq. I4.33J is in the traditional notation in B physics, and has 
nothing to do with the Goldstone field £. 
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Figure 4.5: Diagrams contributing to the one-loop calculation of decay constants. The thin 
(thick) solid lines are the heavy-light pseudoscalar (vector) mesons. The dashed lines are 
Goldstone mesons, and the crosses are the "double poles" which appear in (P)QxPT. The open 
squares are the operators defined in Eq. J4.34t and the dots are vertices from the HM^PT 
Lagrangian. 
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Figure 4.6: Diagrams contributing to the one-loop calculation of the B parameters. The open 
circles are the operators defined in Eq. <4.35t . 



Although this is the first one-loop calculation for these decay constants and B pa- 
rameters in finite volume, some results in the infinite volume limit already exist in the 
literature: fp (s) have been calculated at the lowest order in QCD [61], QQCD [47,48], 
and PQQCD [48], and up to first-order corrections in the chiral and 1/M expansions 
in QCD [60] and QQCD [35]. The B parameters have been calculated only at lowest 
order [48,61]. Our results, as presented in Appendix lC.21 agree with all these previous 
calculations in the appropriate limits. 



4.4.2 Phenomenological Impact 

We have used the one-loop results in Appendix IC.2I to investigate the impact of finite 
volume effects on £. In this work, we only intend to estimate the possible size of errors 
in this quantity, and will leave the actual comparison with lattice data to a future 
publication. Following the usual procedure in lattice calculations for £, we study two 
SU(3) breaking ratios 

£/ = ^ (4.36) 



and 



Cb 



f'B 

Bb 



(4.37) 



in terms of which, £ = £/v£b- Furthermore, we define (£/)fv an d (£b)fv to be the 
contributions from finite volume effects, i.e., those from the volume-dependent part 
in the one-loop results presented in Subsection 14.4. II To be more precise, we use the 
formulae collected in Appendix IC.2I to calculate the volume corrections with respect 
to the lowest-order values of fs s (Bb 3 ) and fs (Bb), then take the difference between 
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the results as an estimate of (£/)fv [(£b)fvL Since these SU(3) ratios are not very 
different from unity (at most ~ 20%), this is a reasonable estimate of these effects. 

Traditionally, many quenched lattice simulations of Bb {s) and fs (s) were performed 
using L ~ 1.6 fm. Therefore we present our estimate for finite volume effects in QQCD 
with this box size. For comparison, we adopt the same volume for QCD. As for PQQCD, 
we work with L = 2.5 fm where most current high-precision simulations are carried 
out [90]. Throughout this subsection, we ensure that the condition M n L > 2.5 holds in 
all the plots presented here. 

We first discuss the procedure in QCD and QQCD. When studying the light quark 
mass dependence of (£/)fv and (£_b)fv> we follow a strategy similar to that in Ref. [75]. 
That is, we use the Gell-Mann-Okubo formulae to express Mk and M v in terms of M„ 
and m,,: 



? +Ml 2 2m 2 x + M 2 

>k g and M * :] 



Mi = ^^- ?- and Mi = -^-! *-. (4.38) 



We investigate the situation where a lattice calculation is performed at the physical 
strange quark mass (m s ) phys , but the up and down quark mass m u = m d is varied. 
By using (M K ) phys = 0.498 GeV and (M 7r ) phys = 0.135 GeV [1], we fix (m ss ) phys = 
2B (m s ) piiys = 0.691 GeV as an input parameter in our analysis. Notice that (m ss ) p h ys is 
not the mass of a "physical" meson, and the subscript just means this mass is estimated 
by using physical kaon and pion masses. To the same order, we can adopt Eq. ( 14. lit 
to write S s = Ai (m 2 ss - M 2 ), and use {m ss ) phys , (M^)^ and physical M Bs - M B = 
0.091 GeV [1] to determine 

Ai = 0.1982 GeV" 1 . (4.39) 

This determines how 5 S varies with M n . We have also tried to use vanishing pion mass 
and Mo s — Mo = 0.104 GeV [1] to fix (m ss ) p h ys and Ai, and the results presented in this 
subsection are not sensitive to this variation from the values quoted above. 

The results for (£/)fv and (£b)fv for QCD and QQCD from this analysis are pre- 
sented in Figs. l4.7H4.10l with two different values for the coupling g (and also 7 in 
QQCD). Here we stress again that the influence on finite volume effects from the pres- 
ence of A and 5 S depends on the size of these couplings, which are not well determined. 
Inspired by the recent CLEO measurement of g in the charm system [91,92], and a 
recent lattice calculation [93], we vary g 2 between 0.3 and 0.5. As for the coupling 
7, which is a quenching artifact and has never been determined, we vary its value 
between g and —g. It is clear from these plots that the finite volume effects are gen- 
erally small in QCD (< 2%), but can be significant in QQCD (~ 3% to ~ 7% for £ B ) in 
the range of M^L where lattice simulations are normally performed. This is clearly 
due to the enhanced long-distance effect arising from the "double pole" structure in 
(P)QQCD, as first pointed out in Ref. [77], and manifests itself in various places, e.g., 
nucleon-nucleon potentials [94] and it— it scattering [77,81,82,95]. 

Although it has been well established that infinite volume chiral corrections are 
smaller in the B parameters than in the decay constants due to the coefficient in front 
of g 2 in the one-loop results, it is clear from these plots that finite volume effects are 
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Figure 4.7: (£/)fv in full QCD plotted against M w with L = 1.6 fm. The pion mass M T 
0.35 GeV corresponds to M^L = 2.8, and M w = 0.5 GeV corresponds to M n L = 4. 



full QCD (£b)fv 
0.0175 



full QCD (f B ) F v 




0.0025 



g =0.3 
L=1.6fm 

A, = 0.00 GeV 
A, =0.10 GeV 
A. = 0.20 GeV 



M„ (GeV) 




g =0.5 

L=1.6fm 
A, = 0.00 GeV 
A, =0.10 GeV 
A, = 0.20 GeV 



M„ (GeV) 



Figure 4.8: (6?)fv in ful1 Q CD Pitted against M n with L = 1.6 fm. 

more salient in £# than in £f. All the quenched lattice calculations for £# have so far 
concluded that this quantity is consistent with unity with typically 3% error. However, 
we find that the volume effects are already at the level of 3— 4% when M n = 0.45 GeV in 
a 1.6 fm box where many quenched simulations were carried out. This error depends 
on both light and heavy quark masses in the simulation, hence is amplified after ex- 
trapolating the result to the physical quark masses. Also, the fact that volume effects 
tend towards different directions in QCD and QQCD when M w becomes smaller indi- 
cates that quenching errors in these quantities can be larger than those estimated in 
Ref. [48]. Since finite volume effects have not been included in the analysis of lattice 
calculations of £b hitherto, one should be cautious when using the existing quenched 
results for this quantity in any phenomenological work. 

For the analysis in PQQCD, we assume that both the valence and sea strange quark 
masses are fixed at that of the physical strange quark. However, we vary the light sea 
quark mass rrij = mi. For this purpose, we define rrijj to be the mass of the meson 
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Figure 4.9: (£/)fv in QQCD plotted against M„, with L = 1.6 fm and different couplings g and 

7. We set a = and M = 700 MeV. 



composed of two light sea quarks. Therefore, 



rn 



03 



(m 2 ss ) phys \m 



m i 



(4.40) 



m r = physical m s 



Also, we can express the mass shifts 8 S and 5 sca in terms of meson masses, 5 S = 
Ai (m 2 ss — m 2 .^ and 5 sca = Ai (m 2 r — M^) by using Eqs. H4.12H and ( 14.13ft with the same 
value for Ai used in Eq. H4.39H . 

The results for the PQQCD analysis are presented in Figs. 14.111 and 14.121 The 
double pole insertions also appear in PQQCD and it is clear from these plots that fi- 
nite volume effects cannot be neglected if one hopes to determine £ to the level of a few 
percent. Especially, in the range of M n and L where current and future lattice simula- 
tions are performed [90], they can already be at about 4%, and the dependence on the 
heavy meson mass is quite strong. Therefore they can become comparable to the error 
presented in the latest review [74], £ = 1.23 ± 0.10 after quark mass extrapolations. 4 

4 Finite volume effects presented in this work are, however, correlated with the errors arising from 
chiral extrapolations. 
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Figure 4.10: (6j)fv in QQCD plotted against M v , with L = 1.6 fm, a = and M = 700 MeV, 
and different couplings. 

4.5 Conclusions 

We have investigated finite volume effects in heavy quark systems in the framework 
of HM^PT. The primary conclusion is that the scales A and 5 S , which are heavy-light 
meson mass splittings arising from the breaking of heavy quark spin and light flavour 
SU(3) symmetries, can significantly reduce the volume effects in diagrams involving 
heavy meson propagators in the loop. The physical picture of this phenomenon is that 
some heavy-light mesons are off-shell in the effective theory, as a consequence of the 
velocity superselection rule, with the virtuality characterized by the mass splittings. 
The time uncertainty conjugate to this virtuality limits the period during which the 
Goldstone mesons can propagate to the boundary. Finite volume effects caused by the 
propagation of the Goldstone mesons naturally affect the light quark mass extrap- 
olation/interpolation in a lattice calculation. On top of this, our work implies that 
they also influence the heavy quark mass extrapolation/interpolation, since the scale 
A varies significantly with the heavy meson mass. The strength of this influence is 
process-dependent, determined also by the relative weight between diagrams with and 
without heavy meson propagators in the loop. The volume effects can be amplified by 
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Figure 4.11: (£/)fv in PQQCD plotted against m rr , with L = 2.5 fm and two different values 



for M„. The pion mass M% = m 2 ss /A corresponds to M^L = 4.4 and M 2 = 



corresponds 



to M^L = 3.1. The mass m rr = 0.197 GeV corresponds to m rr L = 2.5, and m rr = 0.32 GeV 
corresponds to m rr L = 4. 



both heavy and light quark mass extrapolations. Therefore it is important to perform 
calculations to identify these effects in phenomenologically interesting quantities. 

We have presented an explicit calculation in finite volume HM^PT for the B pa- 
rameters in neutral B meson mixing and heavy-light decay constants, in QCD, QQCD, 
and PQQCD. We have used these results to estimate the impact of finite volume effects 
in the SU(3) ratio £, which is an important input in determining the magnitude of the 
CKM matrix element Vu- Within the parameter space where most quenched lattice 
calculations have been performed, we find that, although this impact is quite small 
(<~ 2%) in QCD, it can be significant in QQCD. This is due to the enhanced long- 
distance effects arising from the double pole structure. This error will be amplified by 
the quark mass extrapolations and hence can exceed the currently quoted systematic 
effects. Furthermore, finite volume effects tend towards different directions in QCD 
and QQCD for decreasing M n . This means that quenching errors in £ may be signifi- 
cantly larger than what was estimated before. Therefore one has to be cautious when 
using the existing quenched lattice QCD results for £ in phenomenological work. In 
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Figure 4.12: (£b)fv in PQQCD plotted against m rr with L = 2.5 fm and two different values 
forM^. 

PQQCD, our results indicate that finite volume effects are typically between 3% and 
5% in the data range of future high-precision simulations, and they can be significantly 
amplified in the procedure of quark mass extrapolations. This means that they are not 
negligible in future lattice calculations of £. 
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Chapter 5 

Charge Radii of the Meson and 
Baryon Octets in QxPT and 
PQxPT 



In this chapter, we calculate the electric charge radii of the 577(3) pseudoscalar mesons 
as well as the 577(3) octet baryons in QxPT and PQxPT. The results are needed for 
the extrapolation of future lattice calculations of these observables. We also derive 
expressions for the nucleon and pion charge radii in 577(2) flavor away from the isospin 
limit. 

5.1 Introduction 

The study of hadronic electromagnetic form factors at low momentum transfer pro- 
vides important insight into the non-perturbative structure of QCD. Notable progress 
toward measuring the proton and neutron form factors has been made in recent years 
(see [96, 97] for references), including recent high precision measurements for the pro- 
ton [98]. Experimental study of the remaining octet baryons, however, is much harder. 
The charge radius of the X - has only recently been measured [99]. Although more 
experimental data for the other baryon electromagnetic observables can be expected 
in the future, progress will be slow as the experimental difficulties are significant. 
Theory, however, may have a chance to catch up. 

While quenched lattice calculations have already appeared [100-105], with the ad- 
vance of lattice gauge theory, we expect partially quenched calculations for many of 
these observables in the near future. Lattice simulations employing these approxima- 
tions need to be extrapolated from the heavier light quark masses used on the lattice 
(currently on the order of the strange quark mass) down to the physical light quark 
masses using the appropriate low-energy theories, QxPT and PQxPT (see Sections l2.3l 
andEl. 
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While there are a number of lattice calculations for observables such as the pion 
form factor [102, 106, 107] or the octet baryon magnetic moments [108, 109] that use 
the quenched approximation, there are currently no partially quenched simulations. 
However, given the recent progress that lattice gauge theory has made in the one- 
hadron sector and the prospect of simulations in the two-hadron sector [110-114], we 
expect to see partially quenched calculations of the electromagnetic form factors in the 
near future. 

This chapter is organized as follows. In Section IBT21 we calculate the charge radii 
of the meson and baryon octets in both QxPT and PQxPT up to next-to-leading (NLO) 
order in the chiral expansion. We use the heavy baryon formalism of Jenkins and 
Manohar [33,115], treat the decuplet baryons as dynamical degrees of freedom, and 
keep contributions to lowest order in the heavy baryon mass, Mb- These calculations 
are done in the isospin limit of SU(3) flavor. For completeness we also provide the 
PQ%PT result for the charge radii for the SU(2) chiral Lagrangian with non-degenerate 
quarks in Appendix|Dj In Section IB~3l we conclude. 

5.2 Charge Radii 

In this section we calculate the charge radii in PQxPT and QxPT. The basic con- 
ventions and notations for the mesons and baryons in QxPT and PQxPT have been 
laid forth in Chapter |2j they have also been extensively reviewed in the litera- 
ture [11,17-21,48]. 

5.2.1 Octet Meson Charge Radii 

The electromagnetic form factor Gx of an octet meson (j>x is required by Lorentz in- 
variance and gauge invariance to have the form 

{<Px{p')\Jg m \<Px(p)) = Gx(q 2 )(p + P'r (5.D 

where q^ = (p' —pY and p (p') is the momentum of the incoming (outgoing) meson. Con- 
servation of electric charge protects it from renormalization, hence at zero momentum 
transfer eGx(0) = Qx, where Qx is the charge of 4>x- The charge radius rx is related 
to the slope of Gx(q 2 ) at q 2 = 0, namely 

<r x >=6^G x (q 2 )\ g 2 =0 . (5.2) 

There are three terms in the 0(E 4 ) Lagrangian 

+ia 9 str(L fMU D^D u tf + R^D^D"T<) + ... (5.3) 
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Figure 5.1: Loop diagrams contributing to the octet meson charge radii in PQxPT. Octet 
mesons are denoted by a dashed line, singlets (hairpins) by a crossed dashed line, and the 
photon by a wiggly line. Only the third diagram has q 2 dependence and therefore contributes 
to the charge radius. 
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Figure 5.2: Wavefunction renormalization diagrams in PQxPT. These diagrams, together with 
the third one in Fig. 05.lt . ensure meson electric charge non-renormalization. 



that contribute to meson form factors at tree level. Here L^ v , i? M „ are the field-strength 
tensors of the external sources, which for an electromagnetic source are given by 



L,u, = iV = eQidf.Au - d u A„) + ie z Q z [A^A v 



(5.4) 



Unlike QxPT, where the low-energy constants are unique and have no known connec- 
tion to xPT, in PQxPT the parameters in ( I5.3I > are the dimensionless Gasser-Leutwyler 
coefficients of xPT [8] which can be seen by looking at mesons that contain sea quarks 
only. 

To calculate the charge radii to lowest order in the chiral expansion one has to 
include operators of C in ( 12. 341 ) to one-loop order [see Figs. ( 15. It and ( 15.21 11 and opera- 
tors of ( 15. 3D to tree level. Using dimensional regularization, where we have subtracted 
i + 1 - 7 + log 4tt, we find in PQxPT for the vr+ 



p Qr„2\ 



<V?(« 



—^^F^+F^+agj^q 2 , 



(5.5) 



which interestingly does not depend on the charges of the sea and ghost quarks, qj, qi, 
q r . For the K + we find 



G^ + (q 2 ) 



1 + 
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( 3 -Qjij F uu - 


( 3 -Qjij F uj ~ 


(!- 


Qr 1 r wr 




16vr 2 / 2 
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(5.6) 
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and for the K° we find 



G^tf) 



16vr 2 / 



2 fl 



g - ? 3 7 J Ku + I - + Qji J Kj + ( - + q r J Kr + I - + q r ) ^ 



3 - 9ji + Qr j F us - ( - + g j7 j F js - ( - + q r ) F rs 



(5.7) 



Here qji = qj + <# and we have defined 



Fqq> = ^ log ^21 - m% Q ,T ( -4- ) , (5.8) 

6 M \ m QQ'/ 



where the function .F(a) is given by 



^-U-sJV 1 -^ \_ i+fc+1 + s-3- <6 ' 9) 

The first derivative of Fqqi at q 2 = 0, needed to calculate the charge radii, becomes 

Zj2 F QQ> l? 2 =o = ^g -^- + 1. (5.10) 

Charge conjugation implies 

G£? = -<£?, g£=-G£?, and Cg = -G$, (5.11) 

which we have also verified at one-loop order. The form factors of the flavor diagonal 
mesons are zero by charge conjugation invariance. In the limit rrij — >■ fh, m r — >■ m s we 
recover the QCD result [8,9] as expected. 

It is interesting to note, that duplicating these calculations for QxPT shows that 
there is no meson mass dependence at this order. Specifically we find 

G%tf) = -Gj.(g 2 ) = G Q K+ (1 2 ) = -G Q K -{q 2 ) = 1 + j 2 ^q 2 , (5.12) 

and the form factors of the neutral mesons are zero. Here we annotate the quenched 
constant Og with a "Q" since its numerical value is different from the one in Eq. ( I5.3D . 
Eq. ( 15.121 ) reflects that flavor-singlet loops do not contribute to the ^-dependence at 
this order; thus the virtual quark loops are completely removed by their ghostly coun- 
terparts. This can readily be seen by considering the quenched limit of Eqs. J5.5D — 
( I5.7D . The meson mass independence reveals once again the pathologic nature of the 
quenched approximation and seriously puts into question xPT extrapolations to the 
physical pion mass. 
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5.2.2 Octet Baryon Charge Radii 

The electromagnetic form factors at or near zero momentum transfer that enable the 
extraction of the baryon magnetic moments and charge radii have been frequently in- 
vestigated in QCD [115-124]. There are also recent quenched and partially quenched 
calculations of the octet baryon magnetic moments in Q%PT and PQxPT [12, 13,31]. 
Here, we extend these calculations to the octet baryon charge radii. We retain spin-3/2 
baryons in intermediate states since formally A ~ m n . 

Using the heavy baryon formalism [33, 115], the baryon matrix element of the elec- 
tromagnetic current J M can be parametrized in terms of the Dirac and Pauli form 
factors Fi and F 2 , respectively, as 

(B(P) l^l B(p)) = u(p') {<AFi(<z 2 ) + ^p^ 2 (</ 2 )} u(p) (5.13) 

with q = p' — p. The Sachs electric and magnetic form factors defined as 

G E (q 2 ) = F 1 (q 2 ) + ^ J F 2 (q 2 ) (5.14) 

G M (q 2 ) = F 1 (q 2 ) + F 2 (q 2 ) (5.15) 

are particularly useful. The baryon charge Q, electric charge radius < r\ >, and 
magnetic moment p can be defined in terms of these form factors by 



and n = G M (0)-Q. (5.16) 



Q = G E (0), < r\ >= Q^G E (q 2 : 
Here the baryon charge Q is in units of e. 

Analysis in PQxPT 

Let us first consider the calculation of the octet baryon charge radii in PQxPT. The 
Lagrangian describing the relevant interactions of the Bijk and %jk with the pseudo- 
Goldstone mesons is 

C = 2a (BS^BAf,) + 2/3 (BS^A^B) + J^C [(T ' A V B) + h.c] . (5.17) 

The axial-vector and vector meson fields A^ and V^ are defined by analogy to those in 
QCD: 

A» = - fed"? - e+d^) and V = - Ud^ + ^d^A . (5.18) 

The vector 5 M is the covariant spin operator [33, 115, 125]. The constants a and /3 are 
easily calculated in terms of the constants D and F that are used for the 5C/(3) va i 
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analogs of these terms in QCD. Restricting the indices of Bijk to i,j,k = 1,2,3 one 

easily identifies 

2 5 

a = -D + 2F and j3 = — D + F. (5.19) 

The leading tree-level correction to the magnetic moments come from the dimen- 
sion-5 operators 1 



ie 



2M B 



[fi a (B[S„, S V \BQ) + up {B[S^ S V ]QB)] F^ (5.20) 



which can be matched on the QCD Lagrangian upon restricting the baryon field indices 
to 1-3 

C = — [ti D \x(B[Sv,S v ]{Q,B})+» F \x(B[Sp,S v \[Q,B])\F> a ' (5.21) 

where 

2 5 

A*q = oMD + 2//i? and/i/j = --/xd + /Uf (5.22) 

at tree level. The magnetic moments contribute the so-called Foldy term to charge 
radii via ^(0) in Eq. ( 15.14ft . Likewise, further leading tree-level corrections to the 
charge radii come from the dimension-6 operators 

C = -^ [c a (BBQ) + op {BQB)] v^F^ (5.23) 

x 

and the parameters c + and c_ , defined by 

2 5 

c a = -c + + 2c_ and c/3 = — c + + c-, (5.24) 

are the same as those used in QCD. Here, we take the chiral symmetry breaking scale 
A x ~ 4irf for the purpose of power counting. The NLO contributions arise from the 
one-loop diagrams shown in Figs. ( 15. 3D and 15. 4D . To calculate the charge radii we 
need the form factors F\ to first order in q 2 and ^(0) we find 

< r% > = -j T {Qc-+a D c + ) + ^j^(Qti F + a D tiD) 

X, -D 



16vr 2 / 



r 2f2 2^ 



2 2 

A x log ^ - 5 fa log ^f + 10 x G{m x , A, A*) 



A' 

(5.25) 



Here, we have defined the function Q(m, A, /j) by 



/» 



2 



A , A- VA 2 -m 2 + ie 



G(m, A, /x) = log — - log v . . (5.26) 

M V A 2 - m 2 A + V A 2 - m 2 + ie 
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Figure 5.3: Loop diagrams contributing to the baryon magnetic moments and charge radii. A 
thin (thick) solid line denotes an octet (decuplet) baryon. The last two diagrams in the first row 
contribute to both magnetic moments and charge radii; the magnetic moment part of which has 
already been calculated in [31]. The first diagram in row 1 contributes q 2 dependence only to F\ 
and therefore is relevant for the charge radii. The remaining diagrams have no q 2 dependence. 
These along with the wave function renormalization diagrams in Fig. d5.4t maintain charge 
non-renormalization. 



Table 5.1: Octet baryon tree-level contributions in QCD, QQCD, and PQQCD. 
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Note that in Eq. ( 15.251 1 the only loop contributions we keep are those non-analytic in 
m x . 

The parameters for the tree-level diagrams are listed in Table IBTT1 The computed 
values for the fix, fi' x , and Ax coefficients that appear in Eq. ( 15.251 1 are listed for the 
octet baryons in Tables 15.21 — 15.91 The corresponding values for the AS transition 
are given in Table l5~10l In each table we have listed the values corresponding to the 
loop meson that has mass mx ■ If a particular meson is not listed then the values for 
Px, P'x> an d Ax are zero. 2 



1 Here we use F M „ = d^A u - d u A^. 

2 We have defined the coefficients fix and fix to correspond to those defined in [31] where fi 
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Figure 5.4: Wave function renormalization diagrams needed to maintain baryon electric 
charge non-renormalization. 



Table 5.2: The coefficients f3 x , P' x , and A x in 577(3) flavor p QxPT for the proton. 
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Table 5.3: The coefficients (3 X ,(3' X , and A x in 577(3) flavor PQxPT for the neutron. 
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Analysis in QxPT 

The calculation of the charge radii can be easily executed for QxPT. The operators 
in Eqs. ( 15.211 1 and d5.23l > contribute, however, their low-energy coefficients cannot be 
matched onto QCD. Therefore we annotate them with a "Q". Additional terms involv- 
ing hairpins [11, 13] do not contribute as their contribution to the charge radii is of the 
form (^q/A^) logm-g and therefore of higher order in the chiral expansion. We find 
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(5.27) 



As with the meson case, the diagram where the photon couples to the closed meson 
loop does not contribute to baryon charge radii in the quenched case, cf., Ax is zero. 



aD hd + -rrfs J2x [Pxmx + &^( m x,A, A*)] an d the function f(mx,A, fi) is given in [31]. 
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Table 5.4: The coefficients j3 x , P' x , and A x in 577(3) flavor PQxPT for the ^ 
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Table 5.5: The coefficients p x , P x , and A x in SU(3) flavor PQxPT for the S°. 
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Table 5.6: The coefficients X , P' x , and A x in SU(3) flavor PQxPT for the S" 
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Table 5.7: The coefficients /3 X , 0' x , and A x in St/(3) flavor PQxPT for the 
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Table 5.8: The coefficients (i x , 0' x , 
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Table 5.9: The coefficients 0x, 0' x > anci A * 


- in SU(3) flavor PQxPT for the A. 
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Table 5.10: The coefficie nts /3 X ,/3' X , and A x in 577(3) flavor PQx PT for the AS transition. 
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Table 5.11: The coefficients X and X ' in 5C/(3) flavor Q X PT for the octet baryons. 
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The remaining coefficients appearing in Eq. ( 15.271 ) are listed in Table 15.111 and stem 
from meson loops formed solely from valence quarks. 



5.3 Conclusions 

We have calculated the charge radii for the octet mesons and baryons in the isospin 
limit of PQxPT and also derive the result for the nucleon doublet and pion triplet 
away from the isospin limit for the SU(2) chiral Lagrangian. For the octet mesons and 
baryons we have also calculated the Q^PT results. 

We find that new operators in the QQCD Lagrangian, which are non-existent in 
QCD, enter at NNLO. Hence, formally our NLO result is not more divergent than its 
QCD counterpart. This, however, does not mean that our result is free of quenching 
artifacts. While the expansions about the chiral limit for QCD and QQCD are formally 
similar, < r 2 >~ a + log rriQ + . . . , the QQCD result is anything but free of quenched 
oddities: for certain baryons, £ and H~ in particular, diagrams that have bosonic or 
fermionic mesons running in loops completely cancel so that = 0. In other words, 
< r 2 >~ a + . . . and the result is actually independent of ttiq! The same behavior 
is found for the charge radii of all mesons in QQCD as the meson loop contributions 
entirely cancel. 

PQQCD, on the other hand, is free of such eccentric behavior. The formal behavior 
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of the charge radius in the chiral limit has the same form as in QCD. Moreover, there 
is a well-defined connection to QCD and one can reliably extrapolate lattice results 
down to the quark masses of reality. The low-energy constants appearing in PQQCD 
are the same as those in QCD and by fitting them in PQxPT one can make predictions 
for QCD. Our PQxPT result will enable the proper extrapolation of PQQCD lattice 
simulations of the charge radii and we hope it encourages such simulations in the 
future. 
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Chapter 6 

Electromagnetic Properties of 
the Baryon Decuplet in QxPT 
and PQxPT 



In this chapter, we calculate electromagnetic properties of the decuplet baryons in 
QxPT and PQxPT. We work at NLO in the chiral expansion, LO in the heavy baryon 
expansion, and obtain expressions for the magnetic moments, charge radii, and elec- 
tric quadrupole moments. The quenched calculation is shown to be pathological since 
only quenched chiral singularities are present at this order. We present the partially 
quenched results for both the 577(2) and 577(3) flavor groups and use the isospin limit 
in the latter. These results are necessary for proper extrapolation of lattice calcula- 
tions of decuplet electromagnetic properties. 

6.1 Introduction 

Experiments measuring the decuplet magnetic moments are anticipated in the fore- 
seeable future. The Particle Data Group lists values for the A ++ magnetic moment [1] 
but with sizeable discrepancy and uncertainty, even among the two most recent re- 
sults [126, 127]. A report of the initial measurement of the A + magnetic moment [128] 
has recently appeared, and further data are eagerly awaited. 

Lattice calculations can be used to calculate the decuplet electromagnetic mo- 
ments. While lattice simulations using the quenched approximation have already ap- 
peared [129], there are currently no partially quenched simulations. However, we ex- 
pect to see partially quenched calculations of the decuplet electromagnetic form factors 
in the near future. 

Whatever approximation is used — quenched or partially quenched — , now and in 
the near future all these calculations are performed with unphysically large quark 
masses and therefore must be extrapolated down to the physical light quark masses. 
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For QQCD lattice calculations this extrapolation can be accomplished by the use of 
QxPT. However, the results are not only plagued by quenching artifacts but also 
unrelated to QCD, as explained in Chapter |2] In fact, several examples show that 
the behavior of meson loops near the chiral limit is frequently misrepresented in 
QxPT [2,4, 13,35, 130, 131]. We find this to also be true for the decuplet electromag- 
netic observables. Indeed, to the order we work only quenched chiral singularities 
are present in the quenched electromagnetic moments; the charge radii have no quark 
mass dependence at all. Of course, the unattractive features of QQCD can be remedied 
by using PQQCD and PQxPT. 

The chapter is organized as follows. First, in Section ET2l we calculate the electro- 
magnetic moments and charge radii of the decuplet baryons in both QxPT and PQxPT 
up to NLO in the chiral expansion. We use the heavy baryon formalism of Jenkins 
and Manohar [33, 115] and work to lowest order in the heavy baryon expansion. These 
calculations are done in the isospin limit of SU(3) flavor. Expressions for form fac- 
tors with the q 2 dependence at one loop are given in Appendix lE.il For completeness 
we also provide the PQxPT result for the baryon quartet electromagnetic moments 
and charge radii for the SU(2) chiral Lagrangian with non-degenerate quarks in Ap- 
pendix |K2| We conclude in Section ET3l 



6.2 Decuplet Electromagnetic Properties 

The electromagnetic moments of decuplet baryons in xPT have been investigated pre- 
viously in [132, 133]. Additionally there has been interest in the decuplet electromag- 
netic properties in the large N c limit of QCD [134-137]. In this section we calculate the 
decuplet electromagnetic moments and charge radii in PQxPT and QxPT. The basic 
conventions and notations for the mesons and baryons in PQxPT have been laid forth 
in the last Chapter|2j Additionally the decuplet charge radii in xPT are provided since 
they have not been calculated before. First we review the electromagnetic form factors 
of heavy spin-3/2 baryons. 

Using the heavy baryon formalism [33, 115], the decuplet matrix elements of the 
electromagnetic current J p can be parametrized as 

(T(p')l J P \T(p)) = -u»(p')O^Mp), (6.D 

where u fl (p) is a Rarita-Sch winger spinor for an on-shell heavy baryon satisfying 
v^u,j,{p) = and S^u^p) = 0. The tensor O ppv can be parametrized in terms of four 
independent, Lorentz invariant form factors 

^ ^ »" H<« 2 > + J^fW)} + Mr h^ + ^F^e 2 >} - 

(6.2) 
where the momentum transfer q = p' — p. The form factor F\(q 2 ) is normalized to the 
decuplet charge in units of e: 2*i(0) = Q. At NLO in the chiral expansion, the form 
factor G 2 (q 2 ) = 0. 
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Extraction of the form factors requires a nontrivial identity for on-shell Rarita- 
Schwinger spinors [138]. For heavy baryon spinors, the identity is 



u a (p'){q a 9^ ~ q P <T)up(p) = u a (p') 



a 2 I 

" „a0„.u i n„afl[c« out „ , 



-0 a V + 2g a P[SV, S v ]q u + —-q a q P v 



up(p). 



2Mb M b 

(6.3) 
Linear combinations of the above (Dirac- and Pauli-like) form factors make the electric 
charge GEo(q 2 ), magnetic dipole GMi(q 2 ), electric quadrupole GE2(q 2 ), and magnetic 
octupole GM3(q 2 ) form factors. This conversion from covariant vertex functions to mul- 
tipole form factors for spin-3/2 particles is explicated in [138]. For our calculations, the 
charge radius 



< r% >= 6—^GEoiq 2 ] 
the magnetic moment 



dFi(O) 1 

q 2 =0 



6 <^-T2Ml [2Q - 3F2 ^- G ^ (6A) 



fi = G M i(0)-Q = F 2 (0), (6.5) 

and the electric quadrupole moment 

Q = G E2 (0)-Q = ~G l (0). (6.6) 

To the order we work in the chiral expansion, the magnetic octupole moment is zero. 

6.2.1 Analysis in PQxPT 

Let us first consider the calculation of the decuplet baryon electromagnetic properties 
in PQxPT. Here, the leading tree-level contributions to the magnetic moments come 
from the dimension-5 operator 

C = ti c ^^Qr / )F ltv , (6.7) 

which matches onto the xPT operator [116] 



ieQi 
~Mb~ 



C = Vc^-T^F^, (6.8) 



when the indices in Eq. H6.7I ) are restricted to 1-3. Here Qi is the charge of the ith 
decuplet state. Additional tree-level contributions come from the leading dimension-6 
electric quadrupole operator 

C = -Qc^iT^QT^d^F^. (6.9) 

x 
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Here the action of ' . . . * on Lorentz indices produces the symmetric traceless part of 
the tensor, viz., O^ = 0» v + O"? - \g^ u O a a . The operator in Eq. CHJ matches onto 
the xPT operator [132] 



£ 



x 



(6.10) 



The final tree-level contributions arise from the leading dimension-6 charge radius 
operator 

(6.11) 



£ = c CJT (T°QT tT )v li d v F>»' 

x 



which matches onto the xPT operator 



>~- — <--c j 2 » o",i < J [i u v ± 

X 



(6.12) 



Notice that the PQQCD low-energy constants /i c , Q c , and c c have the same numerical 
values as in QCD. 

The NLO contributions to electromagnetic observables in the chiral expansion arise 
from the one-loop diagrams shown in Fig. 16.11 In addition to the terms describing the 
interactions of the Bijk and 7yj. with the pseudo-Goldstone mesons given in Eq. ( I5.17D 
we need the term 

C = Tti (T^S^A^) (6.13) 

where the axial-vector M l is defined in Eq. A5.18D . Calculation of the diagrams yields 



F 1 (q z )=Q[l 
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(6.15) 
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log- 
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(6.16) 



The only loop contributions kept in the above expressions are those non-analytic in 
the quark masses. The full q 2 dependence at one-loop has been omitted from the above 
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Figure 6.1: Loop diagrams contributing to the decuplet electromagnetic moments and charge 
radii. Octet mesons are denoted by a dashed line, singlets (hairpins) by a crossed dashed 
line, and the photon by a wiggly line. A thick (thin) solid line denotes a decuplet (octet) baryon. 
The diagrams in the first row contribute to the electromagnetic moments and charge radii. The 
remaining diagrams with a photon have no ^-dependence. These, along with the wavefunction 
renormalization diagrams, ensure non-renormalization of the electric charge. 
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Table 6.1: The coefficients A\ for SU(3) for each of the decuplet states in %PT and PQxPT. The index X corresponds to the 
loop meson that has mass m X - Here, we have used q 3 i = q 3 + qi. 
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Figure 6.2: The charge radii of the decuplet baryons in xPT. The contribution from countert- 
erms has been set to zero. The radii (squared) here come from the one-loop diagrams only and 
are plotted in units of fin 2 . For reference we have shown both the proton and neutron charge 
radii [1] (solid and dashed lines, respectively). 



expressions but is given in Appendix lE.il The function TZ(x) is given in Eq. tB.21 . The 
sum in the above expressions is over all possible loop mesons X. The computed values 
of the coefficients A x that appear above are listed in Table l6~Tl for each of the decuplet 
states T. In the table we have listed values corresponding to the loop meson that has 
mass mx for both ^PT and PQxPT. In particular, the %PT coefficients can be used 
to find the QCD decuplet charge radii, which have not been calculated before. Using 
Eq. d6.4l >. the expression for the charge radii is 
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(6.17) 



In the absence of experimental and lattice data for the low energy constants Q c and 
c c , we cannot ascertain the contributions to the charge radii from local counterterms in 
xPT. We can consider, however, just the formally dominant loop contributions. To this 
end, we choose the values C = —2D and H = —3D with D = 0.76 [33], and the masses 
A = 270 MeV, m^ = 140 MeV, and m K = 500 MeV. The loop contributions to the charge 
radii in ^PT are then evaluated for the decuplet at the scale fj, = 1 GeV and plotted in 
Fig. El 
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-^ — r- 

/ / 

-^ — r- 

Figure 6.3: Hairpin diagrams that give contributions of the form ~ y? log m q to decuplet electro- 
magnetic moments in QxPT. The square at the photon vertex represents the relevant operators 
from Eqs. 0E3, EU, and <6~19l . 

6.2.2 Analysis in QxPT 

The calculation of the charge radii and electromagnetic moments can be correspond- 
ingly executed for QxPT. The operators in Eqs. ( 16.71 1. ( 16.91 1 and ( 16.111 1 contribute, how- 
ever, their low-energy coefficients cannot be matched onto QCD. Therefore we anno- 
tate them with a "Q". The loop contributions encountered in xPT and PQxPT above 
no longer contribute because in QxPT A^ = for all decuplet states T. This can be 
readily seen in two ways. The quenched limit 1 of the coefficients in Table 16.11 makes 
immediate the vanishing of Ajr . Alternately one can consider the relevant quark flow 
diagrams with only valence quarks in loops. Due to the symmetric nature of T^ k , these 
loops are completely canceled by their ghostly counterparts. For the charge radii, there 
are no additional diagrams to consider at this order from singlet hairpin interactions. 
Thus in QxPT, the charge radii have the form 

(6.18) 

where the dependence on the quark mass enters at the next order. 

Additional terms of the form ^ logm g involving hairpins [11, 13] do contribute to 
the electromagnetic moments as they are of the same order in the chiral expansion. 
As explained in [139], the axial hairpins do not contribute. In the diagrams shown 
in Fig. 16.31 one sees that there are contributions from the electromagnetic moments of 
the decuplet and octet baryons as well as their transition moments. These interactions 
are described by the operators in Eqs. ( 16.71 ) and H6.9I ) (now with quenched coefficients) 




1 In this case, the quenched limit simply means to remove sea quarks and to fix the charges of the 
ghost quarks to equal those of their light quark counterparts. 
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It is easier to work with the combinations ^ and fij, defined by 



K 



fi% + 2fi% and $ = --1% + n% 



u ^L ± VOf 

(6.19) 



(6.20) 



To calculate the quenched electromagnetic moments, we also need the hairpin 
wavefunction renormalization diagrams shown in Fig. 16.11 These along with the di- 
agrams in Fig. 16.31 are economically expressed in terms of the function 

Y(mi, Ai, n) + Y(m 2 , A 2 , n) - Y(mi,A 2 , n) - Y(m 2 , Ai, fi) 
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where 



y(m, A, //) = A I m 2 - ^A 2 ) log ^- + \m{A 2 - m 2 )K ( — ) 
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(6.22) 



and we have kept only non-analytic contributions. The following shorthands are con- 
venient 

7 %V = J ( m %' m V ' '^)' 



VqV q > 
rAA 
VqV q > 



I(m Vq ,m Vq ,,A,0,n), 
I(m Vq ,m v A,A,{i). 



(6.23) 



Specific limits of the function I Vqr) , appear in [13]. The wavefunction renormalization 
factors arising from hairpin diagrams can then be expressed as 
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(6.24) 



The coefficients B xx i and C xx > are listed in Table l6~2l The sum in Eq. ( 16.24ft is over 

XX' = rj u rj u , rj u rj s , f] s r] s . Combining these factors with the tree-level contributions and 
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Table 6.2: The SU(3) coefficients B xx > and C xx , in QxPT. 
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Table 6.3: The SU(3) coefficients D X x> in QxPT. Decuplet states not listed have D xx , = 0. 
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including the corrections from the diagrams in Fig. 16.31 we arrive at the quenched 
decuplet magnetic moment 
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(6.25) 



and the quenched electric quadrupole moment 
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(6.26) 

In Eq. H6.25H the required values for the constant ao are: ao = \ for £*' + , S*'°, £*>~, 
and H*~, and ao = — | for H*'°. The coefficients D xx > are listed in Table IQ1 If a 
particular decuplet state is not listed, the value of D XX i is zero for all singlet pairs 
XX'. 

The above expressions can be used to properly extrapolate quenched lattice data 
to the physical pion mass. For example, the expression for the quenched magnetic 
moments for the A baryons [Eq. d6.25l ll reduces to 

2 «"?( 1 ^i^¥ (w0,2/ ")- 



/' 



(6.27) 
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In the above expression we need 

fi z A 2 — m x \mx ) 

where the ellipses denote terms analytic in mx ■ Utilizing a least squares analysis and 
using the values [x = 1 GeV and A = 270 MeV , we extrapolate the quenched lattice 
data [129] to the physical pion mass and find for the A resonances 

fj, = 2.89 Q [mat]- (6.29) 

This is in contrast to the value [i « 2.49 Q [/xjv] found from carrying out a xPT-type 
extrapolation [140]. Notice that for many of the decuplet states, in particular the 
A baryons, the quenched magnetic moments and electric quadrupole moments are 
proportional to the charge Q (unlike the xPT and PQxPT results). This elucidates the 
trends seen in the quenched lattice data [129]. 

6.3 Conclusions 

We have calculated the electromagnetic moments and charge radii for the SU(3) decu- 
plet baryons in the isospin limit of PQxPT and also derived the result for the baryon 
quartet away from the isospin limit for the SU(2) chiral Lagrangian. The q 2 depen- 
dence of decuplet form factors at one-loop appears in Appendix IE.1I We have also 
calculated the QxPT results. 

For the decuplet baryons' electromagnetic moments and charge radii, the opera- 
tors, which are included in the QxPT but not in the xPT Lagrangian, enter at NNLO 
in the QxPT expansion. Hence, formally our NLO result is not more divergent than 
its QCD counterpart. This, however, does not mean that our result is free of quench- 
ing artifacts. While the expansions about the chiral limit for QCD and QQCD charge 
radii are formally similar, < r 2 >~ a + filogrriQ + . . . , the QQCD result consists en- 
tirely of quenched oddities: for all decuplet baryons, diagrams that have bosonic or 
fermionic mesons running in loops completely cancel so that j3 = 0. In other words, 
the quenched decuplet charge radii have the behavior < r 2 >~ a + . . . and the result 
is actually independent of itiq at this order. For the quenched decuplet magnetic mo- 
ments and electric quadrupole moments, expansions about the chiral limit are again 
formally similar: \± ~ a + (3 log itiq + ^Jm~Q + . . . and Q ~ a + (5 log itiq + . . . . How- 
ever, quenching forces 7 = and both /3's arise from singlet contributions involving the 
parameter /j, , which is of course absent in QCD. Thus the leading non-analytic quark 
mass dependence that remains for these observables is entirely a quenched peculiarity. 

PQQCD, on the other hand, is free of such excentric behavior and our PQxPT result 
will enable the proper extrapolation of PQQCD lattice simulations of the decuplet elec- 
tromagnetic moments and charge radii and we hope it encourages such simulations in 
the future. 
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Chapter 7 

Baryon Decuplet to Octet 
Electromagnetic Transitions in 
QxPT and PQxPT 



In this chapter, we calculate baryon decuplet to octet electromagnetic transition form 
factors in quenched and partially quenched chiral perturbation theory. Again, we work 
in the isospin limit of SU(3) flavor, up to NLO in the chiral expansion, and to leading 
order in the heavy baryon expansion. Our results are necessary for proper extrapola- 
tion of lattice calculations of these transitions. We also derive expressions for the case 
of SU{2) flavor away from the isospin limit. 

7.1 Introduction 

The study of the baryon decuplet to octet electromagnetic transitions provides impor- 
tant insight into the strongly interacting regime of QCD. Spin-parity selection rules 
for these transitions allow for magnetic dipole (Ml), electric quadrupole (E2), and 
Coulumb quadrupole (C2) amplitudes. Understanding these amplitudes, both in the- 
ory and experiment, gives insight into the ground state wavefunctions of the lowest ly- 
ing baryons. For example, in the transition of the A(1232) to the nucleon, if both baryon 
wavefunctions are spherically symmetric then the E2 and C2 amplitudes vanish. Ex- 
perimentally, Ml is seen to be the dominant amplitude. However, recent experimental 
measurements of the quadrupole amplitudes in the A — ► Nj transition [141, 142] show 
that the quadrupole amplitudes E2 and C2 are likely non-zero. This has revitalized 
the discussion as to the mechanism for deformation of the baryons. Although we ex- 
pect more experimental data in the future, progress will be slower for the remaining 
transitions as the experimental difficulties are large. 

First-principle lattice QCD calculations of these matrix elements can provide a the- 
oretical explanation of these experimental results. In fact, the experimental difficul- 
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ties may force us to rely on lattice data for the non-nucleonic transitions. Recently 
several such lattice calculations [143, 144], which improve upon an earlier one [145], 
have appeared. Unfortunately now and foreseeably, these lattice calculations cannot 
be performed with the physical masses of the light quarks and must be extrapolated 
to the physical light quark masses. 

For future lattice calculations that use the partially quenched approximation of 
QCD one needs to fit PQxPT to the lattice data in order to determine the low-energy 
constants and to actually make physical predictions for QCD. Unfortunately, partially 
quenched simulations for the A — » N^y transition do not exist yet; what does exist 
are simulations [143, 144] that use the quenched approximation. These simulations 
must be fit using QxPT that exhibit the sickness outlined in the previous chapters. 
However, we expect to see partially quenched calculations of these form factors in the 
near future. 

This chapter is organized as follows. In Section 17.21 we calculate baryon decuplet 
to octet transition form factors in both QxPT and PQxPT up to NLO in the chiral 
expansion and keep contributions to lowest order in the heavy baryon mass, Mg. These 
calculations are done in the isospin limit of SU(3) flavor. For completeness we also 
provide the PQxPT results for the transitions using the SU(2) chiral Lagrangian with 
non-degenerate quarks in Appendix \F\ In Section 17131 we conclude. 

7.2 Baryon Decuplet to Octet Transition 

The electromagnetic baryon decuplet to octet transitions have been investigated pre- 
viously in xPT [146-149]. Very recently there also has been renewed interest in these 
transitions in the large N c limit of QCD [150]. Here we calculate these transitions in 
PQxPT and QxPT. While we have reviewed PQxPT briefly in the last section and our 
recent papers [4,5], for QxPT we refer the reader to the literature [11, 17-21,48]. 

Using the heavy baryon formalism [33, 115], transition matrix elements of the elec- 
tromagnetic current J p between a decuplet baryon with momentum p' and an octet 
baryon with momentum p can be parametrized as 

(B(p)\JP\T(p')) = u{p)0^{p>), (7.1) 

where u fl (p) is a Rarita-Sch winger spinor for an on-shell decuplet baryon satisfying 
v^Ufj,{p) = and S fJ, u ll (p) = 0. The tensor O pp can be parametrized in terms of three 
independent, Lorentz invariant, dimensionless form factors [151] 



G -fV*r-^) + |g; 



Gsiq 2 



+ J^i(Q 2 9^-Q fl Q P )S-q, (7.2) 

4M^A v ' 
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where the momentum of the outgoing photon is q = p' — p. Here we have adopted the 
normalization of the G^{q 2 ) form factor used in [149] so that the leading contributions 
to all three form factors are of order unity in the power counting. 

Linear combinations of the above form factors at q 2 = make the magnetic dipole, 
electric quadrupole, and Coulombic quadrupole moments, 

°"« (0) = (I-<4) G ' (0) + l2^° 2<0> - 

Analysis in PQxPT 

Let us first consider the transition form factors in PQ^PT. Here, the leading tree-level 
contributions to the transition moments come from the dimension-5 and dimension-6 
operators 

C = \fl^ T 2M^ ( Bsfl Q TU ) F ^ + \/f QT ^ (^S^) v a d,F ua (7.4) 

X. 

where the action of ' . . . ^ on Lorentz indices produces the symmetric traceless part of 
the tensor, viz., 0^> = O^ + O^ - \g^ v O a a . Here the PQQCD low-energy constants 
Ht and Qt have the same numerical values as in QCD. 

The NLO contributions in the chiral expansion arise from the one-loop diagrams 
shown in Figs. ( I7.lt and d7.2D . However, because of the constraints satisfied by the 
on-shell Rarita-Schwinger spinors, the diagrams in Fig. d7.lt are all identically zero. 
For the calculation of the diagrams in Fig. ( 17. 2D we need the terms in the Lagrangian 
describing the interaction of the B^ and 7^ with the pseudo-Goldstone mesons given 
in Eas. d5~T7t and dB~T3t : 

C = 2q (BS^BA^ + 2/3 (BS^A^B) + J^C [(T U A U B) + h.c] + 27i (T'S^A^) . (7.5) 



We find 
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fi z V m x 



(7.6) 
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Figure 7.1: Loop diagrams that contribute to the transition moments but are zero to the order 
we are working. A thin (thick) solid line denotes an octet (decuplet) baryon whereas a dashed 
line denotes a meson. 



1 



1. 



Figure 7.2: Loop diagrams contributing to the transition moments 
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Table 7.1: Tree-level coefficients a T in X PT, QxPT, and PQxPT. 
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(7.8) 



where the function 1Z(x) is defined in Eq. HB.2H and we have only kept loop contribu- 
tions that are non-analytic in the meson mass mx- The tree-level coefficients ax are 
listed in Table 17711 and the coefficients for the loop diagrams in Fig. ( 17.2ft . 0% and fix, 
are given in Tables 17.21 and 17.31 respectively. In these tables we have listed values 
corresponding to the loop meson with mass mx- As required, in the QCD limit the 
PQxPT coefficients reduce to those of xPT. It is comforting that the one-loop results 
for the Gs(q 2 ) form factor are finite. This is consistent with the fact that one cannot 
write down a dimension-7 operator that contributes at the same order in the chiral 
expansion as our one-loop result for G^{q 2 ). The full one-loop q 2 dependence of these 
form factors can easily be recovered by replacing 



mx 
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x)q' L 



(7.9) 



Notice that the tree-level transitions £*' - — > £^7 and H*'~ — > H~7 are zero because 
they are forbidden by d <-► s [/-spin symmetry [152]. There is also symmetry between 
the S*' + — >■ £+7 and E*'° — > H°7 transitions as well as the E*~ — >■ £~7 and E* ~ — ► H~7 
transitions that holds to NLO in xPT and PQxPT. 



Analysis in QxPT 

The calculation of the transition moments can be repeated in QxPT. At tree level, the 
operators in Eq. ( 17.4ft contribute, but their low-energy coefficients cannot be matched 
onto QCD. Therefore we annotate them with a "Q". At the next order in the chiral ex- 
pansion, there are again contributions from the loop diagrams in Fig. ( 17.2ft . The results 
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Table 7.2: The SU(3) coefficients 0% in X PT and PQ X PT. 
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Table 7.3: The 577(3) coefficients 0% in X PT and PQ X PT. 
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Table 7.4: The SJ7(3) coefficients /3^ Q and /?£ Q in QxPT. 







4 


Q 




Px 


g 








7T 


K 


■K 




x 




A- 


>iV7 


1 
\/3 





U DQ - 


3F^) 







s *,+ . 


-> S+7 





1 

\/3 







-T^(^ Q - 


-3F«) 


s *,o . 


^£° 7 





2^3 







^(^ 9 - 


-3F«) 


s *,o 


-> A7 


1 

3 


~~ 6 


-i (J DS- 


3F«) 


~^ Q - 


3F«) 


s*>- - 


^s- 7 


















"*,0 _ 


.S°7 





1 

V3 







"^(^- 


-3F«) 


"*;~ 


-> 5-7 



















81 



,-X 



T 




,<X 



n— 



• -X- 






Figure 7.3: Loop diagrams contributing to the transition form factors in QxPT. The four 
diagrams correspond to terms involving the parameters A X x>, B X x>, Cxx<, and D X x< in Eqs. 
f7~13l and f7~14l 

are the same as in the partially quenched theory, Eqs. ( I7.6D — < T7THI >. with the coefficients 
0% and /3| replaced by f3^ Q and j3^' Q /(D Q - F Q ), which are listed in Table 17^1 

In addition, there are contributions of the form $ log m q at the same order in the 
chiral expansion that are artifacts of quenching. These come from hairpin wavefunc- 
tion renormalization diagrams and from the four loop diagrams in Fig. A7.3I I. In these 
diagrams the photon can couple to the baryon line via 
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(7.10) 



and via the terms in Eq. d7.4D including their hermitian conjugates (with quenched 
coefficients). 1 It is easier to work with the combinations /x'g and fij defined by 



Mc 



,Q 



,Q 



-IJl + 21% and ^ 
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(7.12) 



Although the argument presented in [139] does not apply to the case of different initial 
and final states, the axial hairpin interactions still do not contribute simply because 
their presence requires closed quark loops. The hairpin wavefunction renormalization 
diagrams have been calculated in QxPT for the baryon octet [13] (Z^) and decuplet [5] 
(ZjO and we do not reproduce them here. We find the hairpin contributions to the 



1 Note that possible contributions from diagrams involving 
are identically zero due to the constraints satisfied by the on-shell Rarita-Schwinger spinors. 



+ $(BQB)\ tfcflUF*"' + c?^(T"Q%)v^F^ 



(7.11) 
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transition form factors to be 
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(7.14) 



and Gf p (g 2 ) = 0. Thus in Q X PT: G9(q 2 ) = Gf Q {q 2 )+Gf p (q 2 ), where the /3 X and f3 x co- 
efficients of G- {q 2 ), Eqs. I7.6M7.8I . are understood to be replaced by their quenched 
values (3 X and [3 X /(D® — F Q ). Above we have used the shorthand notation I nn , = 
I(m Vq ,mr, ql , 0, 0, n), I^ v = I(mr, q ,m Vq , , A, 0, /x), and I Vq \j q , 2 = I(m Vv ,m %l , Ai, A 2 ,/x) for 
the function I (mi, m 2 , Ai, A 2 , /x) that is given by 

y(mi, Ai, /x) + Y(m 2 , A 2 , /i) - y (mi, A 2 , /j) - y(m 2 , Ai, /x) 



I(mi,m 2 ,Ai,A 2 ,/x) 
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(7.15) 
(7.16) 



'^ 

The coefficients Axx'j B xx i, C xx >, and Dxx' are listed in Tables I7TB1 and I7TB1 Note 
that the symmetry between the £*'+ — » £ + 7 and S*'° — » H°7 transitions as well as the 
S*' — ► S 7 and S*'~ — > H _ 7 transitions that holds in xPT and PQxPT is now broken 
by singlet loop contributions. 



7.3 Conclusions 

In this chapter we have calculated the baryon octet to decuplet transition form factors 
in QxPT and PQxPT using the the isospin limit of SU(3) flavor and have also derived 
the result for the nucleon doublet in two flavor PQxPT away from the isospin limit. 
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Table 7.5: The 


SU(3) coefficients A X x* and B^x' in 
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Table 7.6: The 


SU(3) 


coefficients Cxx> 


and L>xx' in QxPT. 
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For the decuplet to octet transition form factors our NLO QxPT results are not 
more divergent than their %PT counterparts: G\, G<i ~ a + j3 log tuq and G3 ~ a. This, 
however, does not mean that this result is free of quenching artifacts. The quenched 
transition moments pick up contributions from hairpin loops. A particular oddity is 
that the quark mass dependence of the £*' and E*~ quenched transition moments is 
solely due to the singlet parameter ^l; even worse, G®(q 2 ) = at this order. These 
transitions thus present extremes of the quenched approximation in agreement with 
the quenched lattice data of [145] where the S*' and E*>~ E2 moments were found to 
be significantly different from the other transitions. In contrast to QxPT results, our 
PQxPT results will enable not only the extrapolation of PQQCD lattice simulations of 
the transition moments but also the extraction of predictions for the real world: QCD. 
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Chapter 8 

Hadronic Electromagnetic 
Properties at Finite Lattice 
Spacing 



In this chapter we augment the electromagnetic properties of the octet mesons as well 
as the octet and decuplet baryons calculated in Chapters [5j-|2l in QxPT and PQxPT to 
include 0(a) corrections due to lattice discretization. We present the results for the 
SU(3) flavor group in the isospin limit as well as the results for SU(2) flavor with 
non-degenerate quarks. These corrections will be useful for extrapolation of lattice 
calculations using Wilson valence and sea quarks, as well as calculations using Wilson 
sea quarks and Ginsparg- Wilson valence quarks. 

8.1 Introduction 

In the previous three chapters, we considered the electromagnetic properties of the 
octet mesons and both the octet and decuplet baryons in QxPT and PQxPT. Owing 
in part to the charge neutrality of singlet fields, the quenched results are not more 
singular in the chiral limit than their unquenched counterparts. We showed, however, 
that despite this similarity, the quenched results contain singlet contributions that 
have no analog in xPT. Moreover, quenching closed quark loops alters the contribu- 
tion from chiral logs. For the decuplet baryon form factors, for example, quenching 
completely removes these chiral logs. Many others have also observed that the be- 
havior of meson loops near the chiral limit is misrepresented in QxPT, see for exam- 
ple [2, 13,47, 130, 131]. On the other hand, PQxPT results are devoid of such complica- 
tions and allow for a smooth limit to QCD. 

Not only are lattice calculations limited to unphysically large quark masses, they 
are also severely restricted by two further parameters: the size L of the lattice, that 
is not considerably larger than the system under investigation; and the lattice spac- 
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ing a, that is not considerably smaller than the relevant hadronic distance scale. To 
address the issue of finite lattice spacing, x?T has been extended (following the ear- 
lier work of [153-156]) in the meson sector to 0(a) for the Wilson action [157] and for 
mixed actions [158]. Corrections at 0(a 2 ) have also been pursued [159, 160]. Correc- 
tions to baryon observables have only recently been investigated [161]. To consider 
finite lattice spacing corrections, one must formulate the underlying lattice theory and 
match the new operators that appear onto those in the chiral effective theory. This 
can be done by utilizing a dual expansion in quark mass and lattice spacing. Follow- 
ing [159, 161], we assume a hierarchy of energy scales 

m q <C A x < - (8.1) 

a 

and ignore finite volume effects. The small dimensionless expansion parameters are 

{m q /A x , 
aA x , (8.2) 

where p is an external momentum. Thus we have a systematic way to calculate 0(a) 
corrections in %PT for the observables of interest. 

In this chapter we investigate the 0(a) corrections to the electromagnetic proper- 
ties of the meson and baryon octets, the baryon decuplet, and the decuplet to octet 
electromagnetic transitions in QxPT and PQxPT. We work up to NLO in the chiral ex- 
pansion and to leading order in the heavy baryon expansion. The paper is structured 
as follows. First, in Section l8~2l we review PQxPT at finite lattice spacing with mixed 
actions. Since the setup for QxPT parallels that of PQxPT, we will only highlight 
differences where appropriate. Next in Section 18.31 we calculate finite lattice spacing 
corrections to the charge radii of the octet mesons to 0(e 2 ). This is followed by the 
calculation of such corrections to: the charge radii and magnetic moments of the octet 
baryons; the charge radii, magnetic moments, and electric quadrupole moments of the 
decuplet baryons; and the decuplet to octet electromagnetic transition moments (Sec- 
tions UlU-EHi. Corresponding results for the above electromagnetic observables in the 
SU(2) flavor group are presented in Appendix lG.il In Appendix IG.2I we determine the 
0(a) corrections in an alternative power counting scheme for coarser lattices where 
e ~ aA x . A conclusion appears in Section 18771 

8.2 PQxPT at Finite Lattice Spacing 

In PQQCD [22-29] the quark part of the Symanzik Lagrangian [153, 162, 163] to 0{a) 
is written as 

C = Q(Hp-m Q )Q + aQa tMU G til/ CQQ, (8.3) 
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where the second term, the Pauli-term, breaks chiral symmetry in the same way 
as the quark mass term. Here, the nine quarks of PQQCD are in the fundamen- 
tal representation of the graded group 577(6|3) [14-16] and appear in the vector 
Q = (u,d,s,j,l,r,u,d,s) that obeys the graded equal-time commutation relation in 
Eq. ( 12.30ft . The quark mass matrix mq is given in Eq. ( 12.451 ) while the Sheikholeslami- 
Wohlert (SW) [164] coefficient matrix for mixed actions reads 

cQ = diag(cvvv s ,c s ,c s , c vv u ). (8.4) 

If the quark Qi is a Wilson fermion [10], then (cq)j = c sw . Alternately, if Qi is of the 
Ginsparg- Wilson variety [165] (e.g., Kaplan fermions [166] or overlap fermions [167]), 
then {cq)i = 0. Since one expects simulations to be performed with valence quarks 
that are all of the same species as well as sea quarks all of the same species, we have 
labeled the SW coefficients in Eq. ( 18.4ft by valence (v) and sea (s) instead of flavor. In 
the limit rrij = m u , mi = m<j, and m r = m s one recovers QCD at O(a). 

In addition to the SW term in Eq. ( 18.3ft . the vector-current operator of PQQCD also 
receives corrections at 0(a). There are three operator structures to consider [168] 

O% = aQQc m Q ^Q 

0% = aD u (QQc 2 ^ u Q), (8.5) 

where D M = D^ — D M and D M is the gauge covariant derivative. The form of the 
matrices Co, c\, and C2 in PQQCD is 

Cj — Qiag yCj , Cj , Cj , Cj , Cj , Cj , Cj , Cj , Cj J , I.O.OJ 

where cj and A are the coefficients of the vector-current correction operator O^ for 
valence and sea quarks respectively. If the vector-current operator is 0(a) improved in 
the valence (sea) sector, then cj = ( cj = ). The operator Oq, which corresponds to 
a renormalization of the vector current, contains a factor ofarriQ that renders it 0(e 4 ). 
Thus contributions to electromagnetic observables from Oq are neglected below. The 
equations of motion which follow from Eq. ( 18.3ft can be used to show that the operator 
02 is redundant up to 0(a 2 ) corrections. Therefore, we need not consider O^ ■ For ease 
we define the matrix product ci : g = Qc\. 

8.2.1 Mesons 

For massless quarks at zero lattice spacing, the Lagrangian in Eq. ( 18.3ft exhibits a 
graded symmetry SU (6|3) x,® SU (6|3) r®U (1)v that is assumed to be spontaneously bro- 
ken down to SU(6\3)v <8> U{l)y. The low-energy effective theory of PQQCD that results 
from expanding about the physical vacuum state is PQxPT. The emerging 80 pseudo- 
Goldstone mesons can be described at 0(e 2 ) by a Lagrangian which accounts now for 
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the two sources of explicit chiral symmetry breaking [153, 157, 158] 

£ = ! — str [D^Y)D^j + A m str (m Q £ + m^E f ) + 0^*0^*0 - Mo$o 

+aA a str(c Q £ + c^£ t ) (8.7) 

where E, $, M, M, and x are defined in Eqs. d2Tl5b . CTi . ( l2~48l l. and <2~49l Ex- 
panding the Lagrangian in Eq. < I8.7D one finds that to lowest order mesons with quark 
content QQ' have mass 

m QQ' = 72 [ X m(m Q + mcy) + a\ a (c Q + c Q >)] . (8.8) 

The flavor singlet field is $0 = str($)/y / 6. It is rendered heavy by the U(1)a 
anomaly and can be integrated out in PQxPT, with its mass /^o taken on the order 
of the chiral symmetry breaking scale, ^0 —* A. x . In this limit the propagator of the 
flavor singlet field is independent of the coupling a and deviates from a simple pole 
form [22,23]. In QxPT, the singlet must be retained. 

8.2.2 Baryons 

At leading order in the heavy baryon expansion and at 0(a), the free Lagrangian 
for the 240-dimensional super-multiplet Bijk and the 138-dimensional super-multiplet 
T* h fields is given by [11, 161] 

C = i(Bv- VB) + 2a M (BBM+) + 2/? M (jBM+B) + 2a M (BB) str (M+) 
+2a A (BBA+) + 2(3 A (BA+B) + 2a A (BB) str (.4+) 
-i (7% • PT M ) + A (7%) + 2 1M ^M+%) - 2a M (7%) str (M+) 
+2 1A (T^A+T^) - 2a A (T%) str (A + ) , (8.9) 

where M + = \ (^ m Q^ + C m Q^,) and A + = \a {^cq^ + £,cq£). The parenthesis nota- 
tion used in Eq. ( 18.9ft is defined in Eq. ( 12.62ft . Notice that the presence of the chiral 
symmetry breaking SW operator in Eq. ( 18.3ft has lead to new 0{a) operators (and new 
constants a A , Pa, &a, 7a, and a a) m Eq. ( 18.9ft . The Lagrangian describing the interac- 
tions of the B^k and T? k with the pseudo-Goldstone mesons is given in Eq. ( 17.5ft . with 
the axial-vector and vector meson fields A^ and V^ defined in Eq. ( 15.18ft . 

8.3 Octet Meson Properties 

The electromagnetic form factor G(q 2 ) of an octet meson <f> has the form 

W)\J»\m)=G(q 2 )(p + p'T (8.10) 
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where q^ = (p' — pY. At zero momentum transfer G(0) = Q, where Q is the charge of 
<fi. The charge radius r is related to the slope of G(q 2 ) at q 2 = 0, namely 



<- 2 >=4 Gtf: 



(8.11) 



Recall, at one-loop order in the chiral expansion the charge radii are 0(e 2 ) (see Chap- 
terUDl. 

There are two finite-a terms in the 0(e 4 ) Lagrangian [159] 






8 ^str( J D M St J D^)str(c S + c^) + a A ^ a 

(8.12) 
that contribute to meson form factors at tree level. The new parameters a AA and a A ^ 
in Eq. ( 18.121 1 are finite lattice spacing analogues of the dimensionless Gasser-Leutwyler 
coefficients 04 and 05 of xPT [8]. The above terms contribute to meson form factors at 
0(e 2 ) but their contributions are independent of q 2 and annihilated by the correspond- 
ing wavefunction renormalization, thus ensuring charge non-renormalization. 

The SW term can potentially contribute at 0(e 2 ) when A+ is inserted into the ki- 
netic term of the leading-order C in Eq. ( I8.7D . Contributions to form factors from such 
terms vanish by charge non-renormalization. Insertions of A+ into the a$ term of the 
Gasser-Leutwyler Lagrangian produces the 0(e 6 ) terms 



C = im^F^ str ({Q + ,A + }D^D u T) + {Q + ,A + }D^ D U Z^ 
+im 2 A x F^ str (q+D^A+D^ + S + CEU + D"s) 
+im 3 A x F^ str (q+D^D^ + Q+D^ D v t) str(^+), (8.13) 



where we have defined Q + = \ {$0$ + £Q£). These terms contribute at 0(e 4 ) to the 
charge radii and can be ignored (see Appendix IG.2I for discussion relating to larger 
lattice spacings). 

Additionally we must consider the contribution from the vector-current correction 
operator Of in Eq. ( I8.5D . In the meson sector, the leading operators Of in the effective 
theory can be ascertained by inserting aA x c\ t Q in place of Q in the operators that 
contribute to form factors. The effective field theory operators must also preserve the 
charge of the meson 4>. It is easiest to embed the operators Of in a Lagrangian so that 
electromagnetic gauge invariance is manifest. To leading order, the contribution from 
Of is contained in the term 

C = ia Ai9 aA x F^ str (ci jS ^S^S t + c liQ d^d u z) . (8.14) 

Thus the correction to meson form factors from Of is at 0(e 4 ). 
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The charge radius of the meson to 0(e) then reads 



2 24Q 1 v-^ „ , mi 

< r 2 >= a 9 jf + ^-^^^bg-f , (8.15) 



where X corresponds to loop mesons having mass nix [the masses implicitly include 
the finite lattice spacing corrections given in Eq. d8.8D . otherwise the expression is 
identical to the a = result]. The coefficients Ax in PQxPT appear in Chapter |5j In the 
case of QxPT, the coefficients Ax = for all loop mesons and there are no additional 
contributions from the singlet field at this order. Thus there is neither quark mass 
dependence nor lattice spacing dependence in the quenched meson charge radii at this 
order. 

8.4 Octet Baryon Properties 

Baryon matrix elements of the electromagnetic current J^ can be parametrized in 
terms of the Dirac and Pauli form factors F\ and F%, respectively, as outlined in Chap- 
ter[5l Recall, that the one-loop contributions in the chiral expansion to the charge radii 
are 0(e 2 ), while those to the magnetic moments are 0(e). 

There are no finite-a operators in Eq. d8.9l > that contribute to octet baryon form 
factors. As in the meson sector, however, the SW term could contribute when A+ is 
inserted into the Lagrangian. Here and henceforth we do not consider these insertions 
into the kinetic terms in Eq. ( I8.9D because their contributions alter the baryon charges 
and will be canceled by the appropriate wavefunction renormalization. 

The SW term, however, does contribute when A+ is inserted into the charge radius 
and magnetic moment terms. For the charge radius, we then have the 0(a) terms 

C = ^{-)^ + ^)^ + ^B kji {Q + ,A+} W B i: > k 'v ll d v F^ 

+^B kji {Q + ,A+} ii 'B* jk v„d v F» u + ^(-)*i'^ + ^B k i i & 

Ax 

+ *L(_)W+ 1 B fc * (Qj'Uf +A%Q^ ) B i ' j ' k v lx d u F^ 

Ax V ' ' / ■ 

+-±- [be (BBQ+) + b 7 (BQ+B)] v^F^ str(A+) + ^- (BB) v^F^ str(Q + A + ), 

X, X. 

(8.16) 
that contribute at 0(e 4 ) to the charge radii and are thus neglected. Insertions of ^4+ 
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into the magnetic moment terms produce 

C = ib[(-)^ + ^ ){llk+v ^B kji [S^S u }{Q + ,A + } kk 'B ijk 'Ff lu 

+ib' 2 B kji [S^ S v \{Q + ,A + }**B i ' ik Fi u ' + ib'^-^'^+^B^iS^, S^Qf^&'^F^ 

+ib' 4 {-) r i^ + ^B kji [S fl , S„] ^A%B^' k F^ 

+ib' 5 (-)^' +1 B kji [S^S u } (Q% ' A j l + A% 'fif ) B i,j ' k F^ 

+i [b' 6 (B[S li ,S„]BQ+) + ti 7 {B[S^S U }Q + B)} F^stx{A + ) 

+ib' s (BIS^, S„]B) F^str(Q + A+), (8.17) 

which are 0(e 2 ) corrections to the magnetic moments and can be discarded [161]. 

Finally we assess the contribution from the operator Of in Eq. A8.5I I. As in the 
meson sector, the charge preserving operators can be constructed by the replacement 
Q — > oA x ci 7 q in leading-order terms. Again it is easier to embed these operators in £ so 
that gauge invariance is transparent. For the charge radius, the leading vector-current 
correction operator is contained in the term 

C = -?- [c A ,a (BBc 1>Q ) + c A)0 (Bc liQ B)] v^d u F^, (8.18) 

which leads to C(e 4 ) corrections. For the magnetic moment operator, such a replace- 
ment leads to 

f ft 

C = - [n A , a {B[S^ S u ]Bc hQ ) + /m,/3 (B[S», S u ]c 1)Q B)] F^, (8.19) 

and corrections that are of higher order than the one-loop results [161]. See Ap- 
pendix 1^2] for results in an alternate power counting scheme. 
To C(e 2 ) the baryon charge radii are thus 



<r E > - 1 ^{Qc-+a D c + ) + ^ T {QiiF + OLD^D) 



16tt 2 / 2 ^ 

J x 



2 2 

A x log ^ - 5 (3 X log ^f + 10 (3' x G{m x ,A,n) 



(8.20) 



and the magnetic moments to 0(e) read 

M L 
4^/ 



/i = (Q n F + a D /jl d ) + — ^^ [Pxmx +/3 x J r (mx,A,IJ,j\ ■ (8.21) 



x 

The a-dependence is treated as implicit in the meson masses. The PQxPT coefficients 
A x , 0x, and f3' x can be found in Tables 15^2145. 101 of Chapter |5] and in Ref. [31] along 
with the functions f(mx,A )( t() and 5(mi,A,/i). The quenched charge radii at 0(e 2 ) 
are similar in form (although A x = in Q%PT) due to the lack of singlet contribu- 
tions at this order. The QxPT coefficients (3 X and j3 x appear in [4, 13]. The quenched 
magnetic moments, however, receive additional contributions from singlet loops. The 
relevant formula of [13] are not duplicated here in the interests of space but only need 
trivial modification by taking into account the a-dependence of meson masses. 
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8.5 Decuplet Baryon Properties 

Decuplet matrix elements of the electromagnetic current J p have been calculated in 
Chapter [H] where it was found that the charge radii are C(e 2 ) at NLO in the chiral 
expansion, while the magnetic moments are 0(e) and the electric quadrupole moments 
are O(e ). At one-loop order in the chiral expansion, the magnetic octupole moment is 
zero. 

There are no finite-a operators in Eq. ( 18.91 1 that contribute to decuplet baryon form 
factors. The SW term can potentially contribute when A+ is inserted into the La- 
grangian. There are three such terms: the charge radius, magnetic moment, and 
electric quadrupole terms. Insertions of A+ into the charge radius term produces 

A x A x 

+^- {T°Q + %) v^F^ str(A+) + ^ {T°%) v^F^ str(Q + A+). (8.22) 

A x A x 

These contribute to decuplet charge radii at C(e 4 ). As in the octet sector, insertions of 
A+ into the magnetic moment term, namely 

C = id' 1 T^ i {Q + ,A + } ii 'T^' jk F pu + id' 2 (-)^' { ^ +r >i' ) T'v i q i ^A J i'TJ'i' k F^ 

+i 4 (T^Q+%) F^ str(A+) + * < (T^T U ) F^ Btr(Q+A+), (8.23) 

produce 0(e 2 ) corrections. Likewise, insertions of A + into the electric quadrupole term 
have the form 1 

c = ^-r^^ i {Q + ,A + } ii 'r^^ k v a d lx F ua 

A x 
d" 



A X 

f- fr^Q+T^) v a d fi F uaS tr(A+) + f 

v x v J A x 



X. 

-3- (r^Q+T^ v a d fi F uaS tr(A+) + -^ (T^T^) v a d^F va str(Q + A +! 



(8.24) 



and produce 0(e 2 ) corrections. All of these corrections are of higher order than the 
one-loop results. 

Finally we assess the contribution from the operator Of in Eq. H8.5I ). The effective 
operators can be constructed by replacing Q by aA x c^Q in LO terms. Embedding these 
terms in a Lagrangian, we have 

C = ^^{T°ci >Q %)v„d v F^ + MaiiA, c ^»c m %)F» v 

X. 



:H; "^c 1]g f})^ a . (8.25) 



A_ 



lr rhe action of * '"' on Lorentz indices produces the symmetric traceless part of the tensor, viz., O^"^ — 
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Each of these terms leads to corrections of higher order than the one-loop results and 
can be dropped. Thus at this order the only finite lattice spacing corrections to decuplet 
electromagnetic properties appear in the meson masses. For reference, the expressions 
are 



-*> = n^r + ^r 

19 + 5C 2 v- , m\ 25 U 2 v-^ , „, A , ,„„,., 

-3li^E^>°^-^I^ ? E^5(m.v,A, fl ), (8.26) 

M B W 2 C 2 M B sr^ 

\i = 2fi c Q - 3 2 2^ Aif(mi,A, /i) - — ^ 2^ ^ m ^. (8.27) 

X X 



and 



/ ^ 2M 2 \ M 2 C 2 v^ . , ml MlH 2 *-^ , „, an 



(8.28) 

The coefficients ylx are tabulated in Table I67T1 Extending the result to QxPT, where 
Ax = 0, one must include additional contributions from singlet loops. With finite 
lattice spacing corrections, the expressions are identical to those in Chapter |6] except 
with masses given by Eq. ( 18.8ft . Thus for brevity we do not reproduce them here. 

8.6 Decuplet to Octet Baryon Transition Properties 

The decuplet to octet matrix elements of the electromagnetic current J M have been 
calculated in Chapter [7J There we found, that at next-to-leading order in the chiral 
expansion, G\(q 2 ) is 0(e) while G 2 {q 2 ) and G 3 (q 2 ) are 0(e ). 2 

There are no new finite-a operators in Eq. ( 18.9ft that contribute to decuplet to octet 
transition form factors. Insertion of A+ into leading-order transition terms leads to 
corrections of 0(e 2 ) or smaller. For completeness the terms are: 

£ = it{B kji S^ &lA%X' jk F» v + it 2 B kji S^A% Q%T^ k F^ 

+it 5 (BS^Q+%) F^ str(A+), (8.29) 

2 Here, we count e ~ A/ Mb. 
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for the magnetic dipole transition; and 

£ = ±LB kji S^Q%A l £T^' i ' jk v a d IJl F va + ^-B kji S^A%G^T v ^ jk v a d il F va 

A x A x 

A x 

+^-{-) vAvi+n i' ) B kji S {lx A%^T^^ j ' k v a d^F va 

A x 



+|l (BS^Q + T^ v a d ll F ua str{A+), (8.30) 



for the quadrupole transition. Finally, insertion of A+ into the PQxPT term propor- 
tional to i(BS^QT u )d a d fl F ua leads to 

C = ^B^S^ Q%A l (V' jk d a d»F v a + ^-B^S^Af Q%T^ k d a d»F v a 



X X 

+ l -k(-^A'nj+V j ')Qkji S ^ A Jj' T i'fk Q«Q» F V a 

X 

+ Tj ( " f 1 ' {r ' 3 +r ' 3 ' ] B kji S^A% Q^'r^' k d a d^F v a 
A x 

+|f (BS^Q+%) d a d»F\ str(A+), (8.31) 

x 

for the Coulomb quadrupole transition. 

Similarly, constructing Of in the effective theory by replacing Q with aA x ci t Q in the 
transition operators leads to terms of at least 0(e 2 ) which are contained in the terms 

C = ia^ T ^(BS^ Q T u )F^ + ^^^(BS^c 1)Q T^v a d^F va 

X, 

+! ^r I Vl ( Bs ^Q TU ) ^ V- (8.32) 

A, 

All of these corrections from effective Of operators are of higher order than the one- 
loop results. Thus at this order, the only finite lattice spacing corrections to the tran- 
sition moments appear in the meson masses. For reference the expressions are 



Ci(0) = ^a T -^HC^Y.Pxj dx(l-|)jP(m x ,xA,/i) 

A X 

+4ttC(D - F)^- V/?| [ dx{l-x)F{m x ,-xk,ii), (8.33) 



x x 
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and 



M 2 I m f 1 T (\ _ T \ 

G 2 (0) = jf-l -4.® T a T + WHCj2PxJ dx V g > g(m x ,xA,») 

-16C(D - F)T] fix / da;i(l-i)6(mx I -iA, J u) I, (8.34) 

x Jo J 

g 3 (o) = -16-rfX// dxx(i-x)(x 



x — .,0 2; m^-xW 

\ m ^n[— ) +C p-F)/3|7ef-— )1 . (8.35) 

3 \m x J V ™x/_ 

The coefficients (3t and /9s are given in Tables l7l2l and rr3l the function lZ(x) is defined 
in Eq. JB.2D . Extending the result to QxPT, where the coefficients are replaced with 
their quenched counterparts p§ and (3®, one must include additional contributions 
from singlet loops. With finite lattice spacing corrections the expressions are identical 
to those in Chap ter[7] except with masses given by Eq. J8.8D . Thus for brevity we do not 
reproduce them here. 

8.7 Conclusions 

In this chapter we have calculated the finite lattice spacing corrections to hadronic 
electromagnetic observables in both QxPT and PQxPT for the SU(3) flavor group in 
the isospin limit and the SU(2) group with non-degenerate quarks. In the power count- 
ing scheme of [159, 161], 0(a) corrections contribute to electromagnetic observables at 
higher order than the one-loop chiral corrections. Thus finite lattice spacing manifests 
itself only in the meson masses at this order. 

In practice one should not adhere rigidly to a particular power-counting scheme. 
Each observable should be treated on a case by case basis. The actual size of a and 
additionally the size of counterterms are needed to address the relevance of 0(a) cor- 
rections for real lattice data. For this reason we have presented an exhaustive list 
of 0(a) operators relevant for hadronic electromagnetic properties. In an alternate 
power counting for a coarser lattice (as explained in Appendix IG.2I ). some of the op- 
erators listed above contribute at the same order as the one-loop results in the chiral 
expansion. 

The corrections detailed in Appendix lG.2l in the baryon sector may also be necessary 
if one goes beyond the heavy baryon limit (that is, including I /Mb corrections). For 
example, in the case of the octet baryon magnetic moment [see Eq. H8.21H at NLO in 
the heavy baryon expansion \x would be known to 0(e 2 ). Thus 0(a) corrections in the 
power counting of Eq. ( 18.21 ) are needed since they are also 0(e 2 ). 

Knowledge of the low-energy behavior of PQQCD at finite lattice spacing is crucial 
to extrapolate lattice calculations from the quark masses used on a finite lattice to the 
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physical world. The formal behavior of the PQQCD electromagnetic observables in the 
chiral limit has the same form as in QCD. Moreover, there is a well-defined connection 
to QCD and one can reliably extrapolate lattice results down to the quark masses 
of reality. For simulations using unimproved lattice actions (with Wilson quarks or 
mixed quarks), our results will aid in the continuum extrapolation and will help lattice 
simulations make contact with real-world data. 



97 



Chapter 9 

Summary 



In this thesis we have presented model-independent analytic results for various ob- 
servables in the one meson and one baryon sectors in the effective field theories QxPT 
and PQxPT. These results are needed to extrapolate lattice QCD simulations which 
use light quarks that are heavier than those in nature. 

In particular, LQCD simulations that use the quenched approximations need to use 
the appropriate low-energy theory, QxPT, to do this extrapolations. Although QxPT 
is the proper theory to extrapolate quenched simulations, it is problematic to draw 
conclusions from the results for real-world QCD. The reason for this is that the flavor 
singlet — the equivalent of the rf in QCD — is not heavy in QQCD and cannot be inte- 
grated out. As a consequence, QxPT results are usually plagued by quenching artifacts 
and found to be more divergent in the infrared limit (m q — > 0) than their x?T coun- 
terparts. More general there is no relation between the low-energy constants in QxPT 
and those in xPT. Hence, extrapolated quenched lattice data is unrelated to QCD. 

Recently, more and more lattice QCD simulations are performed using the partially 
quenched approximation of QCD. The proper method of extrapolating PQQCD data 
to the physical regime is to use PQxPT. In contrast to QxPT, PQxPT does have an 
analytic connection the QCD: In the limit where the sea quark masses are equal to 
the valence (and ghost) quark masses one recovers QCD. In particular, the low-energy 
constants of PQxPT have the same numerical values as their counterparts in xPT. 
Thus, PQxPT not only enables a clean extrapolation of PQQCD lattice data to the 
physical regime but it also accurate physical predictions for the real world: QCD. 

The first half of this thesis (Chapters |U and HJ is concerned with the calculation of 
several observables in the heavy-light meson sector: 

In Chapter|3l we calculate chiral 1/M 2 corrections to the semileptonic B(*> — » Z?W 
decays at zero recoil in QxPT that are due to the breaking of heavy quark symmetry. In 
Chapter HJ the emphasize is shifted to the investigation of finite volume effects in the 
lattice QCD treatment of the heavy-light meson sector. In particular, we investigate 
the role of the vector-pseudoscalar mass splitting, A. We find that finite volume effects 
arising from the propagation of Goldstone mesons in the effective theory are significant 



98 



and can be altered by the presence of the scale A. 

The second part of this work (Chapters [5J-[71l contains a number of calculations of 
hadronic properties in the one-baryon sector: 

Specifically, in Chapter [5] we calculate the electric charge radii of the SU(3) octet 
baryons in QxPT and PQxPT (we also include this calculation for the SU(3) pseu- 
doscalar mesons). We find that in the QxPT calculation new operators, which appear 
because the flavor singlet must be retained, enter at NNLO. Although these do not 
render the quenched NLO result more divergent than its QCD counterpart, quenching 
artifacts do show up: Not only are the low-energy constants different in QxPT and 
xPT, but for certain baryons the diagrams which have bosonic or fermionic mesons 
running in loops cancel so that the quenched result is actually independent of ttiq ! We 
come to similar findings in Chapter|6l where we calculate electromagnetic properties of 
the decuplet baryons in QxPT and PQxPT. Here, the expansions about the chiral limit 
for QCD and QQCD charge radii are formally similar, but the QQCD result consists 
entirely of quenched oddities. In Chapter |71 we determine baryon decuplet to octet 
electromagnetic transition form factors in QxPT and PQxPT and — once again — come 
to similar results: In contrast to the quenched transition moments that pick up con- 
tributions from hairpin loops, the PQxPT is analytically connected to the x?T results: 
The low-energy parameters have similar values in the two theories. In Chapter |8] we 
augment all these calculations to include 0{a) corrections which are due to lattice dis- 
cretization using two different power-counting schemes. Our results are important to 
extrapolate simulations that use unimproved lattice actions (with Wilson quarks or 
mixed quarks). 



This is an exciting time for nuclear physics as for the first time rigorous predictions 
for the structure and interactions of nuclei from QCD using lattice simulations seem 
within reach. This development is caused by the availability of faster computers as 
well as by the conceptual advances in lattice computing algorithms and it will enable 
the simulation of many baryon properties with improved precision in PQQCD. How- 
ever, it will be some time before real simulations with light physical quarks become 
feasible. Until then, the lattice results for all these quantities need to be extrapolated 
to real-world QCD using PQxPT. Furthermore, effects due to finite lattice spacing and 
due to finite lattice volume must be taken into account and included in the PQxPT 
treatment. Many quantities in the one-hadron sector still await such treatment, mak- 
ing this a very exciting area for future research. 

Another thrilling avenue of research is the two-nucleon sector. Although this sys- 
tem is much more complicated than a single nucleon, lattice simulations of this system 
now appear feasible and promise predictions about the nucleon-nucleon interaction 
directly from QCD. The lattice QCD treatment of the two-nucleon sector has just be- 
gun [84-86, 114], is still largely uninvestigated, and offers great opportunity for future 
research. 
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Appendix A 



Baryon Transformations for 
Flavor SU(2\2) and SU(4\2 



Table A.l: Embedding of the baryon doublet and quartet for SU(2\2) V for QQCD. 



qqq 
qqq 
qqq 
qqq 



Doublet 
SU (2) y SU(2) e dim 




20 



Quartet 

SU(2) V SU(2) g dim 



(4,1) 
(3,2) 
(2,3) 



4 
6 
6 




16 



T able A.2: Embedding of the baryon doublet and quartet for SU(4\2) V for PQQC D. 
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SU(2) Y SU(2) B SU(2) g 


dim 


SJ7(2)y 
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4 


qqq* 
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qq s q s 
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8 




(2,3,1) 
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q s q s q s 
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qq s q 
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70 


44 
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Appendix B 



Formulae relevant for B { * 
at Zero Recoil in QxPT 



D 



* 



We list the functions H\, H 2 , F\, H 5 , H 8 , and G5 (some of which have appeared in the 
literature before [45,60]). Here, m = m qq is the mass of the qq light meson in the loop 
where q = u, d, or s is the light (spectator) quark content of the heavy mesons. We have 
calculated loop integrals in d = 4 — 2e dimensions and used dimensional regularization 
with the minimal subtraction (MS) scheme, where 



1/e = l/e-7£ + log4vr + l. 
As a shorthand we have defined the function 



R(x) = v x 2 — 1 log 



\[W 



1 + ie 



+ V 7 ? 



1 + iej 

which occurs frequently. We also need its derivative dR/dx given by 

x 



(B.l) 



(B.2) 



R'(x) 



ar 



1 



R(x) 



(B.3) 



For the calculation of the wave function renormalization contribution we need the 
derivatives of the loop integrals for the diagrams in Fig. 13.11 



ffi(A) 



167T 2 



log- 



m 



1 



1 



M' 



m 



(B.4) 



H 2 (A) 



16vr 2 



16 



10 



_A Z - — m z 



>-. < m 2 1\ 4 A /A 

2(m 2 - A 2 ) log -^ - - + -Ami? — 



+ ( £ A 2 - -m 2 
1 3 3 



TO 



(B.5) 
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and 



Fx(A) 



16vr 2 



10 



.2 4, 



m 



3 AZ " 3^ + ^ " 2A ^ l l0g 7? " ?) + ^ (£ 



+ -(A 2 -m 2 )i?' 



,/AN 



3 \my 

For the loop integrals of the vertex corrections one finds 



(B.6) 



H 5 (A,A) 



16vr 



m 1 

2 1 lo § -2" ~ "7 ~ X 



/» 



A-A 



ml \ m 



(B.7) 



H 8 (A,A) 



16vr 2 



2m 2 



:(a 2 + aA + A 2 



log 



m 1 \ 16 , . 9 . r io 

^-?J + ¥ (A+AA + A 



10 9 m5m 2 - 2A 2 ) / A 

— m 2 H = -R — 

3 3(A-A) U 



m(5m 2 -2A 2 ) /A\ 



3(A-A) 



R - 



m) I 

(B.8) 



and 



G 5 (A,A) 



16vr 2 I 9 
+ 



-(A 2 + AA + A 2 )--m 2 



4 

— J 

3 



m 2 --(A 2 + AA + A 2 ^ 



m 1 
log — - - 
H z e 



A 



2m(A 2 - m 2 ) ( A\ 2m(A 2 - m 2 ) _. 

H = R — = R — 

3(A-A) \mj 3(A-A) \m 



(B.9) 
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Appendix C 

Formulae relevant for HMxPT in 
a Finite Volume 



C.l Integrals and Sums 

We have regularized ultra-violet divergences that appear in loop integrals using di- 
mensional regularization with the MS scheme [see Eq. dB.lH . The integrals appearing 
in the full QCD calculation are defined by 



I- X (m) = ix 



4-d 



d d k 



1 



ira 



H- x (m,A) = fcr-vV)^ 



(2vr) d k 2 - m 2 + it 16vr 2 
d f d d k 



1 . / ' m 1 
- - log — 



kinKy 



(C.l) 



dAJ (2n) d (k 2 - m 2 + ie)(v ■ k - A + ie) 



where 



F- x (m,A) 



3^KA), 



1 / m 

- - log — 



(C.2) 



16tH 



/2A 2 



V 3 

2(A 2 -m 2 ) /A 

+ J -mR — 

6 \m 



m 2 ) A + 



10A 2 Am z 



9 



with 



R(x) = \/x 2 — 1 log 



\/x 2 — 1 + ie 



+ Vx 2 - 1 + 



(C.3) 



(C.4) 



ie 



and /i is the renormalisation scale. For the quenched and partially quenched calcula- 
tions, we also need the integrals 



Av') = ^-d 



d d k 



1 



(27r) d (k 2 -m 2 + ie) 2 
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dm 2 



(C.5) 



and 

(J j (X rv Ki pKip 



H«(m,A) = (gP»- v P v ")^- d JLJ. 



dAJ (2-ir) d (k 2 -m 2 + ie) 2 (v-k- A + ie) 
"' H- x (m,A). (C.6) 



dm 2 



In a cubic spatial box of extent L in four dimension with periodic boundary conditions, 
one obtains the sums (after subtracting the ultra-violet divergences) 

r "»4£/¥ f-rf + <. =; " , » + J ""»»' <a7) 

k 
and 

(9 /" dk$ k p k u 



H(m,A) = (^-^)^E^/ 



27r (fc 2 — m 2 + ie)(v ■ k — A + ie) 
= H(m,A) + H Fy (m,A) (C.8) 

for the full QCD calculation, where the momentum k is quantized according to k = 

(2ir/L)n. Furthermore, I(m) = I^(m)\^ =0 and H(m) = H^(m,A)\ x=0 are the infinite 
volume limits of 2 and H, and (n = \n\) 

im s—^ K\inmV) 



j_ jn7^ e _ nmL _J 15 

47r 2 ^V2n 3 L 3 I 8nmL 128(nmL) 2 



nmL 



(C.9) 



is the finite volume correction to I(m) in the limit mi » 1. The function #pv is the 
finite volume correction to H{m, A) and can be obtained via 

#Fv(m, A) = i [(m 2 - A 2 )K FV (m, A) - 2A J FV (m, A) + iIfvM] , (CIO) 

where JFv{m, A) and Kpyfa, A) are defined in Eqs. ( 14.201 1 and < I4.27I >. 
For QQCD and PQQCD calculations, one also needs 

tv' I \ \~~^ /" rt \ t i \ tn in 

L 3 ^^ 7 27r (fc 2 — m 2 + ie) 2 <9m 2 9m 2 

and 

d f dko kpk v 



W\m,A) ee (f-//)^M/' 



27r (A; 2 - m 2 + ie) 2 (v • A; - A + ie) 

k 

' H(m,A) + 1 —H FV (m,A). (C.12) 



dm 2 dm 2 
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C.2 One-Loop Results 



We collect the results for one-loop corrections to fp, s) */Mp, s) and Bp, s) . For conve- 
nience, we introduce 



and 



C±(m, x) = T(m) ± g 7i(m, x) 



Cl (m,x)= J"' (m) ± g 2 ^' (m, x) , 



(C.13) 



(C.14) 



where the functions T(m), H(m, x), I 71 (m) and H v (m, x) are defined in Eas. dC.7l l. dC.8l >. 
dC.lll l. and dC.121 ). respectively. 
In full QCD, we find 



f P ^fMp~ = J 1 - -^ [9C_(M^, A) + 6C_(Af x , A + 8 a ) + C.(M V , A)] I , (C.15) 



12/ 



f Ps y/Mp~ s = k <| 1 - ^ [3C_(M^, A - 5 S ) +C_(M„, A)] \ , 



(C.16) 



5p = H J 1 " 6^ [3C+(M^, A) +C + (M„ A)] 
^ = ^{l-^[ C+ (M„A)]' 



(C.17) 
(C.18) 



In QQCD, we find 



fp ^/M p - = K \i + — 



|C_(M^, A) + aM ' 3 M °V (M^,A) + 2 57 W(M ff , A) 



(C.19) 



/ftV^^l+or. 



272 



a - . . v . amf s - M ^. 



3 -C_(m ss ,A) + 



CI(m M ,A) + 2g7W(m Sfl ,A) 



(C.20) 



Bp = i{ 1 -^ 



H 3/? /l * 



where 



1 - |) C + (M n , A) - " M ' 3 M ° C!l(M n , A) + 2^^, A) 



(C.21) 



1 - -) C + (m ss , A) 2£_ °_d (m ss , A) + 2 57 W(m ss , A) 



m. 



2M2 - M2. 



(C.22) 
(C.23) 
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In PQQCD, we find 



i 

2/2 



2C_ (m uj , A + <5 sca ) + C_ (m ur , A + 5 SC3i + S s 



+ TT 



,2 

rr 



,2 



27 \ m x - m uu 



3 m^ u — m x ' • > <"~ 

jj ''^uui \--~rr ■■~uui pri 



l(m 2 -m 2 J(m 2 r -m 2 J 



i 2 

"uu 



A 



Cl(m uu ,A) 



" l C_(m x ,A) 
(C.24) 



ipsVMp: 



k{ 1 



2P 



2C_ (m S j, A + <5 sca - 8 S ) + C_ (m sr , A + <5 sca + <5 S - <5 S ) 



'■j:i 



3 m 2 



m 



ss -C_(m ss , A) + — 



/?)" 



m 



x 



27 \ m 2 ^ 



m 



1 (m^- - m 2 ss )(m 2 r - m 2 ss ) ^ 



m- 



m 



Cl(m ss ,A) 



x 



" ] C_(m x ,A) 
(C.25) 
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B Ps 



3/3 
2^2 



i 

T 2 



C + (m uu ,A) 



1 (mjj - m 2 J(m 2 r - mjj , 

Li [jTluui ^) 



>nt 



m 



x 



1 m 2 r - 
3< 



m' 



m 



*C+(m uu ,A) + — 



m; 



m 



.u 



a 



27 \ m 2 



m 



C+(mx,A) 



x 



(C.26) 



3/3 
2^? 



C + (m„, A) - l^L_^X^_^ C %(rn ss , A) 



1 (m 2 - - m 2 ss )(m 2 r - m 2 s ) , 



? 'S.S' 



4 



9 9 

1 m^- — m; Q 

+ o^ ^C + (m ss ,A) + - 

3 m| s — m^ 



2 I m 2 „. - m 2 - 



A 



C+(m x ,A) 



where 



m X = o ( m jj + 2m ?r) 



(C.27) 



(C.28) 



It is straightforward to show that the PQQCD results reproduce those for full QCD in 
the limit rrij = m u and m r = m s . 
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Appendix D 

Charge Radii of the Meson and 
Baryon Octets for Flavor SU{2 



We consider the case of SU(2) flavor and calculate charge radii for the pions and nu- 
cleons. We keep the up and down quark masses non-degenerate and similarly for the 
sea-quarks. Thus the quark mass matrix reads m Q *• ' = diag(m u , md, mj,mi,m u , ma). 
Defining ghost and sea quark charges is constrained only by the restriction that QCD 
be recovered in the limit of appropriately degenerate quark masses. Thus the most 
general form of the charge matrix is 

Q su ^ =diag(|, -£,<&, «,<&,«) . (D.l) 

The symmetry breaking pattern is assumed to be SU(A\2)l SU(4\2)r <& U(l)v — ► 
SU{4\2) v ®U{l) v . 

For the ir + , ir~ , and ir° we find 



<#V) = 1+ ' 



16vr 2 / 2 



o + 1i ) F dd + f o ~ Qj ) F uu - (1 - Qj + qi) Kd 
3 + Qj j F jd - ( - + qi J Fia - ( - - qj J F ju - ( - - q, 1 h)„ 



+ a 9 ^q 2 , (D.2) 

g pq = -c^, and G^ = 0, respectively. 

The baryon field assignments are analogous to the case of SU(3) flavor. The nucle- 
ons are embedded as 

B ijk = -= {eijN k + e ik Nj) , (D.3) 

where the indices i,j and k are restricted to 1 or 2 and the SU(2) nucleon doublet is 
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defined as 



■n 



N=[ P ) (D.4) 



The decuplet field %jk, which is totally symmetric, is normalized to contain the A 
resonances T^ = %jk with i, j, k restricted to 1 or 2. The spin-3/2 baryon quartet is 
then contained as 

T m = A++, T 112 = i=A + , T 122 = i=A , andT 222 = A-. (D.5) 

The construction of the octet and decuplet baryons containing one sea or one ghost 
quark is analogous to the SU(3) flavor case [94] and we will not repeat it here. 

The free Lagrangian for B and T is the one in Eq. ( 15.171 1 (with the parameters 
having different numerical values than in the SU(3) case). The connection to QCD is 
detailed in [94]. Similarly, the Lagrangian describing the interaction of the B and T 
with the pseudo-Goldstone bosons is the one in Eq. A5.17I I. Matching it to the familiar 
one in QCD (by restricting the B^ and %jk to the qqq sector) 

C = 2g A NS»A fl N + ^iV^iVtr^) + g AN (rf A^N^m + h.c.) (D.6) 

one finds at tree-level 

4 1 2 1 

a = -gA + ^gi, P = -gi--gAi and C = -gAN- (D.7) 

The contribution at leading order to the charge radii from the Pauli form fac- 
tor F 2 (g 2 ), involves only the magnetic moments which arise from the PQQCD La- 
grangian [94] 



C 



2M N 



fi a (b[s^s u \bq su{ - 2) ) +H0 (b[s^s v \q su{ < 2) b) 



+ ^str(Q su ^){B[S fl ,S u ]B) 



F^. (D.8) 



Note that in the case of SU{2) flavor the charge matrix Q is not supertraceless and 
hence there appears a third operator. In QCD, the corresponding Lagrange density is 
conventionally written in terms of isoscalar and isovector couplings 

C = 2^- {jtoN[S„, S U ]N + niNiS", S u ]t 3 N) F^ (D.9) 

and one finds that the QCD and PQQCD coefficients are related by [94] 

Vo = t, (/•*<* + M/3 + 2 /u 7 ) , and m = - (2// a - //#) . (D.10) 
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Table D.l: The coefficients j3 x , P' x , and A 


x in 


5(7(2) flavor PQxPT for the proton. 
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(3x 
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A x 


uu 


l(4gi + 2g 1 g A + gt)(2-3q j ) 






( — 27 + TgQjWAN 
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ud 
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-lgf(l + 3 qi ) 
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'27 — Tr1j)9an 
(27 ~ TE1 1 )9an 
\ 27 — q9j)9an 

( _ 2? ~ f)1l)9AN 
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uj 


-\{Ag A + 2g 1 g A +gl){2-Zq ] ) 
-l(4g A + 2g igA +gf)(2-3 qi ) 






f"2ft 


ul 






|-2ft 


dj 








-1-ft 


dl 







Table D.2: The coefficients /3 X , /? x , and A x in 5(7(2) flavor PQxPT for the neutron. 
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A x 



UU 

ud 
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dl 
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- 2giffA 
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|(4gj + 2g igA +^)(l 



3 '* 
li) (g + Tsft + \<U)9an - l + 2 ft + « 

1 + 2* 
1-3j 



- 3 ft) 
-3ft) 



(s + T$1i)9an 
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AJV 



(27 ~~ q1u9an 
(~ 54 ~ ~ml3)9 AN 
(~ 54 Z SJWAJV 



5 -ft 

-!-2ft 
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Likewise, the leading tree-level corrections to the charge-radii come from the La- 
grangian 



£ = 

A 2 



c a (£!£<2^ 2 )) + cp (BQ SUi2) B) + c 7 str(Q^ 2 ))(l!/3) 



v li d v F> u ' 



(D.ll) 



that matches onto the QCD Lagrangian 



£ = -j (co iViV + ci /Vr 3 ^) t^a^"" 
x 



(D.12) 



with 



1 1 

co = - (c a + C/3 + 2c 7 ) , and c x = - (2c a - c-g) . 



Evaluating the charge radii at NLO order in the chiral expansion yields 



(D.13) 



9 6c 3a 

< t% >= -— + 



A 2 2M 2 N 16ir 2 f 



r 2f2 2^ 



A 



2 2 

777 777 

A x log -f - 5 /3 X log -f + WP x g(m x , A, M ) 



(D.14) 

The coefficients c are given by c p = cq + c\ and c n = cq — c\ while a p = /j,q + /j,% and 
a n = fiQ — (Mi. The remaining coefficients are listed in Table ILvfl for the proton and 
Table El for the neutron. 
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Appendix E 

More on the Baryon Decuplet 
Form Factors 

E.l q 2 Dependence of the Form Factors 

For reference, we provide the q 2 dependence of the decuplet electromagnetic form fac- 
tors defined in Section 1(01 at one-loop order in the chiral expansion. To do so we define 



Px 



x(l — x)q 2 



mA 



Then we have 

Fi(q 2 ) 



Q[i- 

3 + C 2 

16tt 2 / 5 

n 2 

24tT 2 / 2 
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Qcq 2 , c c q 2 



2M| 



2A| 
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A I 



S> 



X 



X 



Triog- 



m 
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p 



2m\ 
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Jo 



dxP x log P x 



sr^ A ( 11 2 , ™"Y 5 . ^/A\ 

Y^A X { — q 2 \og^ + -Am x K 

x \ 36 n 2 3 \ m W 



d.x 



10 /m 



'x 



A 2 -— x(l-x)g 2 J lot-, P s 



r , ,5 x(l — x)g 2 



\m x P x ) 



(E.l) 



(E.2) 
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F 2 ( q 2 ) = 2 ^Q-^-'£A x m x j\p x + ^^J2 A 4^ ] ^ 



and 



m 



x 



p 



J dx 2AlogP x -m x P x n(^ r ) I, 



+ / dx 



(E.3) 



4Q[fi c + ' 



+ 



2M|\ M 2 B C 2 



Y^A X ^log^f+ f dx2x{l-x)\ogP x 
x L 6 V Jo 



2Mpi 2 
9vr 2 / 2 



X> {>■ 



X 



o + / waJ»(l — 27 

6 M 2 Jo 



21ogP x 



Am x Px v ( A \ 
A2-m 2 Y P2 l mx pj 



(E.4) 



E.2 Electromagnetic Properties for Flavor SU(2) 

Here we consider the case of SU{2) flavor and calculate the electromagnetic moments 
and charge radii of the delta quartet. We keep the up and down valence quark masses 
non-degenerate and similarly for the sea-quarks. Thus the quark mass matrix reads 
rag ( ' = diag(m u ,m d ,mj,mi,m u ,m d ). Defining ghost and sea quark charges is con- 
strained only by the restriction that QCD be recovered in the limit of appropriately 
degenerate quark masses. Thus the most general form of the charge matrix is 



Q SU{2) =diag( -,--, qj ,qi, qj ,qi 



(E.5) 



The symmetry breaking pattern is assumed to be SU(A\2)l <8> SU(4\2)r <g> U(l)v — ► 
SU(A\2)v <8> U(l)v- The baryon field assignments are analogous to the case of 51/(3) 
flavor. The nucleons are embedded as 



1 



V6 



(E.6) 



where the indices i, j and k are restricted to 1 or 2 and the 577(2) nucleon doublet is 
defined as 

V 



N 



n 



(E.7) 



The decuplet field %jk, which is totally symmetric, is normalized to contain the A- 
resonances T^ = %jk with i, j, k restricted to 1 or 2. Our states are normalized so 
that Tin = A ++ . The construction of the octet and decuplet baryons containing one 
sea or one ghost quark is analogous to the 577(3) flavor case [94] and we will not repeat 
it here. 
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Table E.l: The SU(2) coefficients A x in X PT and PQxPT. 
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3« 
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1 1 



9i 
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' 7 ' 



The free Lagrangian for B and T is the one in Eq. ( 12 .601 1 (with the parameters 
having different numerical values than the 577(3) case). The connection to QCD is 
detailed in [94]. Similarly, the Lagrangian describing the interaction of the B and T 
with the pseudo-Goldstone bosons is the one in Eq. ( I5.17I I. Matching it to the familiar 
one in QCD (by restricting the B^p. and %j k to the qqq sector), 



C 



9AN 



Tt ji A v u N 3 e kl + h.c) + 2g AA Tl ji S^ l T^ ljk + 2g x T% ji S„T 1/jijk tr(A>*), 



(E.8) 



one finds at tree-level C = —g&N and H = #aa, with gx = 0. The leading tree-level op- 
erators which contribute to A electromagnetic properties are the same as in Eqs. ( 16.71 ). 
( 16.91 ). and 16.1H . of course the low-energy constants have different values. 

Evaluating the A electromagnetic properties at NLO in the chiral expansion yields 
expressions identical in form to those above Eqs. H6.15H . d6.16D . and A6.17I ) with the 
577(2) identifications made for C and H above. The SU(2) coefficients A x appear in 
Table lE~H for particular A-resonance states T. [tb] In the table, we have listed values 
corresponding to the loop meson that has mass mx for both x?T and PQxPT. Again, 
the xPT coefficients can be used to find the A-resonance charge radii in two-flavor 
QCD. These have not been previously calculated. 

In addition, however, local counterterms appear, that involve the non-zero super- 
trace of the charge matrix in S77(2|2). Using the general form of the charge matrix 
Eq. dE.5l >. we have an additional dimension-5 magnetic moment operator in PQxPT 



3ie 



that matches onto the xPT operator 

There is an additional dimension-6 electric quadrupole operator in PQxPT 

3e 



£ 



27 A 2 
x 



T^T u n v a d fl F ua str(Q su ^), 



that matches onto the xPT operator 



^ 7f<{li.rpV} 



C -Q^Tflf ^F ra tr(Q s ^ 2 )). 

x 



(E.9) 



(E.10) 



(E.ll) 



(E.12) 
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Finally in PQxPT there is an additional dimension-6 charge radius operator 

C = c 7 |j (T a %) VpdvFV str (Q 5C/ ( 2 )), (E.13) 

A x 

that matches onto 

C = c^TIT^v^F^ tr(Q su ^) (E.14) 

X. 

in xPT. Notice the PQxPT low-energy constants // 7 , Q 7 , and c 7 are identical at tree 
level to those in xPT. 

Inclusion of the above operators leads to tree-level contributions to the A quartet 
electromagnetic properties. Since these contributions are proportional to the super- 
trace of the charge matrix, the corrections are identical for each member of the quartet. 
The charge radius should include an additive correction 

while for the magnetic moment 

5fx = 2/x 7 , (E.16) 

and for the electric quadrupole moment 

M 2 
<5Q = -2 /[ / 7 -4Q 7 ^f. (E.17) 

X. 

Notice that these corrections only affect the counterterm structure of the results. 
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Appendix F 

A — > A^7 Transitions for Flavor 

sua 



We repeat the calculation of the transition moments for the case of SU(2) fla- 
vor with non-degenerate quarks, i.e., the quark mass matrix reads m Q = 
diag(m u , md,mj,mi,m u ,md)- Since defining ghost and sea quark charges is con- 
strained only by the restriction that QCD be recovered in the limit of appropriately 
degenerate quark masses, the most general form of the charge matrix is 

Q SU(2) =&ag(~,-±q j ,q l ,q j ,q l \ . (El) 

The symmetry breaking pattern is assumed to be SU(4\2)l <S> SU(4\2)r <S> U(l)y — ► 
SU(4\2)v <8> U(l)v- The baryon field assignments are analogous to the case of SU(3) 
flavor. The nucleons are embedded as 

B ijk = 4= (€ijN k + dkNj ) , (F.2) 

where the indices i, j and k are restricted to 1 or 2 and the SU(2) nucleon doublet is 
defined as 

N = ( P j (F.3) 

The decuplet field 7^, which is totally symmetric, is normalized to contain the A- 
resonances T^ = %jk with i, j, k restricted to 1 or 2 and Tm = A ++ . The construction 
of the octet and decuplet baryons containing one sea or one ghost quark is analogous 
to the SU(3) flavor case [94] and will not be repeat here. 

The free Lagrangian for B and T is the one in Eq. ( 12.601 1 (with the parameters 
having different numerical values than the SU{2>) case). The connection to QCD is 
detailed in [94]. Similarly, the Lagrangian describing the interaction of the B and T 
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with the pseudo-Goldstone bosons is the one in Eq. ( 17.51 1 that can be matched to the 
familiar one in QCD (by restricting the B^ and %j k to the qqq sector), 



C = 2g A NS»A fl N + 5l 7V^JVtr(^) + g AN [T^A^N^m + h.c 



+2g AA T% ji S^T Uiljk + 2g x T% ji SrT Vtijk tr(Ai>) 



(F.4) 



where one finds at tree-level g\ = —2{D — F), g A = D + F, C = —g A N, and H = g AA , 
with gx = 0. The leading tree-level operators which contribute to A — ► JV7 have the 
same form as in Eq. H7.4K : of course the low-energy constants have different values. For 
transitions no additional tree-level operators involving supertrace of Q SU W appear. 

Evaluating the transition moments at NLO in the chiral expansion yields expres- 
sions identical in form to those in Eqs. d 7 . 6D — < F7T8l > with the SU(2) identifications made 
for C, H, D, and F. For the SU(2) coefficients in X PT one finds /3| = g A /V3 and 
0x = 5/(3y / 3) for the ■k ± . The corresponding values for the case of PQxPT appear in 
TableEH 



Table El: The 5(7(2) coefficients (i x and 0% in PQxPT for A -» iV 7 . 
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Appendix G 

More on Finite Lattice Spacing 
Corrections 

G.l 0(a) Corrections for Flavor SU{2) 

We consider the case of SU(2) flavor PQQCD 1 and summarize the changes needed 
to determine finite lattice spacing corrections to the electromagnetic properties of 
hadrons considered above. For the two flavor case, we keep the up and down valence 
quark masses non-degenerate and similarly for the sea-quarks. Thus the quark mass 
matrix reads 

m Q = diag(m u ,m d ,mj,mi,m u ,m d ), (G.l) 

while the SW matrix is 

c SU{2) = si ag (c v ,c v ,c s ,c s ,c v ,c v ). (G.2) 

Defining ghost and sea quark charges is constrained only by the restriction that 
QCD be recovered in the limit of appropriately degenerate quark masses. Thus the 
most general form of the charge matrix is 

qSU<2) =4i a g(l,-l, qj ,q l , qj ,q l \ , (G.3) 

which is not supertraceless. Analogous to the three flavor case, the vector-current will 
receive 0(a) corrections from the operators in Eq. ( I8.5D of which only the operator 0^ 
is relevant. The coefficient matrix associated with this operator is 

cf (2) = diag(c?,c?,c?,cf,c?,c?). (G.4) 

The 0(a) operators listed above in Sections I8.3H8.6I are the same for the SU(2) 
flavor group, however, the coefficients have different numerical values. Additionally 
there are operators involving str(Q + ^ '). These are listed for each electromagnetic 
observable below. 



J For brevity we refer to SU(4\2) PQQCD as 517(2). The distinction will always be clear. 
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Octet mesons 

In the meson sector, one has the additional term 

C = im 4 k x F^ v sir ( A+D^D v Yt + A + D^Y)D V Y) str(Qf (2) ). (G.5) 



Octet baryons 

In the octet baryon sector, there are terms which originate from A+ insertions 

C = ^[b 9 (BBA + )+b 10 (BA + B)]v^ u F^str(Q s + u{2) ) 

X. 






+-^ (BB) v„d u F^str(Q^^)str(A+) 

Ax 



+i [b' 9 (B[S^, S U ]BA + ) + 6' 10 {BIS,,, S V ]A+B)] F» v str(Qf (2) ) 

+i b' u (B[S^ S U ]B) F^ str(Qf (2) ) tir{A+), (G.6) 

and additional vector-current correction operators 

aCA <"f (nn\ „, a t?w ^wnSU(2)Ju(2)^ ^ taHAy miQ <-, 1kA ^u ai . v . {ri su(2)ju(2)^ 

(G.7) 



C = ^1 (BB) v^F^ str(Q su ^cf U{2) ) + ^1 (B[S^, S U ]B) F^str(Q su ^cT 



Decuplet baryons 

Next in the decuplet sector there are terms that result from A+ insertions 

C = ^ (T°A + %) v^F^ str(Qf {2) ) + ^ (T°%) v^F^ str(Qf (2) ) str(A+) 
A x A x 

+t 4 (T^A+%) F^ str(Qf (2) ) + i d' 6 (T^%) F^ str(Qf {2) ) str(„4 + ) 

+iL (r^A+T v A v a d»F ua str(Q S + U {2) ) 

A x V / 

+^- (^r">) v-d^F^ str(sf {2) ) str(„4 + ) 

(G.8) 
and also further vector-current correction operators 



3a c' 4 „ . 

' — o 

(T&T'A v a d fl F ua str(Q su ^cf U{2) ). (G.9) 



C = ^^ (T°%) v,d v F^ str(Q^( 2 ) C f (2) ) + Sm^ (T^) F"" str(Q^( 2 ) C f (2) 
A x 



A x V )" ""' 
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Baryon transitions 

Finally for the transitions, there are only new A+ insertions 

C = it 6 {BS^A+%) Fi» str(Qf (2) ) + A. (bS^A+T v A v a d fl F ua str(Qf (2) ) 

A x v ' 

+|| (BS^A+%) d a d»F\ str(Qf (2) ). (G.10) 

x 

For each electromagnetic observable considered above, contributions from all 0(a) 
operators in the effective theory are of higher order than the one-loop results in the 
chiral expansion. Thus one need only retain the finite lattice spacing corrections to the 
meson masses and use the previously found expressions for electromagnetic properties 
in SU(2) PQxPT in Appendices El IE! and|B as well as Refs. [94, 161]. 

G.2 Coarse-Lattice Power Counting 

Here we detail the 0(a) corrections to electromagnetic properties in an alternate 
power-counting scheme. We imagine a sufficiently coarse lattice, where ak x can be 
treated as 0(e), so that 2 

m q /A x , 

a 2 A 2 , . (G.ll) 

UVA X 

In this case, there are known additional 0(a 2 ) corrections [159] to the meson masses 
that are now at 0(e 2 ) and must be included in expressions for loop diagrams. The free 
Lagrangian for Bijk and T? k fields contains additional terms of 0(a 2 ) that correct the 
baryon masses, and modify the kinetic terms. Potential contributions due to the latter, 
whatever their form, must be canceled by wavefunction renormalization diagrams. 
The only contribution of 0(a 2 ) could come from tree-level electromagnetic terms but 
these are necessarily higher order. Thus in this power counting there are no unknown 
0(a 2 ) corrections for electromagnetic properties. 

The only possible corrections come from the 0(a) operators assembled above. A few 
of these do contribute at tree level and are spelled out below. 

Octet mesons 

The O(a) corrections to the meson form factors are now C(e 3 ) in the power counting. 
While the meson charge radii at NLO in the chiral expansion are at 0(e 2 ), further 
corrections in the chiral expansion are at C(e 4 ). Thus one can use the 0(a) operators 

z This power counting coupled with the chiral expansion is most efficient for valence Ginsparg- Wilson 
quarks where 0(a) corrections vanish. We thank Gautam Rupak for pointing this out. 



e 2 ~ < 
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to completely deduce the charge radii to C(e 3 ) [apart from C(e 3 ) corrections to the 
meson masses]. These O(a) operators are given in Eqs. ( 18.131 1 and ( 18. 141 1 and yield a 



correction 6 < r E > to the meson charge radii of the form 



5 <r 2 E >=Q 



24aA )i 
1^ 



[c v (2mi + m 2 ) + 3c s m 3 + c^a^g] 



(G.12) 



Notice that there are no corrections associated with an unimproved current operator 
in the sea sector since cf is absent. 

In the case of SU(2) flavor, there is an additional contribution from the operator 
in Eq. dG.51 ). At tree level, however, this operator vanishes. The only correction to 
Eq. AG. 121 ) in changing to SU(2) flavor is to replace 3c s with 2c s which reflects the 
change in the number of sea quarks. 



Octet baryons 

For the octet baryon electromagnetic properties, the 0(a) corrections to the charge 
radii are now C(e 3 ) and can be dropped as they are the same order as neglected I /Mb 
corrections. The magnetic moments, however, do receive corrections from local oper- 
ators. Specifically, the 0(a) operators which contribute to magnetic moments at 0(e) 
are insertion of A+ into the magnetic moment operator given in Eq. ( 18. 171 1 and O^ 
corrections given in Eq. H8.19H . Calculation of these corrections yields a shift in the 
magnetic moments 



Sfj, 



aM B lc v 



Alb' 1 + -b' 4 )-B(2b' 2 + b' 3 -b' 5 ) 



+ 3c s (-Ab' 6 -Bb 7 



+C(c s -c v )q jlr b' 8 + ^ 



-HA, a A - HA,pB 



(G.13) 



where qj\ r = qj + qi + q r . The coefficients A and B are listed for octet baryons in 
Table iGll while C = 1 for all octet magnetic moments and C = for the AS transition 
moment. Notice that there are no corrections associated with an unimproved current 
operator in the sea sector. 

In the case of SU(2) flavor, there are additional contributions given in Eqs. ( IG.6I > 
and dG.71 1. For the proton and neutron, we have 



5v su{2) 



aM B lc v 



A(h\+ l -b'^\ - B {2b' 2 + b' 3 -b' 5 ) + ± {b' 9 + b' w ) 



+ 2c s (±Ab' 6 -Bb' 7 + h' n ) + 



c s q jt + c v 



1 



Qjl 



b' 



2 



1 (I 

-AjJ,A,a ~ B fl A ,l3 + I r - Qjl J MA, 7 



+ -J Qjl VAa > , 



(G.14) 



where q jt = qj + q t . 
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Table G.l: The coefficients A and B for the octet baryons. 
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Decuplet baryons 

For the decuplet baryon electromagnetic properties in coarse-lattice power counting, 
the 0(a) corrections to the charge radii are 0(e 3 ) and the corrections to the electric 
quadrupole moments are 0(e), both of which are higher order than the one-loop re- 
sults. The magnetic moments, however, do receive corrections from local operators. 
Specifically, the 0(a) operators which contribute to magnetic moments at 0(e) are A + 
insertions into the magnetic moment operator given in Eq. A8.23I ) and 0f correction 
operators given in Eq. ( 18.251 ). Calculation of these corrections yields a shift in the 
magnetic moments 



6(i = 2aM B 



\c v Q (2d[ + 4) + c s Q 4 + (c s - c v )q ]lr d' A + c\Q ii A ,c 



(G.15) 



Notice that in 577(3) strQ = 0, hence there is no dependence on c\ in the above result. 
In the case of SU(2) flavor, there are additional contributions given in Eqs. dG.81 ) 
and < IGr.9b . The corrections to the A quartet magnetic moments are then 



Sfi su(2) = 2qMb J v ( 2Qd>i + Qd > 2 + d ,^ + s (Q d / g + d ^ + 



c b qji + c" 



Qjl 



d', 



+c\ [Qva,c + (1 - 3tyj)/4 )7 ] + SciQjlVA,' 



(G.16) 



Baryon transitions 

For the decuplet to octet electromagnetic transitions in coarse-lattice power counting, 
the 0(a) corrections to G^O) and 03(0) are 0(e) which are of higher order than the 
one-loop results. The G\(q 2 ) form factor does, however, receive corrections from local 
operators. Specifically these 0(a) operators which contribute to Gi(0) at 0(e) are the 
insertions of A+ into the magnetic dipole transition operator given in Eq. H8.29H and the 
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vector-current corrections given in Eq. H8.32H . Calculation of these corrections yields a 
shift of Gi(0) 

«5Gi(0) = aM B a T \J^ | c v U +t 2 + t 3 - ^4) + 3c s t 5 + c\ \ia?^\ \ , (G.17) 

where the transition coefficients ar appear in [6]. Again, at this order the result is 
independent of O (a) improvement to the electromagnetic current in the sea sector. In 
the case of SU(2) flavor, there is an additional dipole operator given in Eq. dG.101 ). 
At tree level, however, this operator vanishes. The only correction to Eq. dG.171 ) in 
changing to SU(2) flavor is to replace 3c s with 2c s which reflects the change in the 
number of sea quarks. 
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